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Stall in a Single-Stage Axial Compressor 
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SUMMARY 


The investigation was conducted on a single-stage axial flow 
compressor consisting of a row of 30 inlet guide vanes with solid 
body prerotation followed by a row of 31 rotor blades with con 
stant work along radius, at design, and a row of 34 stator blades 
re-establishing conditions at rotor inlet. The stage has a hub to 
tip diameter ratio of 0.5 with a tip diameter of 0.50 m. and is set 
up between an entrance cone and a diffuser with a throttling de- 
vice at the downstream end. The rotating stall characteristics 
were carefully investigated by pressure measurements performed 
with stagnation and static pressure probes connected to quick 
response pressure gages, as well as by smoke observations. 

The stall occurs toward the tip of rotor blades with closing 
throttle at the peak pressure point. It appears initially in the 
form of a single propagating stall region, which interests three 
blade channels and which rotates at a speed equal to 0.68 the 
rotor speed. 

The number of stalled regions increases upon further decrease 
in flow rate, from one to eight regions. Further throttling causes 
separation at the suction side of the guide vanes; a one-region 
total span rotating stall is then observed on the guide vanes, 
which affects the instantaneous flow rate through the stage. 

The propagation speed of the stall, and to a certain extent, the 
width of the stalled regions, appear to be independent of the 
number of regions for the first kind of stall which originates at the 
rotor blade tips. However, the intensity of stall of the different 
regions is largely and erratically variable from one to another 
cycle of rotating stall, especially for the multi-region pattern 

Stroboscopic observation of smoke filaments shows a_ cyclic 
variation of the air velocity direction at the rotor inlet, with a 
frequency equal to the number of turns per second of rotation 
stall times the number of regions. The variation is governed 
by separation on the rotor blades and increases gradually in 
amplitude with the number of stalled regions 

High-frequency pictures confirm the observations made with 
pressure gages 
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SYMBOLS 


= vector absolute velocity 
= vector relative velocity 
tangential component of absolute velocity 
axial component of absolute velocity 
= tip velocity 
circumferential velocity at radius 7 
= dimensionless axial velocity at radius r, w/l 
average flow coefficient 
dimensionless tangential velocity at radius r 
radius at the tip of the blades 
= dimensionless radius r/r 
airfoil chord length 
cascade pitch 
maximum thickness of cylindrical sections of blades 
direction of absolute velocity from compressor axis 
direction of relative velocity from compressor axis 
local total pressure rise coefficient 
mass weighted average total pressure rise coefficient 
local static pressure rise coefficient 
mass weighted average static pressure rise coefficient 
compression efficiency 
r.p.m, 
absolute rotative speed of rotating stall 
lift coefficient per unit length of blade 
drag coefficient per unit length of blade 


(1) INTRODUCTION 


OTATING STALL in an axial compressor has been 
R the object of extensive experimental investiga 
tionst during the course of the last few years. Its 
characteristics have been described in many instances 
and its mechanism seems to be clearly understood. So 
far, two main methods of approach have been used 
the hot wire anemometer technique and static pressure 
measurements. In order to complete the observation 
already described, an attempt has been made to apply 
two other methods of approach for investigating the 
rotating stall on a single stage axial flow compressor 
stagnation pressure measurements and smoke surveys. 


+ See References. 
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Fic. 1. Schematic diagram of compressor installation and axial 
measuring stations along compressor. 


Both methods complete gach other and, together with 
static pressure measurements, seem to permit a pre- 
cise determination of rotating stall characteristics. 

Smoke observations show the development of stalling 
on the rotor blades quite clearly; they confirm the as- 
sumptions already made for the mechanism of the 
onset of rotating stall. The rotor speed being kept con- 
stant, partial-span stall is first observed at the tip of 
rotor blades, with closing throttle from peak pressure 
point. The number of stalled regions increases upon 
further decrease in flow rate, from one to eight regions, 
while no significant variation is observed in the speed 
of stall. Further throttling causes separation on the 
suction side of the guide vanes, the severity of which 
increases towards the tip. <A one-region total-span 
rotating stall is then observed on the guide vanes, which 
affects the flow rate through the stage. The absolute 
speed of this stall is about three times lower than the 
speed of partial stall on the rotor blades. Further de- 
creases in flow rate cause the disappearance of periodic 
phenomena; random stall is observed until, at a very 
low flow rate, an intense one-region root to tip stall is 
initiated on the rotor blades. Reverse flow is then ob- 
served in the stalled region, which extends through the 
entire stage, up to a certain distance upstream of the 
guide vanes. 

Thence three distinct types of stall are observed; 
multiregions are observed only in the first type, and the 
stall speed seems in this case not to depend on the num- 
ber of regions. As to the instantaneous flow rate at a 


TABLE | 

Distributor (30 blades) 

é 0.50 0.625 0.75 O.875 l 
s/¢ 0.523 0.653 0.783 0.913 1.047 
¢ (mim.) 50 dO 50 50 50 
t/ce (max.) (0.124 0.112 0.0998 0.087 0.0734 
ay? (outlet) 14°39’ 18°24’ 22°13’ 26°09’ 30°25’ 
Cr-(c/s) 0.519 0.656 0.800 0.954 1.126 
Rotor (31 blades) 

s/¢ 0.551 0.703 0.861 1.028 1.207 
c (mim.) Hf) 15 14 15 42 
t/c (max.) 0.150 0.130 0.110 0.090 0.070 
31° (inlet) 31°28’ 39°38’ 47°31’ 5°09" 63 
Bo° (outlet) 24°08’ S710’ 13% «6333 °a9" 18°04’ 
Cr:(c/s) 2.114 1.730 1.404 1.140 0.927 
Stator (34 blade S) 

s /¢ 0.511 0.641 0.768 O.895 1.027 
¢ (mm.) 15 15 15 1 15 
t/c (max.) 0.140 0.123 0.105 O LOSS 0.070 
a;° (outlet) 13°28 16°39" 31°30’ 37°33’ 35°20’ 
ae° (inlet) S202 80°40" S1°07" 53°11" 57°30’ 
Ci:(c/s) 1.736 1.488 1.316 1.162 1.081 
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given point, it remains constant at a given throttle for 
the first type of rotating stall, while it fluctuates period 
ically for the second and third types, although in these 
latter cases the mean flow rate through the compressor 
remains constant. 

Those results may hold only for the present particu- 
lar design, but it must be stressed that this design simu- 
lates internal flow commonly employed in modern high- 
performance turbomachines. 

An attempt has been made of the application of 
existing theories for the prediction of stall propagation 
speed of the rotating stall of the first type. In this case 
the stall interests only a narrow field about the rotor 
row, an assumption which is made in most of the cal 
culations. Sears theory’ yields results in agreement 
with the experiments, so long a suitable lag is intro- 
duced, of which evidence seems to be given by the 
smoke pictures (see, for instance, Plates 1 and 2). 
Some details of the application of the Fabri and Sies 
trunck theory,'' as suggested by Prof. M. Roy, are 


given in the Appendix. 


(2) Test COMPRESSOR 


The stage consists of a row of 30 inlet guide vanes with 
solid body prerotation, followed by a row of 31 rotor 
blades with constant work along radius at design, and 
a row of 34 stator blades re-establishing flow conditions 
at rotor inlet. Design flow angles at the rows from the 
axial direction, as well as the main characteristics of the 
blades, are given in Table 1. The stage has a hub to 
tip diameter ratio of 0.5, with a tip diameter of 0.50m. 
The blading was computed, with the assumption of per- 
fect fluid, by the use of the linearized three-dimen- 
sional flow theory (axisymmetric flow) to give a total 
pressure rise coefficient of 0.637 at a flow ratio of O.4S. 

The cylindrical sections of the blades at different 
radii were computed as cascades of airfoils, by means 
of the conformal transformation of suitable ovals. The 
design allowed sufficient space between the rows for the 
instrumentation required to permit a detailed analysis 
of blade row flow distribution; the cylindrical casing 
was made of transparent perspex to permit visual 
observations. A cross-sectional view of the test in- 
stallation is presented in Fig. 1, where the principal 
dimensions are summarized. 

The air flows into the compressor through a double 
curvature entrance cone. The shape of this cone as 
well as the streamlined shape of the front part of the 
hub have been designed so as to provide a uniform 
flow across the compressor inlet, with a very thin bound- 
ary layer at the walls. The contraction ratio is about 
6.4. <A fine gauze was installed at the entrance cone 
inlet in order to damp out eventual gusts. Two short 
ducts of constant annulus area are seen between the 
compressor stage itself, on one hand, and between the 
compressor stage and exit cone, on the other hand. 
Those ducts are of sufficient length for allowing the flow 
to assume very nearly its asymptotic pattern and per 
mit easy measurements. An interconnected set of four 
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2 (smooth flow) 


5 (stel1) 


6 (flow reestablished) 


PLATE 2. Five region rotating stall (high-frequency pictures), ¢ = 0.37, — = U.88, 2,460 r.p.m., time interval between successsive pic 
tures = 1/450 sec. 
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wall static pressure orifices is located midway through 
the front cylindrical duct, which was calibrated to give 
the total flow coefficient. The exit portion of the in- 
stallation consists of a slightly diverging duct with an 
expansion area ratio of 1.36. The throttling device at 
the downstream end of the exit cone is driven by an 
electric motor operated from the test bed stand. 

In order to permit easy visual observations, the blades 
of the three rows are set on the hubs. The pieces of 
metal required for holding the instrument carriages 
screwed on the outside wall of the casing make only a 
small obstruction. A driving mechanism allows the 
rotation of the casing; the clearance between the tips 
and the casing is 0.2 to 0.6 mm. A 22 hp. electric 
asynchronous motor drives the compressor shaft by 
means of pulleys and belts. 

Most of the measurements were made at a rotative 
speed of 2,460 r.p.m.; hence, compressibility effects 
were negligible, and the blades could be considered rigid 
as far as aerodynamic forces were concerned. The 
power input to the compressor could be found from ac- 
curate measurements of the rotative speed and of the 
driving torque which was determined by means of a 


torquemeter. 


(3) INSTRUMENTATION 


For the evaluation of overall performance and individ- 
ual blade row flow details, in the normal range of op- 
eration as well as inside the rotating stall range, exten- 
sive instrumentation was located at stations 0, 1, 2, 
and 3 as shown in Fig. 1. Stations | and 2 are situated 
midway between adjacent rows; station 0 and 3 are 
situated midway through the front and rear ducts of 
constant annulus area. Static pressure orifices were 
provided also at the outer wall immediately upstream 
and downstream of the rotor blade tips. 


Flow Rate Measurements 


Station 0 was used to determine the radial distribu- 
tion of velocity at the compressor inlet and conse- 
quently, by integration, the flow rate as a function of the 
rotational speed and of the throttle position. To this 
effect a static probe of the conventional type was 
shifted radially and, for steady flow, readings were taken 
on an alcohol manometer, of the pressure difference be- 
tween the pressure transmitted by the probe and the 
stagnation pressure at the hub nose. This last pressure 
was taken as reference pressure in all pressure measure- 
ments. The set of static holes at station 0 were cali- 
brated after the static probe surveys. During opera- 
tion with partial-span stall the net flow rate through 
the compressor and the torque necessary to drive the 
compressor remain essentially constant. 

To permit the correct measurements of unsteady flow, 
(so long there is no return flow), the static probe was 
then connected to a dynamic pressure gage, as de- 
scribed below. Station 0 was also used for the distri- 
bution of smoke through streamlined tubes of con- 
venient diameters in order to get a picture of the flow 


pattern through the stage, and especially upstream 
and downstream of the rotor row. 

When return flow is observed at station 0, at very low 
flow rates, as will be seen further on, the rate of flow 
had to be measured by means of static pressure holes 
drilled on two flanges at the exit of the throttling device 
(those flanges are not shown in Fig. 1). 


Evaluation of Flow Details at Station 1, 2, and 3 in the 

Normal Range of Operation 

Station 1 was used to evaluate both the flow condi- 
tions downstream of rotor and upstream the stator. 
Station 3 was used for the evaluation of flow conditions 
downstream of stator. 

Probe carriages at those three stations would permit 
radial surveys with different types of probes—one yaw 
and stagnation pressure probe of the well known claw 
type, and one static probe. Circumferential surveys 
were achieved by the rotation of the casing. The 
probe carriage could also hold minute streamlined tubes 
for smoke distribution, for carrying out detailed ob- 
servations of the flow pattern through individual row 


Investigation of Flow Details in the Rotating Stall Range 


The same stations as described above were used for 
the performance of pressure and smoke radial and cir 
cumferential surveys, but the instruments used for this 
part of the program are different from the ones used 
for investigations in the normal range of operation. 
Moreover, two more stations, a and b, already de- 
scribed, immediately upstream and downstream the 
rotor, were used for the measurement of static pressure 
by means of small orifices drilled through the casing. 

The stagnation and static pressure probes are made 
of monel tubes of short length (about 60 mm.) and small 
inner diameter (1 mm.). The stagnation pressure 
heads are made as conventional pitot tubes, whereas 
the static pressure heads are of the disc type. Very 
short rubber tubing is used for connecting the probes 
with electrical pressure gages. 


Electrical Pressure Gages 


The electrical pressure gages are of the variable air gap 
inductance type. They answer to the following re- 
quirements: small size, good frequency response to 
several hundred cycles per second, insensitivity to 
mechanical vibrations, sensitivity at low pressures, 
convenient output, simple construction. 

The gage components are arranged in the following 
way: a corrugated thin diaphragm of magnetic ma- 
terial (u-metal), is clamped between two cup-shaped 
case halves (made also of u-metal); each of them has a 
cylindrical core piece through which is drilled a pressure 
opening. An inductance coil is cemented about each 
core piece. Each case half, together with the dia- 
phragm, constitutes, for the coil inside, a magnetic cir- 
cuit with a small air gap (1 10 mm.) between the core 
piece and the diaphragm. 

The pressure openings permit the connection of 
both the air gaps to two pressure probes or one to a 
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Stagnation and static pressure taps in connection with 
pressure gages. 


Fic, 2 


pressure probe on one side and the other to a pressure 
orifice at the hub nose (see Fig. 2). The coils are con- 
nected in adjacent arms of an a.c. Wheatstone bridge, 
which is fed by an audio frequency oscillator. The 
block diagram of the carrier equipment used is shown 
in Fig. 4; the electronic equipment itself is shown in 
Fig. 3. Galvanometers of the moving magnet type, 
of suitable natural frequency and damping, are used in 
the recorder. The static and dynamic behaviors of the 
pressure gages associated to different pressure heads 
have been thoroughly investigated. 

Frequency response curves to direct repeated rec- 
tangular impulses were established with a pressure pul- 
sator, simulating the flow conditions about the dif- 
ferent heads, which are prevailing at the different sta- 
tions when the compressor is stalling. 

A typical static calibration curve is shown in Fig. 7, 
as well as a frequency response calibration curve. No 
acoustical resonances were observed up to frequencies 
of the order of 250 cycles per sec. Satisfactory be- 
havior of the pressure gages was consequently to be ex- 
pected for the investigation of rotating stall of which 
maximum apparent frequency is around 200 cycles per 
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Fic. 3. Electronic equipment, stroborama, and recorder 
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sec., When the rotational speed of the rotor is 2,460 


r.p.m. 


Smoke Generator and Light Sources 


Ammonium chloride was used in the present investi- 
gation, as in the previous ones which the author had 
performed on propellers. This was produced with a 
generator described previously (see reference 12, pp. 
15-48). 

The smoke is distributed in the stream, most usually 
at station 0, with streamlined tubes made of monel 
which the smoke cannot corrode (the blades themselves 
made of duralumin are protected by a thin coat oi 
alumin). 
capacitor of suitable capacity in a tube containing rare 
The high voltage 


Flashes are produced by discharge of a 


gases and a few drops of mercury. 
generator which has been used is the so called ‘‘stro- 
borama™’ built by Seguin in Paris, in which a thyratron 
switch has been substituted for the mechanical one 
This device can be used as a stroboscope with a flash 
frequency ranging between € to 50 per sec. Alterna 
tively, series of 10 to 12 flashes of very short duration 
can be produced, the time interval between each flash 
varying between 1/300 to 1/1000 sec. In order to 
achieve the highest contrast for the observation of 
smoke filaments (the filament lines are as described by 
Bairstow), the light source must be located opposite to 
the observer. Also the use of a concave mirror in 
order to focus the light over the smoke jet is recom- 


mended. 
(4) RESULTS AND DISCUSSION 


Performance Characteristic of the Compressor in the 

Steady State, and in the Rotating Stall Range 

Mass weighted average total and static pressure rise 
coefficients for the rotor are plotted in Fig. 5 against 
flow coefficient; they are based on the difference in pres- 
sure between stations 2 and 1. 

Maximum pressure rise coefficients are seen to take 
place at @ = 0.388. From this point, with further 
throttling WV, seems to be decreasing continuously. 





Fic. 4. Block diagram of the carrier equipment used with the 
pressure gages. 
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The efficiency curve, plotted also on the same Figure 
against ¢, appears to be very flat in a wide range of flow 
coefficient, up to the peak pressure point. 

Rotor cascades radial variation of C, and Cp are 
plotted in Fig. 6 for design ¢ and @ = 0.3SS correspond- 
ing to maximum pressure rise. There is a tendency 
for C, to increase towards the tips of the blades when 
the compressor is throttled. Simultaneously the Cp 
increases rapidly, forecasting the onset of stall, for sec- 
tions toward the tips of the blades. 


Rotating Stall Characteristics 


The characteristics described hereafter are depicted 
after the examination of dynamic pressure records and 
from smoke observations performed as stated in Sec- 
tion (3) above. Those will be discussed below. 

The stall occurs towards the tip of rotor blades with 
closing throttle at the peak pressure point. It appears 
initially in the form of a single propagating region. 
This type of stall starts at @¢ = 0.384 and extends from 
the outer radius to € = 0.75, while toward the inner 
radius the fluctuations in pressure become much smaller 
and irregular. It interests three blade channels at the 
outer radius and contracts inwards. 

Further throttling results in the formation of two 


stalled regions at @ = 0.38, then of three, four, five, 
six, and eventually, seven or even eight regions, around 
@ = 0.35. 


The absolute rotating speed of the stall has been 
measured for all patterns by comparing records made 
with two stagnation pressure heads at station 2, 45 
apart along the circumference. No significant varia- 
tions are observed following the pattern; the rotating 
speed seems to be practically irrespective of the number 
of stalled regions and equal to 0.68 the rotating speed 
of the rotor. This fact is confirmed by stroboscopic ob- 
servation of smoke as will be seen further on. Also the 
width of the regions and their radial extent seems not to 
depend on the number of regions. 

However, the intensity of stall of the different regions 
appears to be largely and erratically variable from one 
to another cycle of rotating stall, especially when the 
number of regions is high (6 to 8); for instance, one or 
two regions may vanish and temporarily modify the 
pattern. Still, on the whole, the severity of stall 
seems to increase with the number of regions. 

As a rule, for the two- and three-region pattern, the 
distribution of the regions is not regular along the cycle. 
But apart from the erratic variations pointed out 
above, the distribution of regions tends to become 
symmetrical over the cycle, when the number of regions 
is above three. 

Simultaneously with the increase in number of re- 
gions, under ¢ = 0.35, the stalled regions extend up- 
stream of the rotor row. The interaction upon the 
flow round the guide vanes is large enough from ¢ = 
0.35 downwards, to cause separation on the suction side 
of these blades, from the leading edge. A one-region 
root to tip rotating stall pattern is then initiated along 
the guide vanes with increasing severity towards the 
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tip, which affects the flow through the entire stage. It 
interests circumferentially half of the blades of the 
distributor row. 

The absolute speed of the stall in terms of the speed 
of the rotor is 0.293. Therefore it is quite distinct from 
the previous speed of stall on the rotor, but it must be 
pointed out that this new speed corresponds to a quite 
different type of stall. 

No change in the characteristics of this new pattern, 
apart from the severity of stall, is observed with clos 
ing throttle, down to @ = 0.245. Further decreases of 
the flow rate cause the disappearance of the periodic 
phenomena. Random stall is then observed, with clos 
ing throttle, to be especially intense over the outer 
parts of the blades of the three rows, until, when @ = 
0.146, a sudden change of the pattern occurs—a one 
region root to tip stall is observed which is initiated on 
the rotor blades and interests half of the blades of the 
rotor row. In the stalled region, separation takes 
place from the leading edge of rotor blades and a re- 
verse flow occurs through the entire stage, which ex- 
tends up to a certain distance upstream of the guide 
vanes. To this new pattern of stall corresponds first 
a net increase of the mean flow rate, @ jumping from 
0.145 to 0.213. Further throttling brings no change 
of this ultimate type of rotating stall, down to ¢ = 
0.186. Inside the range 0.213 ¢ 2 0.186, the abso- 
lute speed of stall remains constant and equal to 0.207, 
the rotational speed of the rotor. 

To summarize, three distinct types of rotating stall 
are observed with closing throttle from peak pressure 
point. The first kind initiates on the rotor blades and 
affects only the outer parts of the rotor blades. The 
stalled region extend steadily upstream upon further 
throttling, until the second type of stall occurs as the 
result of separation over the guide vanes. The ulti- 
mate type of rotating stall is initiated on the rotor 
blades as the result of total span separation over those 
blades, from the leading edge. To the first type corre 
sponds, for a given throttle and at a given point, no 
variation against time of the rate of flow through the 
stage. To the second and third types correspond, for 
a given throttle and at a given point, periodic fluctua- 
tions of the instantaneous rate of flow. Still in both 
those latter cases the mean flow rate through the 
stage remains constant. To each type of stall corre- 
sponds a given propagation speed, which is practically 
irrespective of the number of regions, at least for the 
first type. For the second and third type, a one-re- 
gion Stall only is observed. 


Stagnation Pressure Fluctuations Due to Rotating Stall 


An investigation has first been made of the influence 
of the angle of the pressure head, from the rotation axis, 
on the galvanometer readings. The head was held at 
station 2 and at = 0.964 and was set at 55° from the 
rotation axis, which orientation was found to be cor- 
rect by means of the yaw probe, in the steady state, at 
peak pressure rise (@ = 0.388). Tests were then per 
formed after further throttling, when the one or the four 
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ROTATING STALL IN SINGLE 


region stall was established. Decreasing or increas- 
ing the setting angle would result, after the readings, as 
shown in Fig. 8, in an apparent reduction of the maxi- 
mum pressure rise during the periods when the flow is 
not separated. Consequently, the readings are as- 
sumed to be correct in the unstalied zones, when, at any 
radius, the head is set at an angle as determined by pre- 
vious measurements in the steady state. It was ad- 
mitted that such settings would be still convenient for 
measuring to a certain extent the pressure drops in the 
stall zones. At least the readings in the stall zones 
should permit an appreciation of the severity of the 
stall. 

An example of stagnation pressure records at station 
1 and 2 is given in Fig. 9, where are shown the time 
base (1/100 sec.); the fluctuations in pressure rise at 
station 2, with the stagnation pressure at the hub nose 
as reference pressure; the fluctuations in pressure rise 
at station 1, with the same reference pressure; and 
finally, the rate of flow as measured by means of a 
probe at a given circumferential position at station 0 
(first and second types of stall). 

The readings concerning the pressure fluctuations at 
station 1 are of small amplitude; they correspond, 
with decreasing flow, to a slight rise followed by a drop. 
Fluctuations in the direction of the velocity which are 
quite significant in this field, as will be seen further on, 
may have some bearings on the readings. 

For the first type of rotating stall and for one-region 
stall the pressure drops at station 2 appear to be, on the 
average, equal to about one third of the pressure rise in 
the unstalled zones. The pressure rise coefficient out- 
side the stall zones remains equal to the peak pressure 
rise in the steady state. 

An example is given of the pressure records in the 
case of multiregion stall for eight-region stall, which 
occurs for ¢@ from 0.332 to 0.342. Radial variation of 
relative stagnation pressure rise fluctuations and of 
inaximum pressure rise coefficient at stations 2 and 3 and 
at @ = 0.384 (one region stall), are shown in Fig. 11. 


Static Pressure Fluctuations for the First Type of Stall 


Measurements at the Outer Wall.—The pressure are 
measured with the stagnation pressure at the hub nose 
as reference pressure. At station 2 a drop in pressure is 
observed in the stall zones, while, outside the stall zones 
and for patterns ranging from one- to four-region stall, 
the rise in pressure remains equal to the rise in the 
steady state at ¢ = 0.388. This pressure rise decreases 
slightly with further throttling. At station | a pres- 
sure drop is observed outside the stall zone, which, when 
the number of regions is one or two, is equal to the pres- 
sure drop in the steady state at @¢ = 0.588. In the 
stalled zones, with decreasing flow the fluctuation in 
pressure is composed of a rise in pressure followed by a 
drop. For one-region stall the fluctuations in pressure 
at station 2 are smaller than half the pressure rise out- 
side the stall zone. At station | the increase in pres- 
sure corresponding to the fluctuations, and for the same 
stall pattern, is about one third of the pressure drop 
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fluctuations and of maximum pressure rise coefficient at stations 
2 and 38, 2,460 r._p.m., @ = 0.384 (one-region stall 


outside the stall zone; the decrease in pressure is about 
one quarter of the same pressure drop. For the same 
pattern and outside the stall zone, at station 2 the pres- 
sure rise is nearly equal to the pressure drop at station 1. 

Similar phenomena are observed at stations a and b 
adjacent to the rotor on both sides. However, outside 
the stall zone the pressure differences are smaller 
compared with the differences in pressure measured at 
Inside the stall zones, on the con 
A large part of 


stations 1 and 2. 
trary, the fluctuations are far larger. 
this result seems to be due to a clearance effect. 

An example of the static pressure records at stations 
1 and 2 is shown in Fig. 10. The composition of this 
diagram is similar to that of Fig. 9. In the case of 
multiregion stall, the records concern the six-region 
stall pattern which occurs for @ between 0.356 and 
0.3608. 

Measurements Along the Radius—-The observations 


are very much similar to the ones above described. 


Smoke Surveys and Smoke Pictures 


Photographs have been taken of the filament lines in 
the steady state and in the rotating sta!l state. Those 
were taken with the apparatus already described, by 
means of an ordinary camera for stroboscopic obser- 
vations, and by means of a rotating drum for images in 
rapid succession. A sample of the images in the steady 
state is shown in Fig. 12, corresponding very nearly to 
peak pressure rise. 

All along the radius, vortices are seen shedding from 
the rotor blades, as the circulation in the flow increases 
from hub to tip. Although this is not apparent in the 
picture, those vortices are still present over the stator 
blades and even downstream of the stator. 

Stroboscopic observations of smoke filaments be- 
tween the distributor row and the rotor row and be- 
tween the mean radius and the outer radius, performed 
at a frequency slightly smaller than the absolute rota- 
tional speed of the stall (first kind of rotating stall), 
show up for one-region stall—a nonuniform rotation 
of the filament lines in the sense of the stall rotation. 

No significant variation of the rotational stall speed 
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‘lament lines downstream of rotor, ; 
O.40, 2,460 r.p.m 


is observed with the number of regions. But the am 
plitude of the fluctuations of the filament lines increases 
with the number of regions, to become very large at 
flow rates corresponding to the second type of rotating 
stall. The frequency of fluctuations of the filament 
lines is equal to the rotational speed of propagation of 
the stall multiphed by the number of stalled regions. 
lo these fluctuations correspond obviously a cyclic 
variation of the relative incidence at the rotor blades, 
which is far from being harmonic. 

Plates 1 and 2 are extracted from high-speed films 
taken by means of the rotating drum. In both plates 
the number of frames per second ranges from +50 to 
DOG. 

Plate 1 concerns the second type of stall. The flow 
The fol 


lowing frame, which is not reproduced here, shows a 


is seen to stall severely in the images 5 and 7. 


reversal of the flow through the rotor; as a matter of 
fact, the air appears to be eyeling in front of the rotor 
in the stalled regions. The complete stalling cycle 
appears to last about 2 100 sec. Such a time seems to 
be rather long, but it must be recalled that when the 
compressor 1s stalling severely, the number of stalled 
regions is quite irregular (the effects of the first kind of 
stall appear to be superimposed on the fluctuations of 
the flow rate due to the second kind of rotating stall). 

Plate 2 concerns the five-region rotating stall state 
(first type of stall). Stalling is shown clearly in images 
tf and 5, but comparison with the images of Plate | 
shows the stall to be less severe. The complete cycle is 
shown to last about 1 130 sec.—-a duration which is in 
good agreement with the data of dynamic pressure 
ieasurements. 

In both series of pictures there is some evidence of 


aerodynamic lag. 


(5) CONCLUSIONS 


Investigations have been made of the rotating stall 
flow pattern on a single-stage axial compressor by 
means of dynamic stagnation and static pressure meas- 
urements and of smoke surveys. Both techniques 
vield results which are in good agreement. Smoke sur- 
veys appear to be of great assistance in the explanation 
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of the rotating stall mechanism. Three distinct type 
of rotating stall can occur, but application of the exist 
ing theories should be sought for the first type only 


APPENDIX 


Principle of the Theoretical Computation of Rotating Stall 
Absolute Rotative Speed, After Fabri and Siestrunck 
The calculations of references 10 and 11, which art 

performed in the case of a single row, can be extended 

to the case of a rotor row between two stationary rows 

The following assumptions are made which permit 
calculations to be performed the blade rows ar 
made up of an infinite number of blades, outside the 
stall zones the rotor cascade relative inlet flow angk 
remains equal to the zero lift angle, there is no through 
velocity inside the stalled zones, and the stall zones 
may be represented as wakes in the Helmoltz sense. 

The disturbances at the rotor do not interact notice 
ably on the distributor at the onset of the first kind of 
rotating stall, (an assumption which is in agreement 
with the observations as described above in the present 
paper). 

An expression for the ratio NV is derived, from the 
above assumptions and from the compatibility condi 
tions at the wakes boundaries throughout the rotor 
cascade, as 


) 


¢g tan 6 


hN 1/2 {I 


._— 


with (1 2) pu 


.— 
> 
| 
re 


being the local static pressure rise coeflicient before 
stall, p) the stagnation pressure at the hub nose, ¢ the 
local flow coefficient before stall, and 8, the rotor cas 
cade relative inlet flow angle. 

The above formula applied to the case of the com 
pressor of the present paper, at dimensionless radius 
0.8, yields h N = 0.72, instead of 0.68, as derived 
from measurements for the propagation speed of the 
rotating stall of the first kind. 
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Aerodvnamics of Nonuniform Flows as Related 
to an Airfoil Extending Through a Circular 
Jet 


SCOTT RETHORST** 
California Institute of Technology 


SUMMARY 


A method has been derived for the determination of the down 
ish in the field of an airfoil extending through a circular jet 
Phis solution for values of downwash was then employed in lifting 
surface theory to determine the lift distribution on a wing ex- 
tending through a circular jet 

Phe basic approach for the downwash solution was to con 
sider the flow resulting from the viscous action of the vortex sheet 
constituting the jet boundary The method of solution is essen 
tially based on a division of the flow thus induced by this jet 
boundary into parts which are even and odd with respect to the 
direction of flow. The analysis due to the even part alone is 
similar to previous theories, which in effect disregarded the odd 
part. Such previous results for the lift distribution based on the 
even part alone differ considerably from the experimental values 
rhe results based on the total of the even and odd parts show 
good agreement with the experiments, illustrating the necessity 
f including the odd part of the flow when the segment of the 
wing immersed in the jet is of low aspect ratio 

Phe problem has been solved in parametric form, so that the 
results may be employed to determine the characteristics of any 


geometry wing-jet combination at any jet velocity ratio 


SYMBOLS 


= jet radius 


x, ¥,5 = rectangular coordinates (Fig. 1 
é,7, ¢ = nondimensionalized rectangular coordinates, x/? 
¥/Tor2 
tr, 7,5 = rectangular coordinate running variables 
a, 8, y = nondimensionalized rectangular coordinate running 
variables, £/ro, ¥/To, 3/1 
Kt... = cevlindrical coordinates 
p = nondimensionalized radius, 7/7 
= semiwidth of horseshoe vortex 
= wing span 
S = wingarea 
A = vector velocity potential of flow induced by horsc- 
shoe vortex 
@ = sealar velocity potential of flow induced by jet 


boundary 
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' jet velocity 
I outer or free-stream velocity 
Mu = V./V 
I = circulation strength 
v = vector perturbation velocity 
(,v,w = reetangular perturbation velocities 
( = radial perturbation velocity 
ILA = modified Bessel functions 
i, 28 = coefficients in series representation of @ 
= summation index 
\ = variable of integration 
,d = spanwise locations of horseshoe vortices, nondimen 
sionalized with respect to jet radius 
i = perpendicular distance to a line vortex 
R = three-dimensional nondimensionalized potential dis 
tance = ¥ oc? + (n 8)? + & 
R = two-dimensional nondimensionalized potential dis 
tance = \ n 8B)? + ¢ 
W = number designating a particular horseshoe vortex, 


starting from left wing tip 

= number designating a particular downwash point, 
starting from left wing tip 

dl = differential length element of horseshoe vortex, non 


dimensionalized with respect to jet radius } 


‘ = downwash points 

Q = vortex points 

N = number of horseshoe vortices across total wing span 
F(z, 7) = downwash velocity at any point (¥, 7) caused by a 


rectangular horseshoe vortex of unit semiwidth 
ind circulation strength equal to 4z 
F, = wing downwash coefficient; the downwash at any 
downwash point P?, due to the mth horseshoe 
vortex 
G(é, 7) = downwash velocity at any point (£, 7) caused by the 
jet boundary reaction to a rectangular horseshoe 
vortex of unit semiwidth and circulation strength 
equal to4z 
Gy, = jet boundary downwash coefficient; the downwash 
caused by the jet boundary at any point ?, due to 
the mth horseshoe vortex 
= cvordinates of a point on the wing surface with re 
spect to root quarter chord, based on semiwidth 
of horseshoe vortex 
] = coordinates of a point in the wing surface with 
respect to the origin of a unit horseshoe vortex, 


based on semiwidth of horseshoe vortex 


W = Wronskian of Bessel functions = J(A)A‘(A 
I'(X)K(A) = 1/A 
IT = modified Wronskian VOOR VP) K'O 
Note lim JW= +(1/A)I 
1 >I 


INTRODUCTION 


HE AERODYNAMICS of an airfoil in a uniform flow 
I field is a highly developed subject. If the flow 
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Fic. 1. System of coordinates 


field is nonuniform, the subject becomes much more 
complex. <A general attack on this nonuniform flow 
problem was made by von Karman and Tsien in 1945.* 
The method developed was subsequently utilized by 
Fejer for the case of a flow whose velocity varied 
linearly in the spanwise direction.‘ 

Many examples of nonuniform flow arise in fluid 
mechanics, which would benefit from recognition and 
treatment as special cases of this branch of aerodynam- 
ics. Among these are vectored slipstream, jet flap, 
Coanda effect, supercirculation, thrust augmentation, 
and shrouded propellers. 

In general, the areas of endeavor which have stimu- 
lated interest in the aerodynamics of nonuniform flow 
fall into one of two categories. First, there are the 
many arrangements such as those listed above which 
involve the behavior of airfoils immersed in or adjacent 
to flows of finite extent generated by and external to 
power plants. Second, such flow conditions are en- 
countered within certain power plants, providing, for 
example, the environment of nonuniform flow in which 
the blades of turbomachinery operate. 

This paper is concerned primarily with the first of 
these areas of endeavor. Since the power plants 
which are used to generate these flows of finite extent 
most often involve rotating machinery, this paper 
specifically treats a jet which is circular. 

The problem of determining the characteristics of an 
airfoil extending through a circular jet has itself been 
investigated from several points of view. One is the 
influence of a propeller slipstream on an airplane's 
flight characteristics. Another is the influence of the 
boundaries of a wind tunnel on an airfoil being tested 
therein. 

An analytical treatment of the problem of a wing in a 
propeller slipstream was given by Koning,’ based on 
lifting line theory. The results of Koning’s analysis, 
because of his method of computation, are valid only 
for small increments of the jet velocity. Other early 
estimates®* were based on semiempirical factors, and 
gave satisfactory results within the speed-power range 
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for which they were developed. However, the higher 
ratios of jet velocities attainable through turboprop 
engines have stimulated renewed interest in solutions 
valid at all jet velocities. This interest is further aug 
mented by the development of certain types of vertical 
take-off aircraft, and other configurations where the 
higher jet velocities now attainable produce the domi 
nant effect on performance. 

Graham, Lagerstrom, Licher, and Beane’ recently 
surveyed the slipstream problem, and extended 
Koning’s results in a form valid (within the original 
assumptions) at all jet velocities. Graham's study also 
included an application of slender body theory to this 
problem. These two theories were compared with the 
experimental data of Stuper,"’ which seemed to lie 
between the two theories. 

In most cases of interest the width of a propeller slip 
stream is of the same order as the wing chord immersed 
in it. With the aspect ratio of the immersed segment 
thus of the order of one, it would seem lifting surface 
theory would be required to properly describe the 
actual conditions. The finding of Graham that the 
available experimental data lay midway between the 
two limiting theories, lifting line and slender body, 
would appear to support this view. 

The difficulty in applying lifting surface theory to 
this problem is the determination of the downwash due 
to the boundary of the jet at points downstream of the 
lifting line. However, the downwash aft of the lifting 
line is of interest in certain wind-tunnel problems, and 
has been determined by Lotz'! and Burgers.'* These 
wind-tunnel investigations were situations where the 
outer velocity was zero and the airfoil lay entirely 
within the jet, and as such were special cases of the 
general problem of an airfoil extending through a circu- 
lar jet. A similar investigation for a closed wind 
tunnel was carried out by Eisenstadt.!* 

This paper develops the general case where there ex- 
ists an outer velocity and where the airfoil extends 
through the jet. This permits the calculation of the 
downwash due to the boundary in the entire field aft 
of the lifting line. Any wing problem in principle can 
then be solved by lifting surface methods. 

As an example, the wing problem investigated experi- 
mentally by Stuper"’ is solved by finite step methods, '! 
based on the approximate lifting surface method of 
Weissinger.!° 

The solution presented in this paper has necessarily 
employed certain idealizations of flows occurring in 
practice in order to make the problem amenable to 
analytical treatment. This fact, discussed in detail 
by the writer elsewhere,’ should be kept in mind in 


applying the results of this paper in design. 


ANALYSIS OF DOWNWASH DUE TO THE JET BOUNDARY 

The problem, as illustrated in Fig. 2, is an airfoil 
spanning a circular jet. This circular jet of velocity 
is bounded by an infinite outer flow of velocity J,. 
The finite wing extends symmetrically completely across 
the jet into the outer flow on both sides. 
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AERODYNAMICS OF 


In the analysis, expressions will be developed for 
downwash due to the jet boundary reaction to infini- 
tesimal vortex elements, which may be integrated to 
give the contributions of finite vortices. The airfoil 
itself will then be replaced by a series of such finite 
horseshoe vortices. The case for such a vortex inside 
or outside of the jet will then be solved separately. 
These horseshoe vortices, illustrated in Fig. 3, will be 
represented in width parametrically, so that combina- 
tions may later be superimposed to represent any arbi- 
trary jet-wing combination. The circulation around 
the horseshoe vortex as illustrated is taken as positive. 

The fluid is assumed to be an ideal, incompressible, 
nonviscous fluid, and this leads to the Laplace equation 
as the governing equation. 

The boundary of the jet is a free surface. The 
presence of different velocities on its two sides gives 
rise to a vortex sheet on the boundary. This vortex 
sheet has associated with it a flow system such that the 
required conditions on pressure and continuity of the 
boundary are satisfied. 

Thus the requirements of the jet boundary deter- 
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mine the flow, aside from the singularities represented 
by the wing and its wake. Two conditions must be 
satisfied on the boundary in a linearized small per- 
turbation theory. First, the pressure of the inside flow 
(region 7) must at the boundary equal the pressure of 
the outside flow (region 0). Otherwise, the acceler- 
ation of the fluid would be infinite across the boundary. 
In terms of the perturbation velocity, v, whose axial 
and radial components are u and u,, respectively, this 


condition is written 
V tt; = Votte (1) 


The second boundary condition is one of continuity, 
and states that the boundary must consist of stream- 
that is, the inside and outside flow must be 


For a circular jet, this 


lines 
tangent at the boundary. 
tangency condition is expressed as 


Vou, = Vit, (2) 
The perturbation flow v is composed of two parts, 


that due to the wing and that due to the boundary. 
The flow due to the wing, or due to the system of horse- 
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Fic. 5. Two-dimensional vortex system of a pair of sym 


metrically spaced horseshoe vortices inside jet 


shoe vortices used to represent it, is known, and may be 
readily determined by the Biot Savart law. The essen- 
tial problem then is to determine the unknown flow due 
to the jet boundary. The sum of the two flows then 
constitutes the total perturbation flow. 

We note the single elementary horseshoe vortex of 
strength I may be decomposed into two halves, one 
even (1) and one odd (2), each of strength I 2, as illus 
trated in Fig. +. These two systems are independent, 
and satisfy the boundary conditions independently. 

This decomposition of the vortex into its even and 
odd parts with respect to x will facilitate the solution 
of the problem, as the even part is two-dimensional, 
and may be readily solved by two-dimensional poten- 
The solution of the odd part, while more 
anti- 


tial theory. 
complicated, may be facilitated by employing 


symmetry conditions. 


Horseshoe Vortex Inside Jet 


The elementary horseshoe vortex inside the jet is 
first decomposed into its even and odd parts. Each of 
these halves satisfies the boundary conditions sepa- 
rately, so we may adopt whatever method of solution is 
most feasible for each case. 

Even Syste m—The even system, as illustrated in Fig. 
This plane flow may then be 
potential 


t, is two-dimensional. 


most siunply by two-dimensional 


In terms of our boundary conditions, the 


solved 
theory. 
condition (1) has vanished with the x component of the 
perturbation velocity. However, the radial boundary 
condition (2) still applies. 

This requirement, for the case of a vortex located in- 
side the jet, may be met by an image system consisting 
of vortices located outside the jet at the inverse points, 
as shown in Fig. 5. 

It has been shown (see reference 5, p. 391) that, when 
a real vortex is located inside a jet, the contribution to 
the motion inside the jet due to the boundary is that 
represented by this fictitious image system outside the 
boundary, modified in intensity by the factor + [(17)° — 


V,2) (Vj)? + V,*)].. Thus when the outer velocity V, is 
zero, the image strength is undiminished. As J, ap- 


proaches |’;, the image strength vanishes with the 


boundary. 
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Also, the contribution of the boundary to the motion 
outside the jet, due to a real vortex located inside the 


jet, in the presence of a boundary, is that of the real 


vortex itself, modified in intensity by the factor 
—((V; — V,)?/(V;? +. V,? When the outer ve 
locity |’, is zero, this boundary contribution is un 
diminished, and completely cancels the real vortex t 
inside. As the outer velocity I’, approaches 1';, the ¢ 
boundary vanishes, and the effect of the real vortex 
inside is undiminished. 

I:ffect inside jet: For a real horseshoe vortex of 
strength [ lying in the plane zs = 0, of small but finite 


span extending between two parametric nondimen 
sionalized points c and d, the even part consists of two 


4 through 


line vortices of strength [ 2 
The image line vortices then extend 


extending to + 
these points. 
through the inverse points of c and d—namely, | c and 

| d, respectively. 

Since the image system replaces the boundary, the 
downwash due to the boundary, for the even system 
may be obtained as the downwash due to the image 
system. This may be obtained by a straightforward 
application of the Biot Savart law. 

For two line vortices through points ¢ and d, as illus 
trated in Fig. 5, the downwash induced at any point 
in the jet n, in the plane of the vortex due to the image 
vortices, may be found as follows. 

For a doubly infinite vortex of strength [ 2 the in 


duced downwash is 


wy) = (1 /2)/2xh | 


where /i is the perpendicular distance from the point 
n to the vortex. Hence, for the image pair located 
at the inverse points of c and d, where the first sub 
script on w(n) refers to the location of the effect point; 
the second to the location of the real vortex, 


r | | | | 
lrr, L(1 d) — n (l/c) —nJ 


and for the image pair located at the inverse points of 


| l l 
imr, | (1/d) + 7 (l/c) +7 


Hence, for a pair of finite horsehoe vortices, symmet 
the 


—cand —d, 


Wjj(n) = 


rically spaced with respect to the axis of the jet, 
downwash due to the image of the even part, induced 
at a point 7 inside the jet, taking into account the 


image strength factor, is 


Effect outside jet: Here the downwash induced at 
point » outside the jet by the boundary due to a pair 
of real line vortices of strength T 2 extending to 4 


through points c and d is 


motion 
side the 
the real 
factor 
iter v< 
is un 
vortex 
V;, the 
vortex 


rtex of 
It finite 
dimen 
of two 
hrough 
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vstem, 
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Ww, (N lr tar.) )[1/(n — « — 


ud that induced by the vortex pair at c and —d is 
Won) = (T/4ar,) , [1/9 + d - |l(n+c)]} 
Hence for a pair of finite horseshoe vortices, sym 
metrically spaced with respect to the axis of the jet, 
the downwash induced by the boundary, due to the 
even part, at a point 7 outside the jet, taking into ac 


count the 1mage strength factor, 1s 


n | - J | | ['/4rr,) X 
i {l/@— @ — [1/(7 — d)] 4+ 
lintd|— l/m@tol | 


Odd System—An elementary horseshoe vortex inside 


the jet will induce a velocity perturbation, the vector 
potential of which in a uniform flow we represent by 
A. In addition, the requirements of the boundary in- 
duce a further flow, the potential of which we represent 
by Q. 

The total perturbation velocity v for the odd system 
2) of Fig. 4is thus the sum 


v VYxA+Vo¢ 5 


The part of the perturbation velocity due to the horse 
shoe vortex itself, v4, represented by the vector potential 
A, is known. The part due to the boundary, Ug, repre 
sented by the potential ¢, is unknown. 


total perturbation velocity v must satisfy the boundary 


However, the 


conditions (1) and (2). Hence we may use these 
boundary conditions to find the unknown potential ¢, 
and from this potential ¢ we may then find the bound- 
ary induced downwash for the odd system. 

We may represent the potential ¢ conveniently by a 
series of the form for p < | 

aT =. 

o - > sin (2p + 1)0 X 


Tv” p=0 


elopti(A)L pti Ap sin &\ dX 0) 


and for p> | 


or 
b> = - >) sin (2p + 1)0 X 
Tw | 0 


, 
Boyt (A) Kopt+i1(Ap) sin EAdX (7 
Jo 

The potential ¢ is antisymmetric with respect to the 
vy plane, since the downwash w = O@ Os is symmetric 
with respect to this plane. Hence for the potential we 
use only a sine @ series. Also we use only the odd 4 
terms as the even @ terms are antisymmetric with re- 
gard to the z plane, whereas the flow is symmetric about 
this plane. 

Also from Fig. +, the downwash w in the odd system 
2) is seen to be odd with respect to x. Hence ¢ is also 
odd in x. We may thus represent ¢ by a sine series in 
v, or in € in nondimensionalized form. 

With this form of the potential ¢ established, we must 


find the unknown coefficients .1(A) and B(A), employing 
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the known velocity due to the vector potential A of the 
horseshoe vortex itself. 
The relation between the vector potential A and its 


resulting \ elocity field V4 iS 


where .1,, .1,, .1, are the rectangular components of the 


vector potential A. For a rectangular horseshoe vortex 
lying in the plane z 0 as illustrated in Fig. 4, the po 
tential of the bound vortex is and the potential of 


the trailing vortex its .! 


For the bound vortex .- |! | 0, and 
Us, = 1, = —(01,/02 
0.1, Or)(Or O 0.1, O08) (0@ O 
or 
It —(sin 6 r,) (0.1, Op cos 6 7,p) (0.1, 06) () 


which we write as the total « perturbation due to the 
vector potential A, since the trailing vortex does not 


contribute to this axial velocity component. Simi 


larly, 
}'Vs 0 
and Wy Rv, 0.1, ON 
so M w sin 4 sin @ ) 0.1, O€ LQ 


Phe radial perturbation velocity component 1 
will also have a contribution from the trailing vortex 
for which .1, = -1 0 


For this case, 


4, LV, () 
VA, jv, O1,/O2 

sin 6(0.1,- O7 +. (cos 6 7 0.1. O86 
Wa, = Bey = —24,/2y 

cos 4(0.1 Or T sin o/r O | Ot 
But My v cos 6 w sin 6 
Hence, Ur, 1 rp 0.1, O86 1] 


Now combining Eqs. (10) and (11) we have the total 
radial perturbation velocity due to the vector potential 


A 
u,, = sin 6 r,) (0.1, O0&) + (1 rp) (0.1, 08 12 


These expressions (9 and (12) give the axial and 
radial velocity perturbations due to the horseshoe vor 
tex whose vector potential is A, in terms of the -!, and 
.1, components of the potential, corresponding to the 
bound and trailing elements of the horseshoe vortex. 

We now proceed to determine these required com 
pt ments of the vector potential. 

Bound vortex elements: We consider an elementary 
horseshoe vortex of strength [ lying in the plane z = 0 
of small but finite span extending between two para 


The odd 


metric nondimensionalized points c and d. 
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Fic. 6. Cosine law geometry. 


part is then of strength ! 2. However, as can be seen 
from Fig. +, the bound vortex is of full strength I’. 

Now the flow due to a horseshoe vortex of strength I 
may be represented by the vector potential 


A= (I 4r) S (dl R’) (13) 


For a differential element of the bound vortex, dl = 
jd6, and R’ = Vi? + (yn — 8)? 4+ ¢2, 4, = A; = 0, and 


dA, = (de tar) {1 VE +- (n a 8)? + £7) 


Now we may represent the reciprocal distance by 


(see reference 16, p. 75) 


l | 
R’ VB&+(y7-B7P +2 


9 Pf 
= | K, [AV (n — B)? + €&7| cos EXdd 
0 


Te 


Furthermore, from Fig. 6, 
R= V(y — B) + f° = VB + p*? — 2Bp cos 0 


Now applying Gegenbauer’s addition theorem for 
K,(AR) (see reference 16, p. 74, and reference 17, p. 
t+) for the case inside the jet, where | 8) < | p 


K (AR) = 1,(\8)K,(Ap) + 


2 ¥ cos poT,(AB)K,(Ap) (14) 
p=l 


For |p| < |8)|, interchange p and 6. Thus |8 p| < I 


1/R’ = (2/z) T,(A\B)K,(Ap) cos EX\dX + 
0 


p> cos pé-(4 7) | 
p=l 


] (AB) K,(Ap) cos ddr 
J 
For a pair of horseshoe vortices symmetrically spaced 
with respect to the jet axis, each of span from ¢ to d, 
the terms for odd p cancel and those for even p combine, 
SO 
ie - ; 
Ay: = = | J | 1,(xa)dat K,(Ap) cos EA\dX + 
redo VW. | 
2 >. cos 2p6 } } I.,(3)dB K»,(Ap) cos eM (15) 
p=1 JV e < 


Trailing vortex elements: Having found A,, we may 
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determine .1, by noting that the divergence of the vec 
tor potential A is zero-—i.e., 


VA=0 (16 


Since we are considering a vortex in the plane z = 0 


1, = 0. Hence, this divergence may be written 


(0OA,/Ox) + (OA,/Ov) = 0 (17 


The expression obtained above Eq. (15) for 4, may 
be differentiated with respect to y, and the result inte 
grated with respect to x, adding a constant of integra 
tion f(p, 8). The result obtained, after simplification 
through the use of recurrence relations for the trigono 
metric and Bessel functions, is 


r > r 
A, = . p cos (2p + 1)0 | J | (/ (AB) + 
Tv” p=0 l 


JV e 


Top 2 (dB) BS Key 


(Ap) sin EAdX (1S) 


In this expression for .1, we have put /(p, 0) = 0, 
since /(p, 9) is even in £, whereas .!, is odd in &, as seen 
from Eq. (13) and Fig. 4. 

As a verification of this expression for .1,, it can be 
shown (see reference 1, p. 21) that the radial velocity 
obtained from it at & = 4 approaches the correct 
value independently obtained from the two-dimensional 
theory indicated by Fig. 7. 

Application of boundary conditions: We now have 
expressions for .1, and -1,, which are known functions. 
We have series representations for the potentials @ of 
the flow due to the boundary, unknown to the extent of 
a coefficient in each of the two expressions. By em- 
ploying the two boundary conditions (1) and (2), we 
may determ:ne these two coefficients and hence deter 
mine the potentials ¢ for the flow due to the boundary 
in the two regions. 

The first boundary condition which applies on the 


surface p = 1 is 


which on employing Eq. (9) becomes 











Fic. 7. Complex potential geometry 
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";[—(sin 0/7r,) (0.1, 0p) — (cos 0 7r,p) (0.1, 00) + (0g 7,0) ]; = 
—V,[—(sin 0 7,) (OA, Op) — (cos 6 7,p) (0.1, 00) + (Og 7,0 


Now, since the form of the vector potential depends only on the location of the real horseshoe vortex and not on the 


effect region, we may write 


»>—-V, + V.{ [—(sin 9 7,) (0.1, Op) — (cos 6’r,p) (0.1, 06)]; — V,[O@ 7,0 —|,|0@ r,0¢ 
On substituting the potentials from Eqs. (15), (6), and (7), and putting p = 1, we have 
; a aor £4 F } of: sin 6 T 
s—V;+ Voi — = T,(\B)dB¢ AK,'(A) cos &AdvA — 
r TT JO { e f Yo T° 
. * cos A r . . 
2 >> cos 2p6 / } T.,(d3)d3 AK2,'(A) cos EAdX + 2 55 (2p) sin 2 pé | I.,(AB ds! x 
p=l Jo t. ly To p= Jo t. y 
a ae: , 
K.,() cos EAdA | — J - S sin (2p + 190 Aoy+(A)Top4i(X)A CoS EdD 
7, FT" 4 0 J 0 a 
-flears., j ; 
= J ~- £2, SR (2p TT 1)6 B pt (A K +] A) A cos tAdX 
ke t 0 J 0 


Now we divide through by —(I° 7*r,), take the Fourier cosine transform with respect to &, and divide through by 4. 


In this expression and the following analysis, the omission of the argument of a Bessel function or the coefficients A 


and B means the argument ts X: 


\V; — Vet ain 6 ‘ 1,(\3)d8 K,’ — 2 > sin 6 cos 20 ‘ I.,(AB ds, Kz,’ + 
e p l e 


P 2p a 
2 S cos @ sin 2p60 T.,(AB dg ~ K ] + 1V,|2 > sin (2p + 196.1 / 
1 ; = ; f \ ai 0 ; 


/ 


Recurrence relations among the trigonometric and Bessel functions simplify this expression to 


2 “&* Tayt (AB 
i hes , = . - pr i\é 
) | = Vt > (2p + 1) sin (2p + | 6k pti | d AB) + 

d p 0 « d AG 

2V, >> sin (2p + 1)0As,+1 J2,41 = 2V, dO sin (2p + 1)0Bz,+1K 

This equation must hold for all values of 6, hence we may equate coefficients of sin (2p + 1)6, where p runs from 0 
to Also, since now all Bessel functions and the functions .! and B are of order 2p + 1, we drop the subscript 
temporarily: 
iV; — Vig (2/4) (2p + 1K [Z(AB) ABld(AB) + 2V;.17 = 21,BA 


The second boundary condition which applies on the surface p = | is 


which on employing Eq. (12) becomes 


—V,[(sin @/7,) (0.1, O&) + (1 7,p) (0.1, 08) + (0g 7,0p)]; = 
—I’,[(sin 6 r,) (O.1, O£) + (1 4r,p) (0.1, 08) + (0g r,Op 


which, as before, may be written 


1—V, + Vit [(sin 6 7,) (0.1, 0&) + (1 rp) (0.1, 00 — |',/0¢ r,0p —|',[0@ r,0p 
and putting p = 1, we obtain an expression which 
Further simplification, 


On substituting the potentials from Eqs. (15), (18), (6), and (7), 
may be simplified as before, in this case taking the Fourier sine transform with respect to &. 


as previously, yields a second equation for the unknown coefficients 


] 20 


1, — Vif (2A) (2p + DR’ [J(AB) AB] d(AB) + 2V,Al’ = 2V,;BR’ 


Results: We now have two equations—(19) and (20)—-which are solved for the two unknown coefficients .1 and B, 


using 
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i 
ed IK’ —I'K = —(1/n 
‘ia 
Substituting the values of .! and B thus obtained in the expressions for the potentials, Eqs. (6) and (7), for the cas¢ 


of a horseshoe vortex inside a jet, where all Bessel functions are of the 2p + 1 order, 











or =. "7 VP? — VotKK'I (Ap) sin Add (® I(AB 
b= y (96 + 1) sin (26 + 1)6 / ( 
— ~ =|} . \=/ PRN PS “— a ad ; ) 
7" 0 J0 V,2l'K — V;7zK \ Jar AO 
aT a ; ij Bat | K (Xp) sin dX (7° [(AB) . 
Pp = ee (2p i 1) sin (2p a | d > dena See ae l ad AG ‘ 
T° 0 J rn | I K = | IK ail A e e) 
Horseshoe Vortex Outside Jet 
The elementary horseshoe vortex outside the jet 
may be treated in a manner identical to the vortex 
inside. 
Even System—As before, the even system, as illus fr 
trated in Fig. 8, is two-dimensional. A similar image / ; 
e e 2° . \ 
system may be used, consisting of image vortices located \ 
inside the jet at the inverse points, as illustrated in * . oe t 7 o——-? ° va 
Fig. 9 -d -c 1 Cc d 
g. 9. i zh 
. ia meg | 
a es 
x Y 
~~ C 
™ ry. . ° 
Fic. 9 Pwo-dimensional vortex system of a pair of syn 


metrically spaced horseshoe vortices outside jet 


For this case, when a real vortex is located outside 

the jet, it has been shown (see reference 5, p. 389) that 

Y the contribution to the motion outside the jet due to 

the boundary is that represented by the fictitious image 

system inside the boundary, modified in intensity by 

the factor —[(V;? — V,?)/(V;? + V,?)]. Thus when 

the outer velocity |”, is zero, the image strength is equal 

in magnitude and opposite in sign to the real vortex 

outside the jet. As J’, approaches J’), the image 
strength vanishes with the boundary. 





c The contribution of the boundary to the motion inside 

the jet due to a real vortex located outside the jet, in 

| (1) Even the presence of a boundary, is that of the real vortex 
itself, modified in intensity by the factor — [(I’) — V 

(I> + 1,”)|. When the outer velocity I, is zero, this 

| boundary contribution is undiminished, and com 

pletely cancels the real vertex outside. As J, ap 

t “Yo proaches |’;, the boundary vanishes, and the effect of 


the real vortex outside is undiminished. 
Effect inside jet: The downwash induced at a point 


lp, n inside the jet by the boundary due to a pair of real 
line vortices of strength I 2 extending to + © through 
Pere y points cand d is 
ie w,(n) = (I /4ar,) ) [1 (d — n)] — [1 /(e — n)]f 
(2) Odd and that induced by the vortex pair at —c and —d is 
wj(n) = (1 4ar,) > [1 (d + n)] — [1 (ec +n) ]j 


Hence for a pair of finite horseshoe vortices symmetri 
—" ; cally spaced with respect to the axis of the jet, the down 
Fic. 8. Even and odd parts of elementary horseshoe vortex ee , 

outside jet. wash induced by the boundary, due to the even part, at 
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, point y inside the jet, taking into account the image points of cand d is 
strength factor, is ; ‘ 
W,,( [ tar,) ;—1 In ld + | In want 
n) = -—[(V; - ] | + | l'/4ar,) X : 
’ and for the image pair located at the inverse points ot 
} l = 1 -_ | c= T 
; ) ' —cand —d, 
lid +n] —il (e+ yn) fj 21 
: : : : rot | ] 
I:ffect Outside Jet: For a pair of real line vortices Wy (N : 1 
Tr + (1 n ld 
extending to +o through points c and d, the image | ” 


line vortices will extend through the inverse points of Hence for a pair of finite horseshoe vortices, sym 


ind d—namely, 1 ¢ and | d, respectively. metrically spaced with respect to the axis of the jet 


The downwash induced at a point 7 outside the jet the downwash due to the image of the even part, it 


due to the image vortex pair located at the inverse duced at a point n outside the jet, taking into account 


the image strength factor, is 





Odd System—As before, we represent the total perturbation velocity v for the odd system (2) of Fig. S by the sum 
v=VxA+Vo ) 


lhe expressions for the known perturbation velocities due to the vector potential A are unchanged, as given in 
Eqs. (9) and (12 

Bound vortex elen:ents: We consider the bound vortex as before, interchanging p and 8 in Gegenbauer’s addition 
theorem in Eq. (14 Hence, for |8 > |p 


K,(AR) = 1,(\p)K,(48) + 2 ¥ cos pol, (Ap) K,(d8 lta 


rhus, for |8/p| > 1, 


| 2 FF : a tr 
,= T,(Ap)K,(A8) cos &\dd 4 > cos pé - T,(Ap)K,(A8) cos EXdX 
R Tr J0 l Tr J) 
For a pair of horseshoe vortices symmetrically spaced with respect to the jet axis, each of span from ¢ to d, the 
terms for odd p ¢ ancel and those for even P combine, so 
‘ty . . ) a . j . , ’ ] 
= K,(AB dp [,(Ap) cos —XdX + 2 >» cos —2pé K..(\8 ds ].,(Ap) cos EXdd | 23 
mLJo (. Jo . : 
Trailing vortex elements: As before, we obtain .1, from .1, (see reference 1, p. 3S), employing the divergence 
theorem (16) and (17). For the present case of the vortex outside of the jet, we use Eq. (25) for .1 The result 
obtained is 
iY z . ni . ’ . | 
A — Dd cos (2p + 1)0 [A2,(AB) + Aap 2(AB dp i Ap) sin EAdA 24 
T” p=0 Jo 


This expression has been similarly verified (see reference 1, p. 41) by showing that the radial velocity obtained 


fromit at = + © approaches the correct two-dimensional value. 

| pplication of boundary conditions: The procedure used above where the horseshoe vortex was inside the jet 
may now be repeated, using the new expressions for .1, and .1, for the present case of the vortex outside the jet (see 
reference 1, p. 44). The same boundary conditions (1) and (2) apply, as do the series representations (6) and (7 
of the unknown potential ¢ of the boundary induced downwash. 

Results: This procedure yields the following expressions for the potentials for the case of a horseshoe vortex out 


side a jet, where all Bessel functions are of the 2p + 1 order 


2T \- ‘ aan eA 4 l ‘|}7(Ap) sin dX (7° K(AB 
Pio = - 2, (2p + 1) sm (2p + 1)0 a = ae | d(X\o Ob 
0 N V,Al’K — V;7IK J r Jar Xd 


Te : co Py? — VY 2tTT’K (Ap) sin dA (7° A(AG 7 
Do >. 2p + 1) sm (2p 1)@ Sea ; a a( XS ib 
Tv 0 ; J | | A = | TK \ e " AO 


Summary of Results 


Even System—The previously derived expressions for the downwash induced by the boundary due to the even 
system are now collected for convenience, writing 1, 1’) = u, whereO <u < 1. The first index refers to the effect 
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region, and the second to the location of the horseshoe vortex. 


Vortex inside jet 


| l— yp? I Pr | 
wi, n) = . ~ — 3a) 
_ 1+ pw*4ar,L(1/d) — 7 Lijec) =— ¥ (ld) + (l(c) +n 
(1 — p)? I ] l l I 
Wj  (n) = — = | _ + _ | (4a 
eren l + |g lor, y= € es d UI ia d U | i c 
Vortex outside jet 
(1 — yp)? I | l l | 
W jo n= - fe = — = (oti 
even l — un irr, d =? = n ad + | ie U) 
1—yp I l l I | 
s i canton _ 4 - + (22: 
even 7 ym n—(l/d) yn-(le) n+ (le) rt 


Odd System—We may now similarly collect the four expressions for the potentials of the flow induced by the 
boundary due to the odd system. These are designated by ¢,,,, where the first index refers to the location of the 


effect point, and the second to the location of the horseshoe vortex. All Bessel functions are of the 2p + 1 order, and 
are of argument A unless otherwise shown. 


Vortex inside jet 


2r re V2 — ViAtKK'T (Xp) sin EddX (2 I(AB) 
¢, = —~ DX (2p + 1) sin (2p + ie i 3 es | d(8) (Ga) 
Tw” p=0 0 HW r J ch vi 
oT : i ViV; K (Xp) sin Ady (° [(AB) ” 
%; =~ dX (2p + 1) sin (2p + 1)0 = <= i d(XB) (Va) 
TT” p=0 0 Al nN J cr / 
Vortex outside jet 
or ; eT Voy. T(Xp) sin —Ek\dX (** K(AB) 
oj = > > (2p + 1) sin (2p + ef = — | a= Mi | d(\B) 6b) 
Tw” p=0 0 All d cr AL 
2r "IVP — V24lI’K(Ap) sin Edd (7 K(AB 
do = = > (2p + 1) sin (2p + 1)0 [ al = eee | d(AB) (7b) 
Tw” p=0 J0 I DN J cr AB 
where We ti(A) = Ve2’K — VPlK’ (25) 


The downwash due to the jet boundary, at points in : 
For 6 = 0, where 0¢ 02 = +(0¢/r08), cos (2p + 1)@ = 


the plane of the airfoil, may be obtained for each of the : ; ‘ ner : 
+1 for all p, and this expression in front of the ¢,,, 


four potentials by taking the derivatives 


integrals is positive. Similarly for 6 = 7, where 0¢ Oz 

Weary = (OGyn/ 02) = = —(0¢/r00), cos (2p + 1) 6 = —1 for all p, and this 
+ (d¢,,, 708) | + when 6 = bs | expression is again positive. 

— when @ = 180 Hence in the plane of the wing, for all », writing r = 


The derivative with respect to 706 of the expression in YoP, and noting that in the plane of the wing p = 7, the 
front of the integrals of ¢,,,, which includes the angular expression in front of the ¢,, integra’s becomes 


function in each case, is 


+(20 /n*rn) D> (2p + 1)? (26) 
p=0 
(21 /rr) >> (2p + 1)? cos (2p + 190 

p= 0 Employing this result, rearranging the integrand for 

In the plane of the wing, which we have chosen as convenience, and, as before, writing 1”, I’; = u, we ob- 

the xy plane, @ assumes only the values of 0 and 7. tain the downwash induced by the boundary due to the 

odd system: 
Vortex inside jet 
a OS . KK’‘I(nX) sin f\dX (°% J(AB) 
Wij ao 1) = > > (26 + 1)? | , = . | d(AB) (Ge) 
si Tron p=0 Jo [1/A,(1/p?) - 1}| —JK rn J cr nV 


K(md) sin EXdX (7° TAB) _ 
d(AB) (We) 
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[he four downwash expressions for the odd system 
may be investigated in the limit as the variable of inte- 
gration \ > by use of the asymptotic series for the 
Bessel functions. 

It can be shown (see reference 1, p. 59) that for large 
\ the integral expression (6c) for w,; will converge at 
least to the order | \* for the case where both the effect 
point » and the vortex point d are at the jet boundary. 
The integral will converge exponentially to a stronger 
degree as either the effect or vortex point is taken farther 
from the boundary. 

The remaining three integrals may be investigated in 
the limit as \ ~ © in a similar manner, and in each 
case they may be shown to converge with the same 
rapidity. 

Asymptotic Correlation of Downwash of Odd and Even 
Systems—Each of the four expressions for the down- 
wash due to the odd system should in the limit as & > 
approach the corresponding expression for the even 
system. that 
each such integral in the limit does indeed demonstrate 


It can be shown (see reference 1, p. 61 
this behavior, thus verifying the expressions obtained 
for the odd system. 

EXTENDING 


TO AN ARBITRARY WING 


THROUGH A JE1 


APPLICATION 


We have determined a basis for finding the downwash 
at any point &, 7 in the plane of the wing due to the 
boundary reaction to a pair of finite horseshoe vortices, 
symmetrically spaced with respect to the axis of the 
jet. The total downwash at such a point due to such a 
pair of horseshoe vortices will then be the sum of this 
boundary induced downwash and the downwash due 
to the horseshoe vortices themselves. 

The downwash due to a field of finite horseshoe vor- 
tices has beey conveniently tabulated in several 
forms,’ '* and so facilitates the calculations of span 
loadings by finite step methods.'* We will therefore 
compute the boundary induced downwash in a similar 
manner so that it may be conveniently combined with 
the referenced tables for the horseshoe vortices to cal- 
culate the span loading for an airfoil extending through 
a jet. 

With the total perturbation downwash thus deter- 
mined over the entire x-y plane, and in particular over 
the field aft of the lifting line, any wing problem in 
principle can be solved by lifting surface methods. As 
an example, we employ the approximate lifting surface 
method of Weissinger’ as a convenient method of cal- 
culating the span loading. 

The lifting surface method of Weissinger in effect 
takes into account the linear variation of downwash 
along the chord. For such a linear variation, the 
wing characteristics may be predicted exactly for in- 
finite aspect ratio by satisfying the boundary condition 


on the wing surface at the 3.4 chord point. Such a 
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Fic. 10. Arrangement of wing-jet vortex pattern (with notation 


used to determine jet boundary downwash coefficients 


linear variation of downwash is the case for both a flat 
plate airfoil and a parabolic are airfoil, and hence would 
be expected to be a good approximation for most con 
ventional airfoils. 

The lifting line approximation assumes no variation of 
Thus the Weissinger lift 
the linear 


downwash along the chord. 
ing surface theory, in taking into account 
variation of downwash along the chord, is in effect 
considering the next term in a Taylor's expansion of 
downwash around the | 4 chord. If further accuracy 
were required, say for a more complex cambered air 
foil, a second term representing the curvature of the 
downwash could be employed. 

In references 14 and 19 the wing span is subdivided 
into a finite number of horseshoe vortices, the spans of 
which are nondimensionalized with respect to the wing 
span. For the jet problem, it is more convenient to non- 
dimensionalize the lengths with respect to the jet 
radius. The wing data and the jet data may thus be 
conveniently combined when the wing span is an in 
teger multiple of the jet radius; for noninteger multiples 
the method may be adjusted at the wing tips to provide 
suitable data there. 

To provide a general and flexible method useful for an 
arbitrary jet-span combination, we arrange the wing 
as illustrated in Fig. 10 (see also reference 14, Fig. 1, 
and reference 19, Fig. 6). The wing is replaced by a 
system of .V horseshoe vortices along the quarter-chord 
line, the spans of which are so chosen that their semispan 
sis one-fifth of the jet radius r,. Thus five such horse- 
shoe vortices are located inside the jet. 

Sufficient additional horseshoe vortices of the same 
semispan s are placed outside the jet to completely 
represent the wing span. Thus any arbitrary ratio of 
jet diameter to wing span arrangement may be repre 
sented, simply by employ:ng a sufficient (odd) number 
of horseshoe vortices to represent the wing span. 
Reference 14 includes smaller span vortices at the tips 
which may be used if a further refinement in accuracy 
is required. 

An equal number of downwash points is taken along 
the three-quarter chord line, at the midpoint of each 


elementary horseshoe vortex. The downwash velocity 
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used to determine lift distribution). 


from the total vortex system and the boundary is then 
equated to the component of local velocity normal to 
the wing chord at each such downwash point. Appli- 
cation of this tangent flow boundary condition for a 
symmetrical loading provides a set of (V + 1) 2 simul- 
taneous equations in the (V + 1) 2 unknown circula- 
tion strengths across the span, since the circulation 
strengths of symmetrical vortex pairs will be identical. 
Solution of this set of equations provides the span 
loading. 

li the ratio of wing span to the jet diameter is so 
large as to produce an unwieldy number of simultaneous 
equations, the number of horseshoe vortex pairs and 
hence the number of equations may be reduced by em- 
ploying larger horseshoe vortex spans outboard where 
their effect on the jet boundary is less pronounced. 


Downwash Due to the Wing 
The downwash velocity in the plane of a rectangular 
horseshoe vortex due to the vortex itself is given by the 

expression 
w(x, 9) = (1 /4m) (F(X, ¥) 5] (27) 


where F is defined in reference 14, and the f and j are 
expressed in horseshoe semispans. The values of 
F(x, ¥) are conveniently tabulated in references 18 and 
19. 

Distributing an odd number .V of horseshoe vortices 
having .V downwash points across the wing span (Fig. 
10), the downwash velocity at any of the downwash 
points P, resulting from the V horseshoe vortices is 


\ 
w(t,, ¥,) = (1/47) x (Te/ Slr, 9,) (28) 
n= 1 
where Xx, = p, — Pal 99 
P (Zo) 
Y= d- Gn$ 


and the ¥, 7, p, g notation, and the vy, m subscript nota- 
tion, as used in this section, is illustrated in Fig. 11. 
The v and » are the effect and vortex points, respec- 
tively, counted from the left wing tip. 


Since b 2 = Ns, we may write (30) 
. 

w(X,, ¥,) = (N/4r) DS [P,,/(b/2) F(R, ¥,) (31) 
n=1 


For cases of symmetrical loading and geometry, the 
downwash at any of the downwash points P,, becomes 
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with the coefficient F,,, being given in terms of geomet 


rical distances measured in horseshoe semispans, 


~~, 
F{(p, — Dy), (> + Un) | (33 
and the prime on the summation sign indicating that 
for the (V + 1) 2 term, representing the single horse 


shoe vortex at the center of the wing, 


Fyun+i/2} = Fl(p, — pr), (Gq, — 9)] 34 


This condition may be applied at each downwash 
point to provide a set of [LV + 1) 2] simultaneous 
equations in the [(V + 1)/2] unknown circulation 


oy) 
a 


strengths from Eq. 


wel 


Downwash Due to the Jet Boundary 


The expressions for the downwash due to the bound 
ary have been obtained as follows, where the first sub 
script refers to the location of the effect region, the 
second to the location of the real vortex, and the even 
or odd notation refers to the component of the horse 
shoe vortex: 


Even System Odd System 


Expression Eq. No Expression Eq. No 
W jj 3a 6c 
Ty (4a tt (7c) 
(2la Ww; 6d 
22a (7d 


“ ~ 


These expressions are all nondimensionalized with 
respect to the jet radius 7,, We may combine this 
boundary downwash with the horseshoe downwash at 
any particular point, taking note of the different co 
ordinates of the same point in the two systems, due to 
their being nondimensionalized on different lengths, and 
having different origins. 

Furthermore, by employing the relation 


¥, = Bs (36) 


we may write the boundary downwash expressions in a 
form where they may be conveniently combined with 
wing downwash expressions—1.e. 

w(t, n) = (T'/4r) [G(E, n) 5 (37 
where we write the G(£, 7) as a function of £, 7 to em- 
phasize their being nondimensionalized with respect to 
the jet radius 7,. By this relationship we define the 
jet boundary downwash coefficient G,,,(£, 7) as 


G,,(&, 7) = (4rs/T) Wy» \ 


Str 


n (38) 


The quantities G,,, may thus be written in terms of 
the boundary induced downwash expressions pre- 
viously derived, with the aid of Eq. (36), as shown be- 


low. 
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Che jet boundary downwash coefficients (even system) are 


Vortex inside jet 
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Vortex outside je t 
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fa Oe =; ; ~ + ‘ | 21b 
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Tr S1+ pty — (1/e n— (ld n+ (1/d "+ (/o)J 
Note that 
G in =G ou 
eren ere? 
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The jet boundary downwash coefficients (odd system) are 
Vortex inside jet 
S q KK'I(mX) sin &\ dd mer TONG 
fs (t,n) = - +k (2p + 1)? , ; : ) =f d(A\g be 
oun ITN p=0 Jo [1/,(/u?) — 15] — ALK! Ja Ag 
S a oan l K(nd) sin tkdX (7° [AB 
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Vortex outside jet 
S al : l I(md) sin EAdX (7% K(AB ; 
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ona OTN p=0 J0 es r; (1 a) = wi IK nN J or Ag 
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All Bessel functions are of 2p + 1 order and, unless otherwise shown, are of argument X. 

These jet boundary downwash coefficients may be employed in a manner similar to the wing downwash coefficients. 
Thus, for an odd number NV of horseshoe vortices having V downwash points across the wing span (Fig. 10), the 
downwash velocity at any of the downwash points P, resulting from the boundary reaction to the V horseshoe vor- 


tices 1s 
; By Gyn (41 


for symmetrical loading, the prime on the summation again indicating the unique character of the center horseshoe 
vortex. Here the G,,, is the sum of the even and odd parts. As before, this expression may be writteri 


[ NV Cae fT. 
WE 1, = } By G, t 
tor n=1 b/2 1,2 V+1)/2 
Total Downwash Due to a Wing Jet Combination 
For small angles, the tangent-flow boundary condition 
w(x,, ¥,) = V, sin a, = Va, (43 


may be applied at each of the downwash points. Thus, combining Eqs. (35) and (42), and using Eq. (43), we have, 
as our final set of (NV + 1) 2 simultaneous equations in the (V + 1) 2 unknown circulation strengths, 


N \ u 1) 2 if 
| Vea = WwW, = by : | Fn ni Gul | (44 
de gat 0/2 Te e V+1)/2 


This set of equations is written in terms of the wing span b, the wing downwash coefficient F,,, (see reference 14 


and the jet boundary downwash coefficient G,,,. 
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The values of F are obtained from reference 18. The 
values of G are computed from the formulas for jet 
boundary downwash coefficients derived in this paper. 
Both the even and odd G coefficients must be used in 
the above equations. 

To facilitate the employment of the method derived 
in this paper for the determination of characteristics of 
an arbitrary wing extending through a jet, the G co- 
eflicients, both even and odd, were programed in a 
parametric manner for computation by electronic data 
processing equipment. The required coefficients may 
thus be readily obtained for any arbitrary wing-jet 
geometry and for any desired jet velocity ratio. 


COMPARISON OF THEORY WITH EXPERIMENT 


General 

The present theory provides a means to determine 
the lift distribution over a wing extending through a 
circular jet. In order to provide a comparison of the 
theory with experimental data, measured values of the 
lift at successive spanwise points must be employed. 
There is a dearth of such experimental data, particularly 
at higher values of the jet velocity ratio, where the 
effect wou!d be expected to be large, and which would 
accordingly be of interest. 

One such set of experiments is that reported by J. 
Stuper,'”’ employing velocity ratios that are fairly low, 
yet probably adequate to check the theory. Fig. 16 
of Stuper’s paper shows the experimental results for a 
jet velocity ratio uw = 0.735, and for several angles of 
attack, and also shows the comparison with Koning’s 
theory.® 

Stuper employed a rectangular wing 20 by SO cm., 
thus having a geometric aspect ratio of +. His use of 
circular end plates 32 cm. in diameter gave his wing 
an effective aspect ratio of 5.25 (see reference 12, p. 
211). The jet was 12 cm. in diameter, so the ratio of 
jet diameter to wing chord was 0.6, which is thus the 
aspect ratio of the wing segment immersed in the jet. 

Graham (see reference 9, Fig. 5.5) reproduced the 
data of Stuper’s Fig. 16, and included a comparison 
with slende: body theory. As observed previously, 
within the jet the experimental data seemed to lie be- 
tween the two limiting theories, lifting “ine and slender 
body. Outside of the jet, neither theory represented 
the experiments very well, both lying below the experi- 
mental points. Also neither of these theories gave the 
sharp discontinuity in loading found experimenta'ly at 


the jet boundary. 


Method of Computation 


In order to provide a comparison of the present 
theory with Stuper’s experiment, the Weissinger appli- 
cation of the theory developed in the previous part 
of this paper was employed. The wing-jet geometry 
reported in Stuper’s paper was represented by the ar- 
rangement of horseshoe vortices illustrated here in 
Fig. 10. 
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The segment of the wing in the jet was represented by 
5) horseshoe vortices. Stuper’s effective aspect ratio 


of 5.25 gave an effective span of 105cm. Hence, 


Nip = (6/2) corcaten/%, = Be.0/6 = 8.75 (45 


which determined the spanwise extent of the vortex 
pattern as shown in Fig. 10, thus representing the wing 
by a total of V = 45 horseshoe vortices. 

To use the Weissinger lifting surface method, the 
downwash due to each of the horseshoe vortices, both 
from the vortex itself and from the boundary reaction 
to the vortex, must be determined at the downwash 
points located at the intersection of the centerline of 
each such elementary horseshoe vortex and the 3/4 
chord line of the wing. 

The F downwash coefficients, representing the down- 
wash due to the wing itself, are obtained for the vortex 
pattern of Fig. 11 from the tab es of references 18 and 
19. In using these tables, the value of the 3 4 chord 
line must be expressed in terms of the semiwidth s of 
the horseshoe vortex used. On this basis, 


Sx 2/3 = ,—10 [(b/ 2) 5] = 
= 


(—10 X 45)/52.5 = —8.75 (46) 


is the constant value of the 3/4 chord line used in the 
Stuper experiment. 

The required values of the F downwash coefficients 
at this constant value of ¢ are most conveniently ob- 
tained from the tables in reference 18. In those tables, 
the distance in the spanwise direction from the line of 
symmetry of a given horseshoe vortex, expressed in 


terms of the semiwidth s of the horseshoe vortex, is 
y = |2(v — n)| (47) 


The values of the F wing downwash coefficients, for 
the case f = —S.75 used in the Stuper experiment, 
have been taken from reference 18 and are presented in 
Table 1. 

These F downwash coefficients are of the order of 
unity, and hence were computed in the tables of refer- 
ence 18 to 3 decimal places, in order to provide accu- 
racies of 3 significant figures. Hence the G downwash 
coefficients, representing the downwash due to the jet 


TABLE | 

Wing Downwash Coefficients 
(% = —8.57, see reference 18 

0 +4 (14 34 —() OO 
2 —1.320 36 —0.002 
$4 —().255 38 —0.002 
6 —(). 104 $() —() OO2 
§ —0.055 $9 —() O01 
10 =—0.033 14 —0.001 
12 —(), O22 16 —0 OO1 
14 —(0.016 iS —O0O.001 
16 —0.012 50 —0.001 
IS —0.009 52 —0. 001 
20 —(). 007 54. —0.001 
yi 4 —() O06 56 —(). OO] 
24 —0.005 58 —0.001 
26 —(0. 004 60 —(). 001 
28 —0.003 62 —(0.001 
30. —0.008 64 —().001 


32. —0.002 66 —0.000 
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ted by 
ratio 
TABLE 2 
Ls Jet Boundary Downwash Coefficients 
2 & = —1?/3; pw = 0.735 
‘ortex n B Geve G n 3 eC. . G 
wing 0 0 0). 024" 0.021* > 0 0.8 0.008 0 008 
: 0.4 0.048 0.042 1.2 0.015 0.011 0.011 
O-8 0 O48 0.0438 1.6 0 OO6 0 O04 
|. the 2 0.005 0.006 2 0 0.004 0.002 
1.6 0.008 0.0038 ,e 0.002 0.001 
both 2 () 0 002 0 002 28 0.002 0.001 
ction 94 0.001 0.001 3 9 0 O01 
as > me 0.001 0.001 3.6 QO OO] 
wash 3 2 0.001 +0 0.001 
ne of 3.6 0.001 44 0 O01 
3/4 0.4 0 0.0247 0.021* +8 0.001 
0.4 0.052 0.047 2 4 0) 0.0017 0.001" 
OX 0.066 0.061 (). 4 0 OO] 0.001 
own- 1 2 0.008 0.009 O.8 0.002 Q OO? 
1.6 0.004 0. 004 : 2 0.009 0.006 
ortex 2 0 0.002 0.002 16 0 004 0.002 
. and 2.4 0 OO] 0.001 2 0 0 OO? 0 OO] 
2.8 0.001 0.001 9 4 0. O02 0 OO] 
hord 3.2 0.001 2 8 0.001 
s of 3.6 0.001 . 0.001 
0.8 0 ().025* 0). 022* 3.0 0.001 
0.4 0.069 0.064 (0.063 +.0 0.001 
O.8 0.2389 0.238 (0.179 2 § 4 0.001 0 OO] 
. 0.031 0.0384 (0.024 O88 0. O01 OOO] 
(46 1.6 0.007 0.007 12 0.006 0.003 
3s © 0 0038 QO OO8 1.6 0 O08 0.001 
| the 2.4 0.002 0.002 2 0 0. O02 0.001 
2.8 0.001 0. O01 2 4 0 OO] 
zs. 2 0.001 O OO] 2 8 0 OO] 
ents 3.6 0.001 oe 0.001 
$0) 0.001 
ob- 3.2 0.4 0 OO] 

: 2 0 0.003" 0.003" O 8 O OO] QO OO] 
ples, 0.4 0.008 0.009 1? 0 O04 0002 
e of 0.8 0.031 0.0383 (0.023 1.6 0 OOP? QO 001 

. -. 0.180 —(). 182 | 0.110 2 0 0. O01 
1 in 1.6 0.033 0.028 (—0.026 9 4 0.001 
2.0 QO O15 —() O11 > = 0.001 
9 4 -() QOS -0 OOS 
; ‘ 3.6 0.4 QO OO] 
(47) 2 0 O06 0 003 
, 5 0. OO4 0 OO2 U.6 0.001 
> 2 2 : 3 0. O03 QO O00] 
° ,.6 0 OO8 0 OO] a . 4 
lor 10 0.002 0.001 1.0 U.UUS ied 
2 0 0 OO] 
ent, $ 4 QO. 002 QO. OO] x 
ij 18 0.002 0.001 -.4 0.001 
oe Se 0.001 $+.) 6.3 0 O01 
5.6 0.001 1.2 0. O02 0 OO] 

. 6.0 0.001 1.6 0.001 
- Ol 6.4 0.001 20 0.001 
fer- 6.8 0.001 » 4 0.001 

ee 0.001 is : ’ 
Ccu- 7.6 —() OO] +4 I - J. Ue mre 
ash 8.0 0.001 1.9 a 
a > 1) (nH) 
jet 1.6 () O OOL* 0. OO2* 1s 19° 0 00° O01 
Q).4 QO OO4 0.004 “ = 
= = 1.6 0 OO] 
Os QO O07 0.007 
Q- P > () 0 OO] 
ee 0.035 0.030 0.028 - 
1.6 0.012 0.009 5.2 e- 0.001 
2 0 0 OO6 QO OO4 1.6 0.001 
2.4 0 OO4 -() OOD mage” . 
) ] me 
28 -() 003 -(). 001 o.f -_- 0.001 
3.2 -(). 002 0.001 1.6 0. 
3.6 0.002 -() OO] 6.0 . 0.001 
$0 0.00] 1.6 QO OO] 
1 4 0.001 . ‘ 
tz 0.001 6.4 1.2 0.001 
5.2 —0 001 68 1.2 0.001 
5.6 0.001 — :* 0.001 
6.0 0.001 sldhcos oe 
2.0 0 0.001" 0. O00O1* 4.0 Bae 0.001 
0.4 0.002 0 O00”? S00 . = 0 OO] 
* All values are based on a pair of symmetrically spaced horseshoe vortices of span 43 = 0.4, except for the case 8 = 0, where Ad = 
a . . . ty ¥ . pé . . - I . ‘ 
1.2 rhis situation results from the representation of the wing by pairs of symmetrically spaced horseshoe vortices everywhere except 


it the center of the wing, which is represented bv a single unique vortex 
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boundary, must similarly be computed to 3 decimal 
places, as the F and G coefficients are to be summed at 
each downwash point. 

The G downwash coefficients were computed from the 
expressions derived in this paper. In using these ex- 
pressions, the value of the 3 4 chord line must be ex- 


pressed in terms of the jet radius 7, used. On this basis, 


& = x/r, = (—10/6) = —17/; (48) 


is the constant value of the 3 4 chord line used in the 
Stuper experiment. 

The definite integrals in the G, ~ functions were evalu- 
ated by Filon’s Simpson's rule, employing electronic 
data processing equipment. The computing was pro- 
gramed so that yu, &, 8, and n were parameters, and in 
addition the mesh of the variable of integration \ and 
the upper limit of the integral could be specified. Four 
values of the summation index p were employed, which 
included Bessel functions through order seven. 

Computations were carried out only for the constant 
values of & = —1* 3; and w = 0.735, corresponding to 
Stuper’s experiment.’ For this case trial runs showed 
that an interval of AX = 0.5, together with an upper 
limit of \ of about 12, provided sufficiently rapid con- 
vergence to give the required accuracy of 3 decimal 
places, although an upper lim t of X = 25 was used for 
the most slowly convergent cases where both the vortex 
and the effect point were adjacent to the jet boundary. 

For the cases where both the effect point and the 
vortex point were adjacent or close to the jet boundary, 
the convergence with respect to p was inadequate to 
give the required accuracy of 3 decimal places for the 
four values of the summation index pemployed. There 
were eight such cases. 

These eight cases were then extended by employing 
the Aitken V* process (see reference 20, p. 126) to im- 
prove the rate of convergence of the sequence of partial 
sums of the series. This process was chosen because it 
is exact for a geometric series, which the present sums 
of integrals approach asymptotically for large values 
of the summation index p. 


23 
DS Fin + Gin + Gin} 
=] win 
In this equation the F,,, coefficients are as defined in 
Eqs. (33) and (34). Table 1 may be used together 
Table 2 


may be used together with Figs. 10 and 11 to provide 


with Eq. (47) to provide the F coefficients. 


the G coefficients. 

The 23 simultaneous equations w.th the 23 unkuown 
circulation strengths T’,, so constructed were solved by 
an elimination method. The circulation strengths 
obtained were converted into local lift coefficients, 
based on ,, by the expression 


Cr = licat’ (1/2) pVi2¢ = p( VT )icar/(1/2)pV,2c (50) 


These local lift coefficients, for the wing and jet com- 
bination, were computed as outlined above for the case 
The local lift coefficients 


of an angle of attack of 12°. 


. 
rs = 


The computed values of the G jet boundary down 
wash coefficients for the case of § = —1* 3 and uw = 
0.735, used in the Stuper experiment, are presented in 
Table 2. 
were extended by the Aitken V* process, the extended 
values are those presented in the column, and the par 
tial sums of the first four terms from which the extended 


For the eight cases where the G,,¢ coefficients 


values were obtained are shown in adjacent parentheses. 

The G coefficients of Table 2 were at this point avail- 
able for addition to the F coefficients of Table 1 in 
accordance with Eq. (44), to provide a set of linear 
simultaneous equations from which the unknown circu 
lation strengths I’,, could be determined. 

It will be recalled that the formulation of the problem 
leading to Eq. (44) was based on symmetrical loading 
across the wing span. Hence the circulation of sym- 
metrical vortex pairs was identical. This permitted 
a reduction of the number of unknown circulation 
strengths from N = 45 to (VN + 1)/2 = 23, the odd 
number arising from the unique center vortex. 

Hence a solution for the symmetrical Stuper wing 
was obtained by a set of 25 simultaneous equations in 
the 23 unknown circulation strengths of a half of the 
wing on one side of the centerline. In constructing 
these equations, the F coefficients for symmetrical 
vortex pairs were added together. 

The G coefficients were derived for symmetrical vor 
tex pairs, and so are already in the proper form for use 
in the 25 simultaneous equations. However, the G 
coefficients were computed separately for their even 
and odd parts, and so the sum of these two parts was 
used at each downwash point. 

The common factors in Eq. (4+) may be transposed 
to the other side of the equation, giving, for the Stuper 
wing outside the jet, at 12° angle of attack, 


[4(b/2)/N]Va = [(4e X 1.722)/45] X 
98.425 X 0.209 = 9.915 ft.? ‘sec. (49 


Inside the jet, where the velocity is increased by a 
factor of 1.36, this term becomes 13.485 ft.? see. 
Hence Eq. (44), for the case of the Stuper wing at 12 
angle of attack, may be written explicitly as 


9.915 (vy = 1-20) 


13.485 (v = 21-23) |, — 1-9 


were then linearly reduced for the angles of S° and 4 

From these values, for each angle of attack, were sub- 
tracted the local lift coefficients of the wing alone. 
The differences, which are the increase in local lift 
coefficients due to the jet, were added to the experi- 
mental lift coefficients of the wing alone obtained by 
Stuper (see reference 10, Fig. 16), and are shown in 


Fig. 12. 


Results 


The theory of this paper has been applied to the 
wing-jet combination investigated experimentally by 


Stuper in reference 10. The correlation, presented in 
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Fig. 12, has been found to be excellent for that portion 
of the data presented by Stuper which is felt to be valid. 

Che present theory agrees well with the experimental 
data of Stuper in the jet, which is the region of greatest 
effect and hence of principal interest. The strong de- 
crease in the experimental lift observed immediately 
outside of the jet is explained by Stuper as being due to 
the friction boundary layer of the fan enclosure. 

rhe large increase in lift obtained experimentally in 
the region farther from the jet may be explained by the 
The jet 


has a free boundary as it passes over the wing, and a 


test setup shown by Stuper in his Fig. 20. 


free boundary for a distance of one quarter of the chord 
upstream from the leading edge. But for about the 
next two chord lengths upstream the jet has a solid 
boundary, in the form of the fan enclosure. The pres- 
ent theory, as well as Koning’s theory,° is based on the 
boundary of a free jet, and would not be expected to 
agree with the results of a solid boundary. 

The effect of this solid boundary upstream will re- 
verse the sense of circulation of the image vortices for 
the upstream region, and similarly reverse the direc- 
tion of the downwash induced by such image vortices 
in their representation of the boundary. From the 
geometry of the wing-jet arrangement, it is seen that 
the angle subtended by the fan enclosure, and hence 
the relative effect of the solid boundary, will have 


a maximum at about six jet radii outboard from the 
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spanwise lift distribution of wing extending through circular 
jet. 

jet centerline along the 3 4 chord line. The expert- 
mental points of Stuper show such a maximum, as seen 
from Fig. 12. 

Furthermore, the lift increment obtained by the 
theory is the difference of opposing effects. _Examuna- 
tion of these individual effects for the present case of 
the Stuper experiment has shown the lift increment 
to be a relatively small difference of large effects. Hence 
the large effects at appreciable distances from the jet 
boundary obtained experimentally by Stuper are felt 
to be understandable in view of the experimental test 
setup he employed. 

The strong decrease in experimental lift observed 
immediately outside of the jet could also be partly 
caused by the solid boundary upstream, although the 
dynamic pressure drop in the transition region between 
the jet and the undisturbed flow caused by the friction 
boundary layer of the fan enclosure would seem to offer 
an adequate explanation. 

The Koning theory differs from the present theory 
essentially in that Koning considers only what in this 
paper has been termed the even part, and disregards 
the odd part. Hence the difference in the results ob- 
tained by the present theory must be attributed to the 
inclusion of the odd part. The several theories are 
compared in Fig. 13, which includes Koning’s theory, 
as computed both by Stuper (see reference 10, Fig. 16) 


and by Graham (see reference 9, Fig. 5.5), as well as 
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slender body theory, as computed by Graham. It is 
seen that the present theory is in substantial agreement 
with the valid experimental points in the jet, and lies 
between lifting line theory and slender body theory, as 
expected. 

As may be observed in Table 2, for the present case of 
Stuper’s experiment the jet boundary downwash co- 
efficients for the odd part were of the same order as 
those for the even part, particularly for the spanwise 
points of major contribution near the jet boundary. 
Hence the agreement of the present theory with the 
experimental results of Stuper must be attributed to 
the inclusion of the odd part. 


CONCLUSIONS 


A strong correlation of the theory with the available 
experimental data has been demonstrated. The rapid 
approach of the downwash induced by the jet boundary, 
caused by the odd parts of horseshoe vortices repre- 
senting a wing, to its asymptotic even value down- 
stream of the lifting line emphasizes the importance of 
the odd part in obtaining this correlation. 

The wing segment immersed in the jet in the Stuper 
experiment, used for comparison with the theory, was 
of aspect ratio = 0.6. For this aspect ratio, and for the 
spanwise points of major contribution near the jet, 
the boundary induced downwash at the 3 +4 chord 
points used in the Weissinger analysis was almost the 
Trefftz plane value. The conclusion is clear that for 
such low aspect ratio immersed segments the use of 
lifting line theory introduces serious error, and that an 
approximate lifting surface theory such as that of 
Weissinger is required. 

For wing segments immersed in the jet of aspect 
ratio of the order o one, the downwash over the wing, 
and in particular at the 3 4 chord point, will be in gen- 
eral intermediate between the value at the Trefftz plane 
and half that value at the lifting line. This variation 
is large, and the lift distribution predicted by the theory 
is quite sensitive to the value of the downwash used. 
The present theory provides a means to determine the 
proper downwash, and more accurate predictions of 
lift distribution may be obtained thereby. 
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The Hydrodynamics and Heat Conduction of 
a Melting Surface 


GEORGE W. SUTTON* 
General Electric Company 


SUMMARY 


In previous investigations of surface melting the assumption is 
usually made that, upon melting, the molten liquid is immedi 
itely removed. In this paper, this assumption is relaxed and 


the governing differential equations are presented. Solutions 
are Obtained for the aerodynamic melting of Pyrex glass at a 
stagnation point, giving the velocity and temperature profiles 


of the molten film, and its thickness. The results indicate that 
the presence cf the film cannot be ignored because the melting 
rate is more rapid than that predicted by simple theory In 
spite of the additional complexity of the model, the melting 
rate is found to be linearly dependent upon the heat-transfer 
rate, and this relation is relatively insensitive to Mach Number, 
altitude, mass-transfer rate, and interface shear stress for the 


range of values considered 


SYMBOLS 
C,, = specific heat at constant pressure 
F = similarity stream function 
k = coefficient of thermal conductivity 
K = density ratio across normal shock wave 
m = masstransfer, —pi 
P, = Prandtl Number, C,y/k 
p = pressure 
g = heat-transfer rate 
R, = stagnation point radius of curvature of body 
7 = temperature 
u,v = fluid velocities 
v, ¥ = coordinates 
“ = viscosity, ml 7 
T = shear stress 
© = dimensionless temperature ratio, (7 — 7))/T* 
€ = (for planar flow 
€ = 1 for axially symmetric flow 
p = density 


= dimensionless distance 


Subscripts and Superscripts 
interior of surface material 
m = melting point 
= outer edge of gaseous boundary layer 
= ahead of normal shock wave 
= gaseous boundary layer - liquid film interface 
= stagnation point 
liquid-solid interface 


= reference 


INTRODUCTION 


A LARGE NUMBER of investigations have been con- 
ducted recently to determine the aerodynamic 
Received June 6, 1957 
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heat transfer to vehicles traveling at high speeds. The 
these studies indicate that at altitudes at 
that is, below 250,000 ft. 


especially in the region of the 


results of 
which air is continuum 
the aerodynamic heating, 
stagnation point, is sufficiently great that some method 
of heat absorbtion or cooling is required. However 
if the flight time is sufficiently short, as is the case for 
rapid exit or entry into the earth’s atmosphere, it may 
be possible to allow a portion of the surface to melt 
and vaporize. 

In order to design such a system accurately, it is 
necessary to predict precisely the rate at which the 
melting surface recedes. Any such theory should 
consider the effect of the molten liquid film which 
previous analyses of the problem of surface melting 
have ignored. ! 

The object of the present paper is to determine 
quantitatively the effect of such a molten liquid film 
in laminar flow at a stagnation point. The general 
equations are presented, and the results of some typical 
calculations are given in which the viscosity of the 
melt is allowed to be a function of temperature. The 
assumption is made that the system is in quasi-steady 
state; this will not affect the generality of the con 
clusions. 

The coordinate system is shown in Fig. 1; it is at- 
tached to the interface between the molten material 
and the gaseous boundary layer. It is assumed that 
the rate of change of the vehicle shape due to melting is 


small. 


BOUNDARY CONDITIONS 


The assumption is made that the molten surface 
material is acted upon by aerodynamic heating and 
skin friction 7; In addition, there may be mass 
transfer into the gaseous boundary layer, caused by 
evaporation or by pyrolysis of the molten liquid film. 
In the former case, mass transfer is controlled by the 
vapor pressure of the material and by diffusion in the 
boundary layer; in the latter case, the ratio is fixed by 


the chemistry of the material. It is convenient to 


choose the following boundary conditions at y 0 
v(0) = —mi/p, (Ou/Oyv) (0) = —7;/ui, T(O) T; 
(1) 


Two additional boundary conditions are necessary. 
They may be chosen in the interior of the surface ma- 


terial as follows. For a material with a definite melt- 
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ing point the tangential velocity is zero at the melting 
point, hence, 
at y Vs: u o. ff pe (2) 
However some materials, such as glasses, do not 
exhibit a definite melting point. The boundary condt- 
tions for this class of materials are as follows: 


lim u(y) 0, lim 7(y) T» (3) 


It is implicit in these boundary conditions that the 


material is semi-infinite. 


(GOVERNING EQUATIONS 


If the Reynolds Number based on a typical dimen- 
sion is large, then the governing equations are identi- 
cal to the usual laminar boundary-layer equations. 


These equations are: 


Continuity 


(0/Ox) (ure®) + (0/Oy) (z7re°) = O (4) 
Momentum 
pu(Ou,/Ox) + pv(Ou/ Oy) —(Op/ox) -- 
(0/Ov) [u(Ou/Oy)]) (5) 
0 Op Oy 


Energy 
pC,|u(OT'/Ox) + v(OT/Oy)|] = u(Op/ox) + 
(0/dy) [R(O7T'/Ov)] + u(Ou/dOy)* (6) 


If the Reynolds Number is of order of magnitude 1, then 
the convective terms may be omitted from Eq. (5); how- 


ever, the convective terms remain in Eq. (6) if the 
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Prandtl Number of the molten liquid film is large. For 
most materials, and especially the example considered 
herein, the variations in density, specific heat, and 
thermal conductivity are small in comparison to the 
variation of viscosity with temperature hence, they 
are assumed constant. Thus, the Prandtl Number will 
vary with temperature. The following incompressible 
transforms which include the Manger transform® are 
next used to convert Eqs. (5) and (6) into ordinary 
differential equations 


p 92)1/2 € 
6 [p/(25) UAX)o (X)¥ (S) 


The following is used as the stream function: 
y (25)'/"F(¢) (9) 


where pure Oy Oy, pure’ —(Ow/ Ox) (10) 


Also, let 0 (1. — 7,)/T* (a3) 


Eqs. (5) and (6) then become 


\omentum 
[(u/u*) Fee le + Fee + Bl(p./ p) — F; Q (12) 
Energy 
FO, + (R/C,u*)Ore + (ue?/CpT*) X 
[(u/ u*) Fee? — B(p./ p)F, QO (13) 


where 6 = 25(0u, OS), u,is 12 at an axially symmetric 
stagnation point and is zero for zero pressure gradient. 
For stagnation point calculations, the last term in Eq. 
(13) may be neglected because u,7 << C,7*. In the 
remainder of this note, only stagnation point flow on 
a blunt body will be considered for two reasons: (1) 
it is of foremost technical interest, and (2) the flow 
is most likely to be laminar at this point. The relations 
between the physical variables and the dimensionless 
variables at a stagnation point are given by the follow 


ing relations: 


ad — |(e + 1)u*(Ou, Ox), “Y) (14) 

F, uu, (15) 

F.. —(u*/p) [Ce + 1)u*p(Ou,/Ox),° : “(Or Ox) (16) 
0 (T — T>)/T* (17) 

0; —(1/k7*) [(e + 1)p(Ou,/Ox) ux g (18) 

c [(e + 1)p(Ou,/ Ox) u* 1)! “y (19) 


The appropriate boundary conditions for Eqs. (12) 
and (15) are obtained by the substitution of Eqs. (1), 
(2), and (3) into Eqs. (14)-(18). The melting rate is 
given by lim o(y). 


APPLICATION TO A PARTICULAR CASE--THE MELTING 
OF PyREX GLASS 


It is obvious that one has an almost endless variety 


of materials which may be analyzed. However, for the 
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purpose of analysis, they are characterized by the 
Prandtl Number at some reference temperature, and 
the manner in which the viscosity varies with temper- 
ature. Thus, one is still faced with an infinite set of 
parameters. However, if one is satisfied with knowl- 
edge concerning the effect of the film, then only three 
cases need be analyzed Prandtl Number much greater 
than unity, equal to unity, and much less than unity. 
The viscosity behavior of most liquids resembles an 
Arrhenius function and is, therefore, not a major 
determining factor. 

Since it is the purpose of this paper to determine the 
maximum effect of the molten film, a material which 
possesses a large Prandtl Number in the range of tem- 
peratures of interest and with no definite melting point 
was sought as the example. Due to the availability of 
data concerning its high-temperature properties, Pyrex 
glass (Corning 7740 glass) was chosen. Its thermal 
conductivity was estimated from the theory of normal 
liquids’ as k 1.71 X 10-3 B.t.u./fit. “F. sec. 

For the viscosity, it was found that the following ex- 


pression fitted published data’ quite well: 


uw = 0.0672e(8720/7T)'°!? Th. ft.-sec. (20) 


The density was taken as 131 Ib. ft.* and the specific 
heat as 0.29 B.t.u. Ib. In the case of Pyrex, evapor- 
ation or pyrolysis of its constituents represents the 
mass transfer. 

The equations were solved by means of a REAC 
differential analyzer. Preliminary results indicated 
that the terms F.* and FF;,- were negligible, and these 
terms were then omitted from the momentum equation. 

Two more parameters must be chosen. The density 
ratio p./ p was taken as 1.37 X 107‘; this corresponds 
precisely to Mach 18 at 90,000 ft. altitude. At this 


altitude, the density ratio does not vary appreciably 
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Fic. 3. Typical variation of surface temperature with non 
dimensional heat-transfer rate 


with Mach Number; hence the example may be con- 
sidered to be hypersonic flight at that altitude. The 
temperature in the interior of the glass, 7), is assumed 
to be zero. 7* has been taken as 4,000°F., but this 
does not affect the results; any value could have been 
chosen. 

The equations were solved for the following boundary 


conditions at the interface: 
$000, 3,S00, 3,600, 3,400, 3,200, 3,000°F. 


(uu u*)F-. Q). —5 x 10 —10 x 10 


1 


F, 65x 10“, 0) xX 


RESULTS 


A typical set of profiles are shown in Fig. 2. One 
immediately notices that all of the curves, except F;,, 
resemble exponential functions. This is to be expected 
from the nature of the equations. It is also interesting 
to note that the thickness of the molten flowing region 
is about one-tenth that of the heated region. Notice 
also that, at the point where F; Q, the value of 0, is 
only 70 per cent of its value at the interface. This 
means that only 70 per cent of the energy originally 
transferred to the interface is conducted into the non- 
flowing solid interior. Thus the presence of the liquid 
film cannot be ignored in this case. 

For the conditions investigated it was found that, 
as expected, the interface temperature increased mono 
tonically with interphase energy transfer, but varied 
slightly with interface shear stress and interphase mass 
transfer. A typical set of curves showing the interface 
temperature is shown in Fig. 3 in terms of the dimen- 
sionless variables. 

Now a high interface temperature has the effect of 
markedly decreasing the viscosity of the glass so that 
the liquid film flows more readily under the action of a 
given pressure gradient and shear stress. The decrease 
in viscosity thus causes the surface to recede at a faster 
rate, as shown in Fig. 4. Hence, as the heat transfer 
to the surface is increased, the rate at which the surface 
recedes is increased. This effect is shown in Fig. 5, 
where it is seen that the relation between heat-transfer 
and melting rate is almost linear. Note that the sur- 
face shear or interphase mass transfer have negligible 
effect on the melting rate for the range of values investi- 
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gated. It is not known whether this is merely coinci- 
dental for this particular material or if this effect is 
general. The maximum value of the tangential ve- 
locity component was about 2 X 107! u,. The non- 
dimensional film thickness was almost constant, 5.5 + 
10 per cent. The local tangential shear stress increases 
linearly with y from its value at the gas-liquid inter- 
face. Of course, this is meaningless beyond y,, for in 
that region the equations of elasticity apply, not the 
Navier-Stokes equations. 

Calculations were also performed for hypersonic 
flight speed at altitudes of 70,000, 120,000, and 160,000 
ft. At each of these altitudes, the relation between 
melting rate and heat-transfer rate was also nearly 
linear, and was also insensitive to the interface shear 
stress. The linear relation for each altitude was al- 
most identical to that shown in Fig. 5. 

The calculations were repeated for common lime 
glass, whose high temperature viscosity is much less 
than Pyrex. For this glass, it was found that the 
melting rates, as given by lim v(y), were five times 


v— 
greater than for Pyrex for the same energy flux to the 
interface. Thus, viscosity plays an extremely impor- 
tant role in determining the melting rate. 
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ASYMPTOTIC FORM OF ENERGY EQUATION 
For ¢ greater than ¢;, the value of F is virtually con 
stant. Hence, let 


lim F(c) =F (2] 


—> 


F(¢ > 5) 


The energy equation then becomes 
FO, + (1/P,*)0¢, 0 (22) 


This equation is easily solved for 0(¢) and 0,(¢) in terms 
of a multiplicative constant. Their ratios are given by 


6-(¢)/O(¢) = —F.P,* (¢> &) (23) 


Now, at ¢; it was found that 0(¢;) was equal to 0.62 + 


In addition, F.. is nearly proportional to 


Thus, the ratio 0;(¢5;)/0-. is almost 


3 per cent. 
0;, from Fig. 5. 
a constant. The value of the constant was 0.75 + 
0.05. Thus, the heat conducted into the unmelted 
solid is an almost constant fraction of that conducted 


into the molten film at the gas-liquid interface. 
DISCUSSION 
It is of interest to compare the results to the simple 
melting theory,' which may be expressed as follows: 


lim v(y) —[qi/Cppo(T; — To) ] (24) 


or, in terms of the dimensionless variables, 


lim F(¢) = —(T*0,./P,*T;) (25) 


os 


Eq. (24) cannot be applied to this example, because it 


does not predict the interface temperature 7. Even 
when the correct values of 7°; is inserted into Eq. (24), 
it predicts a value of lim F(¢) which is too small by 23 


3,000°F. and 28 per cent too small 


This is to be expected since most 


per cent for 7, 
for 7; = 4,000°F. 
of the material which flows has a temperature less 
than 7. If, on the other hand, one uses 7; instead of 
7; in Eq. (25), then the melting rate predicted by the 
simple theory is larger than that predicted by the 
present analysis. 

Eq. (19) shows that the thickness of the liquid film 


varies with (Ou, Ox) Since® 


(Ou,/Ox), = (va /R,) [2K(1 — K)]'’" (26) 


the thickness of liquid film depends on the Mach Num- 
ber, altitude, and size. For a melting test of a model 
in a rocket exhaust the predicted thickness of melt is 
about 0.016 in., which compares favorably with ex- 
periment; the corresponding value for a vehicle for 
which the Mach Number is 18 and R, l ft. is 0.071 
in. at 90,000 ft. altitude. 

In the above problem it has been assumed that the 
system is in a quasi-steady state so that all time de- 
rivatives are zero. In an actual vehicle the surface 
shear stress, pressure gradient, and heat flux may vary 
However, for a heating rate of 1,000 


(Continued on page 36) 
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Svmmetric Bodies in Ducts 


PHILIP LEVINE** 
Wright Air Development Center, ARDC 


SUMMARY 


The velocity potentials and stream functions for a point source, 
doublet, ring source, ring doublet, disc source, disc doublet, and 
ring vortex are derived for the case of axially symmetric incom- 
constant 
the 


pressible potential flow in an infinitely Jong duct of 


diameter. Several examples are carried out to illustrate 


nature of the flow about simple body shapes 


SYMBOLS 


a = radius of a ring or disc source, sink, or doublet 
é = natural logarithm base 
l = stream function 


= eigenvalue 
kind, re- 


J, Y = Bessel functions of the first and second 
spectively 

m = strength of a source or sink form 

VW = strength of a doublet form 

7 = dimensionless radial coordinate; ratio of radial dis 
tance to duct radius 

S = ratio of 2m/l 

l = superimposed uniform flow velocity 

= dimension'ess axial coordinate; ratio of axial dis 

tance to duct radius 

o = velocity potential 

¢,, ¢, = axial and radial velocity components, respectively 


= absolute value sign 


Subscripts 
’ = partial differentiation with respect to 7 


= partial differentiation with respect to 2 


INTRODUCTION 


/ PROBLEM of determining the complete flow 
field about axially symmetric bodies in ducts 
occurs widely. 
analysis of the performance of axisymmetric nozzles 
axisymmetric 


Typical instances are the design and 


and diffusers utilizing center-bodies, 
wind-tunnel arrangements, ring-type flame holders, 
and turbomachinery elements generally. In view of 
the numerous practical applications, it was decided to 
initiate a study of the problem for the case of incom- 


pressible flow. 
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Little theoretical work has been published concerning 
the flow about axisymmetric bodies which fill a sizable 
A solu 


tion for a point source in an infinitely long duct of 


portion of the cross-sectional area of a duct. 


constant diameter has recently been published. 
However, a numerical error occurs in the expression for 
the velocity potential of a point source. The expres- 
sion is rederived here, as it forms the basis for the de 
velopment of a general approach to the study of the 


flow about axisymmetric bodies in ducts. 


ANALYSIS 


The problem to be considered is the case of steady, 
incompressible and potential flow about axisymmetric 
bodies located on the axis of an infinitely long tube of 
constant diameter. 

The differential equation for the velocity potential 
of an axisymmetric incompressible flow is well known, 


and may be written as* 
o,, + (1 rid, + o.. = O (] 


Fortunately, this equation is amenable to solution by 
the method of separation of variables. Since Eq. (1 


is linear, solutions can be built up by superposition 


Point Source on the Tube Axis 


An expression for the velocity potential of a point 
source will be developed first. Once it is in hand, the 
solutions for a sink and doublet follow directly. 

Consider a source of strength m, as shown in Fig. 1, 
located ats = ry = 0. The boundary conditions which 


must be satisfied are: 


(I) g= + é, = GQ 
(11) z= +00; od: = 2m 
5= —O; od. = —2m 
(Til) v = Il; 6, = 0 
itV) 3 = @: @ =0 forr #0 


(For detailed discussions of the boundary conditions 
and other aspects of the following text, the reader 1s 
referred to references 1, 3, 4, and 5. 

Proceeding by the method of separation of variables, 


a solution to Eq. (1) exists m the form 


@ = 2mz + (Ax 
[BiJo¥jr) + BeYo(ir)] (2) 


To satisfy conditions (I) and (II), it is apparent that 


either A, or A» is zero (depending on whether = 1s 








34 JOURNAL OF THE AERONAUTIC. 


rei 

















r=| 


Fic. 1. Coordinate system and streamlines of a point source. 


positive or negative), and B, = 0. Condition (III) 
requires that 


Ji(j) = 0; 7#0 


Therefore, the values of (7) are the nth-ordered roots of 
J,, where J; is the Bessel function of the first kind of 
first order. To satisfy condition (IV), let 


@ = +=2mz + 


” 


Ane” 7" To( jar) (3) 
1 


Introducing condition (IV) and the orthogonality 
property of Bessel functions, one can solve for the 


coefficients, 


_9 “1 
A, = — | (¢:)- 


° J eB ot Jo(Jmr dr 
JIm- 0” (Im )J0 


Taking the size of the source to be very small (close to 
a mathematical point), and utilizing the Gauss theorem, 
the coefficients can be reduced to 


Am = (—2m)/jmJo?(jm) (4) 


Substituting Eq. (4) into Eq. (3) yields the velocity 
potential for a point source, 


—jniz 


@ = +2mz — 2m > Jo(jnr) (5) 


n=] Heda In) 


The stream function corresponding to the above velocity 


potential is 


e Jn\z 


y= +mr* + 2mr - * = Saar (6) 
n 1 InJ07(In) 
(+) forz > 0 


(—) fors < 0 


sign convention: 


At z = O, the stream function is discontinuous, so that 


approaching z = 0 from the left (see Fig. 1), Y = —m, 
while approaching z = 0 from the right, Y = +m. 


Doublet on the Tube Axis 
To derive the velocity potential for a doublet, one has 
only to place a source and sink of the same strength a 
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small distance apart on the axis, and find the resultant 
potential under the condition that 

M = lim m Ag 


Azc—O 


where .J/ is the strength of the doublet. The velocity 


potential for a doublet is then 
i 


@¢=2M+2Mu >> — Jsi» , 
1 Jo7(7 5 


The corresponding stream function is 


y = —2Mr © —— Jljv (S 
n=1 Jy*(J, 


Ring Source 

The velocity potential of a ring source can be found 
by using the same approach as taken for a point source. 
The boundary conditions put forth for the point source 
are valid in this case also, except condition (IV), which 
must be modified as follows: 


(IVa) = 0; @. = OV for rxa 


The total strength m of the ring source is taken to be 
uniformly distributed about the circumference, so that 
the source flow emitted by the ring is 47m. Therefore, 
Eq. (3) is valid for a ring source as well as a point source. 
To obtain the coefficients A,, condition (IVa) and the 
Gauss theorem must be utilized, finally yielding 


Am = ((—2) /jnIo2( jn) Jo jn) 


The velocity potential for a ring source is 


e ‘ . . 
@ = +2mz — 2m >> ——— Jo(jnr)SJo(jna) =(9) 
n=1JnJo"(Jn) 


The corresponding stream function is 


y 2 
y = tmr- + 2mr > : —— Jd gp Nigh ja) (6 
n=1JnJo°(J 
sign convention: (+) fors > 0 
(—) fors < 0 
Atz = 0, Y = +m, depending on whether the plane 
of the ring source is approached from the right or left. 


Ring Doublet 

The velocity potential for a ring doublet and the 
corresponding stream function can be obtained by the 
same technique as used for the point doublet. The 


resultant expressions are 


¢=2M+2uM>D : ——Ty(jyr)Jo(jna) (11) 
n=1. 0” (Tn) 


y= -2Mr > 


Jigar) Jo( jn) (12) 
1 JSo°( Jn) . ; 
Disc Source 


- . . - . . - 
f The boundary conditions for the velocity potential of 
a disc source are the same as for the point source, ex- 
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sultant , Which is now 


cept condition (IV 
IVb) zs= 0: d, = 0 for r,>a 


[he total strength m of the disc source is taken to be 
uniformly distributed over the disc, so that the source 
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flow emitted by the disc is tam. Eq. (3) is valid for 
this case also, and the coefficients A,, evaluated using 
7 condition (IVb) take the form 
‘ 
Am = [(—4m) Gjn?Jo?(jm) |\Ji(Jma 
rhe velocity potential and stream function are as fol- 
lows: Lai 
(S 
1 CO 
@ = +2mz2 - i —= Ih 77) ja) (12 
a y» BP ] (Jy 
lor os y Fic. 3. Doublets in a uniform flow 
ound y = +mr? 4 > — Ii(jur)STi(jnad (14) 
urce. a Cre sne? ; Pat ; 
se density of the distribution is equal to the strength of 
Ee . . . . - 
hich sign convention : +)torz>0 the vortex. Thus, by Eq. (15), one has 
—- for z < 0 
‘ = - > ¢ weil ° = 
Atz=0,y +m, depending on whether the plane Q 2ra*kh + 4raK 2, ; 124 Jo(Gnt) Si Gna ‘ 
; ; a ' ‘ 1 InJo-(J 
of the dise source is approached from the right or left. 
0 be ca ea init . 
hat Disc Doublet v= drake Ly acy Siar) Sil Jue IS 
, Following the approach for the point doublet, the 
¥ rh ' w! A is the strength of the vort 
‘e ° ° ° ° ° re S * Stren’? oO > VOrTtex. 
wee. velocity potential and corresponding stream function _ : S . _ 
the for a dise doublet are 
APPLICATIONS 
a if : ° . eqe . . 
@¢=2M+2M > Ty jun) Ii (jn 15 Yo facilitate the application of some of the results 
tInJv°(J obtained above, the velocity components and stream 
—4Mr Q function for a point source of one-half unit strength 
, = E i» y ( : ° » e 
, I =i Ino) /i( Jur) Ja(Fut 16 have been evaluated and tabulated in reference 3. 
(G) , . Jr - ; ‘ : Sakis 
5 The tables of J)(jr) and J,(jr) of reference 6 expedite 
Ring icaidiens the evaluation of the series for the other source-sink 
forms. 
The velocity potential for a ring vortex is readily 
(0) obtained by employing the property that a ring vortex Point Source in a Uniform Flow 
is equivalent to a uniform distribution of doublets over -_ ' , : ee 
: isbeeaues ; Che velocity potential of a point source in a uniform 
the surface bounded by it. The axes of the doublets ‘ eae : . a ; 
: : flow in a duct is found by adding the velocity potential 
are taken normal to the surface everywhere, and the ; ee : i a 
: Us for a uniform flow to Eq. (5). To obtain the cor 
ie responding stream function, one has to add the term 
Ur? 2 to Eq. (6). When dealing with semi-infinite 
bodies, a constant 2m must be subtracted from the 
stream function downstream of the source in order to 
maintain the usual convention that y 0 is the body 
e : 
. streamline. 
Bodies were calculated for several values of LU m, 
e ees 
and are shown in Fig. 2. Important features of the 
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Fic. 2. Point sources in a uniform flow 


body shape can be found quickly by employing the 
continuity equation. In the plane of the source, one 


has, by continuity, 


so that 


ro = 2m/l 


The radius of the body downstream at 2 = +© can 
also be found through continuity as 
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fag” = aro / tl -F To) 
rao? = (4m/U)/[1 + (2m/U)] 


It should be noted that there is a departure from the 
usual result that the velocity upstream at infinity is just 
that due to the uniform flow superimposed on the source 
flow. In this case the upstream velocity is less than 
the superimposed value by 2m. Similarly, the velocity 
downstream at infinity is greater by 2m. 


Doublet in Uniform Flow 


The velocity potential for a doublet in a uniform flow 
can be obtained by adding the velocity potential ls 
for a uniform flow to Eq. (7) and the corresponding 
stream function is found by adding Lr? 2 to Eq. (S). 

Bodies calculated for several values of .J/ LU’ are 


shown in Fig. 3. 


CONCLUDING DISCUSSION 


The examples presented above serve to illustrate the 
forms of bodies obtained by using sources or doublets 
in a uniform ducted flow. Application of the results 
obtained to complex body shapes is discussed in refer- 
ence 3. A preliminary attempt to incorporate the 
velocity distributions found from a body identical in 
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shape to that of the hub of an axial flow compressor 
indicates that these results can be used to estimate the 
inlet velocity profile, hub blockage effects, and other 


three-dimensional effects. 

In building up bodies of arbitrary shape, one can 
of course utilize any or all of the above source 
sink distributions through superposition. However, it 
should be pointed out that the results above can easily 
be extended to include the case of varying duct wall 
contour, by using source-sink forms with a = 1 [see 
Eqs. (9) and (10)]. 
diffusers, and ducted bodies can be studied. 


Thus, many forms of nozzles 
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Hydrodynamics and Heat Conduction of a Melting Surface 


(Continued from page 32 


B.t.u./ ft.*-sec., simple melting theory indicates that 
close to steady-state melting is achieved one to two 
seconds after the start of melting; this has been verified 
experimentally. Hence the method can be used for 
the case of slowly varying external conditions. 

In order to apply the results to an actual case, mj, 
u;, 1, g:, and 7; must be matched with the boundary- 
layer solution. However, calculations have been per- 
formed which employed estimates of the aerodynamic 
heat-transfer rate. The results indicate that the melt- 
ing rate is sufficiently small so that this method of 
surface protection may be practical for certain applica- 


tions. 


CONCLUSIONS 


For the case considered, it appears necessary to in- 
clude the effect of the liquid film in order to accurately 


9 


determine the rate at which a surface recedes due to 
melting from aerodynamic heating. Additional solu- 
tions are required to determine the general behavior 
of melting films in order to solve the transient case, 
which is the one of greatest technical interest. 
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The Crippling Strength of Compression 
Elements 


GEORGE GERARD* 
New York Unwersity 


DUMMARY 


his paper consists of a review, primarily for the air-frame 
designer and stress analyst, of a new method of crippling strength 
inalysis. Available test data on various aluminum, magnesium, 
ind titanium alloys and steels are presented in terms of general 
ized parameters, based on the overall geometric configuration 
f the element and two mechanical properties of the material 
the elastic modulus and compressive yield strength. The method 
applies to both formed and extruded elements and _ stiffened 
panels 

The predictions of the method are checked by available test 
data on short-time crippling at elevated temperatures. Simi 
larly, the method is extended to the creep crippling problem by a 
suitable definition of the time-dependent vield strength under 
creep conditions 

Finally, the efficiencies of various materials in crippling appli 
cations under elevated temperature conditions are reviewed and 
simple test methods are proposed for rapidly determining the 
crippling strength of new materials under various environmental 


conditions 


SYMBOLS 


Wrought Alloy Designation* 


Alloy Former New 
\luminum 145 2014 
R301 2014 
17S 2017 
248 2024 
52S 5052 
61S 6061 
75S 7075 
Magnesium FS-1, FS1 \Z31B 
0-1, 01 AZSOA 
Temper Designations* 
\lloy Former New Conditioning 
\luminum -0 0 annealed 
-T T3 heat treated 
-T -T4 extruded and heat treated 
-T -T6 heat treated and aged 
-RT T36 heat treated and rolled 
1/4+H H32 one-quarter hard 
Magnesium h -H24 hard rolled 
HTA -T5 extruded and aged 
* Notre: The alloy and temper designations used herein refer 


to the designation in use at the time the test data were obtained 


Notation 
A = area of element; area of stiffener plus sheet of width, 
b,, in.? 
b = total width of stiffened panel, in. 
b. = width of sheet between stiffeners, in 
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neering. 


= web width, in 
E = modulus of elasticity, psi 
= cladding thickness as a fraction of total plate thick 
ness; form factor 
= number of cuts plus flanges 


m = exponent 

N = compressive loading, lbs. /in 
R = radius, in 

t = thickness, in 

t; = time for creep crippling, hours 
t] = function of time 


weight per unit length, Ibs. /in 


ll 


a@ = crippling coefficient 
8 = crippling coefficient for two corner elements 
8 = generalized crippling coefficient based on ¢ 
” = cladding reduction factor 
n = |1 + 3(¢,,/¢,,)f 1 + 3f 
Ga = applied creep stress, ksi 
o = cladding vield strength, ksi 
o = buckling stress, ksi 
¢., = compressive yield strength, ksi 
a; = crippling strength, ksi 
¢ = density, pci 


Subscripts and Superscripts 
= skin 
= stiffener web 
= effective or weighted value 


INTRODUCTION 


Siw INTRODUCTION and proper selection of many 
new materials in air-frame design under room tem- 
perature conditions and both short time and creep con- 
ditions under elevated temperatures has placed an 
urgent requirement on the acquisition of strength 
data. One of the most basic of these data in air-frame 
analysis and design is crippling strength. Conse 
quently, it appeared that a fruitful approach would be 
a comprehensive review of available test data to ar- 
rive at a semiempirical method of sufficient generality 
and simplicity that it could be applied with reasonable 
confidence to the prediction of the crippling strength 
of materials, both old and new, under various condi- 
tions of temperature and time. 

The purpose of this paper is to review a new method 
of crippling strength analysis which has resulted from 
an extensive investigation of room temperature crip- 
pling data on flanges, plates, angles, Z, channels and 
multicorner sections and stiffened panels, both formed 
and extruded.''*}* The review herein consists pri- 
marily of a presentation of the test data in terms of 
certain generalized parameters which reflect the over- 
all geometric configuration of the element and only 
two mechanical properties of the material, & and ay. 
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The aim is to present the method and data for use of 
the designer and stress analyst and therefore the de- 
tails of the development of the method are omitted. 
In cases where the development is of interest, refer to 
references | and 2. 

Most crippling methods in current use involve sub- 
division of the cross section of the formed section into 
a series of plates’ or angle sections.® The crippling 
strength of the section is obtained as a weighted 
average of the crippling strengths of the subdivided 
elements. The method described herein is empirical 
in nature, as are the other methods, but does not re- 
quire the subdivision of the cross section into elements 
for analysis. The only geometric parameter re- 
quired for the analysis is the developed length-thickness 
ratio, A f*. Ina sense, the method is a generalization 
of the work of Schuette,® Gallaher,’ and Needham’ on 


formed sections. 
(1) ROOM TEMPERATURE 


THEORY AND EXPERIMENTS ON FLANGES AND PLATES 


It is a well known phenomenon that the failing 
strength of a flat plate can appreciably exceed the buck- 
ling stress when the latter occurs in the elastic range. 
This behavior results from the boundary constraints 
at the unloaded edges of the plate which permit signifi- 
cant tension stresses to develop in the plane of the 
plate after buckling. 


Crippling of Flanges 


A hinged flange is the simplest element, from the 
standpoint of analysis, which can carry loads con- 
siderably in excess of the elastic buckling load. As 
compared to a column which is unsupported along the 
unloaded edges and thus deflects appreciably at 
buckling, the flange is supported along one unloaded 
edge. This boundary constraint tends to stiffen the 
flange during postbuckling rotation and permits post- 
critical loads to be carried. 

Because of the nonlinearity of the postcritical be- 
havior and the additional nonlinearities of plasticity, 
failure analyses of buckled elements are complex. For- 
tunately, Stowell® has succeeded in conducting a failure 
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TABLE | 
Test Data on Extruded Angle and T-Sections 


1 /o 
Ref Shape A /t? Gey Ts/o vt- \ E ) 
Angle Sections 
y LS 208 21 $3 ksi 0.754 0.67 
4 LS 209 me 13 0. 737 0.70 
4 LS 210 22 13 0.764 0.70 
) 6 15. 1 12.4 0.45 1.44 
T- See , 
| LS 236 38.5 13 0.721 0.82 
H LS 247 50 13 0. 680 1.06 
| LS 248 50 13 0.652 1.06 
} LS 249 39 $3 0.750 0.83 
9 7 19 6 1H). 4 0.65 1.09 
In the extrusion fillets for angles, 0 R/t 1/2; tor T- 


sections, 1/2< R/t< 1 
® Values of &;/¢-, correspond to an L'/p = 20 


analysis of a flange with a simply supported, straight 
unloaded edge. These results are of considerable im- 
portance because they supply quantitative data as to 
the failure mechanism. Furthermore, the analysis 
provides a reference frame from which semiempirical 
methods can be devised to treat other cases of failure 
for which complete analyses do not exist. 

The very significant physical fact which is brought 
out by this analysis is that the edge stress is intimately 
associated with failure. Apparently failure occurs 
when the stress intensity at the edge reaches a value 
approximately equal to the compressive yield strength. 
This is indicated by Stowell’s analysis where it is shown 
that, for eight widely different hinged flanges, the edge 
stress intensity is a constant within one per cent of the 
compressive yield strength. 

The results of the analysis which are important from 
a design standpoint are shown in Fig. 1. It can be 
observed that Stowell’s theoretical curve is in good 
agreement with test data on 245-T4 aluminum alloy 
cruciform extrusions. In such sections, the opposed 
flanges constrain the simply supported unloaded edges 
to remain straight. For equal flange angles, however, 
the simply supported unloaded edges warp in the plane 
of the flange after buckling. 
data for extruded angles given in Table 1 lie below 


Consequently, the test 


Stowell’s theoretical line which pertains to an undis- 
torted unloaded edge. The test data for extruded T- 
sections given in Table | and shown in Fig. 1 follow 
the trend line for the cruciform sections. 

The choice of parameters used for the coordinates 
of Fig. 1 is given in some detail in references 1 and 2. 
For purposes of completeness however, a brief develop- 
ment follows. Schuette® has shown quite convincingly 
that the results of crippling tests on Z, channel, and H 


extrusions may be represented in the following form: 
Oy Ocr = AlOy CO )’ 1) 


In fact, Stowell utilized the stress parameters of Eq. 


(1), & oe and o.,/o-, in plotting the results of his 
theory. 


By using the relation for the critical stress 
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Ge = KE (t/b)? 2) 
Eq. (1) can be converted into 
&,/o-, = aK'~"[(t/b) (E/a,,) 2 PO 7" (3) 
In a simplified form 
S)/Ocy = Bl (t/b) (E/o,,)'/?]” (4) 


Eq. (4) contains the parameters used in Fig. 1 and gives 

the crippling strength directly in terms of the geometric 

properties and mechanical properties of the section. 
From Fig. 1, the following values are obtained for 


flange elements for &;/ a-y < Q.S: 


unloaded edge m 3 
straight 0.40 0.67 5) 
. = Bae (. 
distorted O.S5 0.56 


The empirical values for the straight unloaded edge 
flange fit the theory and test data reasonably well. 
The values selected for the distorted edge case are 
based on the few test data for angles as well as con- 
siderable other evidence presented in the subsequent 
sections. 

For flanges that cripple in the inelastic range é, o., > 
0.8, test data® conclusively indicate that buckling and 


failure are essentially coincident. Therefore, 


oS; = oc¢; for &,/o4>0.8 (6) 


Crippling of Plates 


As indicated previously, the failure analysis of 
buckled plates is mathematically complex. 
variational principles in conjunction with calculations 


Mayers and 


By use of 


performed on an electronic computer, 
Budiansky” were able to obtain results related to fail- 
ure of buckled plates. 

The load-carrying capacities of flat plates with un- 
loaded edges constrained to remain straight were com- 
puted for 24S-T3 plates that buckle elastically at pre- 
scribed fractions of the compressive yield strength. 
Although the average compressive stress did not have 
a maximum value in the range of end shortening con- 
sidered, the curves were very flat from a strain of ap- 
proximately 0.006 on. Thus, the average compressive 
stress at a strain of 0.010 was taken as an indication of 
failure. 

The results of this analysis, together with experi- 
mental data on failure of flat plates, are shown in Fig. 2. 
A very significant fact brought out by the theoretical 
analysis is that the load-carrying ability of plates with 
straight unloaded edges may be significantly higher 
than plates with unloaded edges which are free to warp. 

In discussing the failure of flat plates and flanges, it 
is necessary to identify differences in boundary condi- 
tions which by themselves may not cause a difference 
in buckling stress: (a) the unloaded edges do not 
warp and remain straight in the plane of the plate (v = 
const., Fig. 3a), and (b) the unloaded edges are free to 
warp in the plane of the plate (0, = 0, Fig. 3b). 

Plates tested in V-groove supporting fixtures are free 
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to warp. In addition, out-of-plane displacements at 
the unloaded edges are not entirely prevented after 
buckling due to lateral shortening in the central region 
of the plates, as shown in Fig. 3. Similarly, the post 
buckling behavior of compressed square tubes closely 
follows that of plates tested in V-groove supports 
For the tubes, both warping and out-of-plane displace 
ments of the corners are evident after failure. 

The parameters used in Fig. 2 are the same as thos¢ 
used in Fig. 1 with the exception that the g term is 
introduced in the abscissa coordinate. As il'ustrated 
in Fig. 4, g represents the number of flanges and cuts 
required to reduce the section to a series of flange ele 
ments. This scheme is quite arbitrary and has the 
virtue that results on many different sections can be 
reduced to one crippling equation by its use, as will be 
shown herein. 

The results for angle sections when presented in terms 


of 6 ¢ are already in the desired form, since, for the 
flange, 
b/t = A/2t 1 gi ‘ 
Thus, Eq. (4) becomes 
Ty Ocy = Bl (gt? A) (E/ o-,)'/?]" g 


The values for m and 6 given in Eq. (5) now apply for 
the flange for the g-correlation scheme. 

For the V-groove plates and square tubes, g 3 
and 12, respectively, as shown in Fig. 4. The results 
shown in Fig. 2 indicate the following values for plate 


elements: 


unloaded edges m B, 


straight 0.40 0.65 


distorted 0.85 0.56 


It can be observed that these values for plate elements 
are in substantially good agreement with the values 
given by Eq. (5) for flange elements. 

Crippling of H Sections 


Rather extensive test data on H extrusions of four 


aluminum alloys and one magnesium alloy® '* are 
available to test further the g-correlation scheme. As 
shown in Fig. 4+, g = 7 for an H extrusion. The test 


data correlated according to the parameter A_ g/* are 
illustrated in Fig. 5. It can be observed that satis 
factory correlation with Eq. (S) is obtained by utilizing 
the 2, (5) straight unloaded 
edges, for G, oy < O.S. 
yond this cutoff Eq. (6) applies. 

The value of m = 0.40 obtained for the H sections is 
the same as that for flanges and flat plates with straight 
Evidently, the opposed flanges of the 


value for cruciform, Eq. 
As discussed previously, be 


unloaded edges. 
H section prevent any warping of the flanges at the 
corner. 

From the results presented in Figs. 1, 2, and 5, it 
appears that the g-correlation scheme permits an 1n- 


tegrated approach to the crippling strength analysis of 
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FF IG. 3. Stresses and displacements of flat plates after buck 
ling under conditions of uniform end shortening (reference 11) 

i) Straight unloaded edges, (b) stress free unloaded edges free 
to warp in the plane of the plate 
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CRIPPLING STRENGTH OF 


the extruded and machined flanges and plate elements 


considered. 


FORMED AND EXTRUDED SECTIONS 


rhe crippling strength of formed and extruded sec- 
tions 1s associated with the short column region (L’  p 
pproximately less than 20) where the strength does not 


vary significantly with the length of the stiffener. The 
various factors influencing crippling of flange and plate 
Additional 


factors such as the increased stress-strain characteristics 


elements also pertain to formed sections. 


in the corner region as a result of forming influence the 
crippling strength of formed sections. Also, sections 
are often formed of Alclad sheet and it is necessary to 
consider the effect of the cladding upon the crippling 


strength. 


Formed Angles 


Needham’ has conducted tests on crippling of equal 
flange angles formed in the -T condition to an inside 
radius approximately equal to 3/, and the test data 
are shown in Fig. 6 according to the parameters of Eq. 
(8S). The empirical constants for Eq. (S) which fit the 


data within +10 per cent limits are 


material m 8 
248-T3 O.S5 0.66 
Alclad 24S-T3 0.85 0.66 (10) 
Alelad 75S-T6 O.S5 0.63 


Che 8, values differ somewhat due to different strain- 
hardening characteristics in the formed corners; an 
effect discussed in a subsequent section. 


Two Corner Sections 


For two corner sections such as Z and channel sec- 
tions, the buckling coefficient is a function of b; } 
Consequently, it was found! that the development 
leading to Eqs. (4) and (S) does not hold for two corner 
sections where the flange elements apparently influence 
significantly the crippling strength. The analysis of 
a large mass of test data presented in reference 1 indi- 
cate that the following form of Eq. (8) applies for two 


: (=) - (11) 
A \oe, 


7 correlated ac- 


corner sections: 
C7/ Cc = 


The test data for several materials® 
cording to Eq. (11) are shown in Fig. 7. The following 
values of the crippling coefficient, 8, fit these data 


within +5 per cent limits: 


material 3 
17S-T 1.00 
248-T +. 00 (12?) 
FS-1lh 3.48 


The lower value of 8 for the FS-1h magnesium alloy is 
attributed to lower strain hardening characteristics in 
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Fic. 8. Crippling data on formed Z and channel sections 
references 4, 5 a) Alclad sections, (b) unequal flanged Z 


and J sections 


the formed corners, as discussed subsequently. Beyond 
the cutoffs shown in Fig. 7, Eq. (6) applies. 

In reference | a mass of test data on extruded Z and 
channel sections with filleted corners of R ¢ | corre- 
lated according to Eq. (11) with 8 = 3.64. The corner 
fillet apparently infiuences the crippling strength some- 
what by delaying the extension of the buckles through 


the corners which occurs at failure. 


ABLE 2 
Test Data on Two-Corner Elements 
(ee) 
Ref Shape t A /t? og To t? \nFE 
5) El 0.0252 136.7 42ksi 0.376 22.10 
in 
E2 0.02538 87.2 42 0.538 14.10 
E3 0.0504 67.6 {2 0. SSS 10.92 
E4 0.0255 95.5 42 0.500 15.45 
} X785-C 0.025 105 $] 0.478 16.52 
0.040 64.5 }] 0. 639 10.68 
0.064 410.3 4] 0.920 6.93 
0.025 103 51 0.431 16.67 
0.040 64.5 51 0.573 11.47 
0.064 $0.3 5 0.853 7.45 
} XLS-165 0.025 99.4 41 0.505 15.94 
0.040 62.1 }] 0.648 10.28 
0.025 99.4 5] 0. 502 17.12 
0.040 62.1 51 0.680 11.02 
j LS-160 0.025 SO.8 41 0.578 12.97 
0.040 50.5 }] 0.790 & 36 
O.OR25 SO.8 5 0.507 13.92 
0.040 50.5 51 0.704 S.US 
4 LS-161 0.025 93.6 41 0.525 15.03 
0.040 58.5 4] 0.715 9 68 
0.025 93.6 5l 0,455 16.15 
0.040 58.5 35l 0.596 10,40 
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(reference 16). 


Alclad Correction 

The cladding which may be present on formed sec- 
tions acts to reduce the value of the buckling stress and 
therefore its principal effect is reflected in the cladding 
reduction factor, 7. Needham? has presented a con- 
siderable amount of test data on Alclad channel sec- 
tions and thus the data of Fig. 7 on 24S-T Z and 
channel sections can be used as a standard of compari- 
son with Needham’s Alclad 24S-T3 data to substantiate 
the method of correcting for the cladding. 

The cladding correction factor 7 can be computed 
from the following formula, given in reference |: 


= [1 + 3(¢a/oer)f]/( + 3f) (13) 


3 


where f is the ratio of total cladding thickness to total 
thickness (f = 0.10 for Alclad 24S-T3 and f = 0.08 for 
Alclad 75S-T6). 

Needham’s test data are presented in the upper por- 
tion of Fig. 8 and compared with the value of 8 = 4.00 
obtained from the 24S-T data of Fig. 7. It can be 
seen that the data fall within the +10 per cent scatter 
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TABLE 3 
Variation of 8 With @,,/o-, Effect 


Material R/t a 8 3/3.2 Soures 
Z and Channel Sections 
Machined / 
MST 3A1-5Cr 0 1.00 3.20 1.00 Fig. 9 
RC-130B 
\ll extrusions 1 fillet 1.00 3.64 1.14 Ref. | 
Ale. 75S-T6 3 1.05 (est 3.42 1.07 Fig. 8 
FS-1h 1.3 1.15 3.78 1.18 Fig. 7 
17S-T 3 LZ $00 1.25 Fig. 7 
24S-T 3 1.35 $00 1.25 Fig. 7 
Ale. 24S-T3 3 1.35 (est 41.00 1.25 Fig. 8 
RC-130A 3 >. 40 Fig. 9 
C-110M } 3. 54 Fig. 19 
RC-55 } 3.95 Fig. 9 
RC-70 } 3.95 Fig. 9 
Angles R/t Tey / Oo Bg 8,/0.56 
248S-T extrusions (0-1/2 1.00 0.56 1.00 Fig. | 

fillet 

248-T3 3 1.35 (est 0.66 1.18 Fig. 6 
Alc. 24S-T3 3 1.35 (est 0.66 1.18 Fig. 6 
Ale. 75S-T6 3 0.638 Fig. 6 
C-110M } 0.63 Fig. 18 


+10 per cent limits (Figs. 6 and 8) the method of 
correcting for the Alclad coating appears to be satis 
factory. 

The Alclad 75S-T6 data in Fig. 8 fall consistently 
below the Alclad 24S-T3 data. These data are fitted 
within +10 per cent limits by Eq. (11) with a value of 
6 = 3842. The lower strain hardening characteristics 
of Alclad 75S-T6 appear to be responsible for the lower 
value of £. 

In the lower portion of Fig. 8, test data of Table 2 
on unequal flange channels’ and J sections* are pre 
sented in terms of the parameters of Eq. (11). It can 
be observed that the available data correlate within 
+1() per cent limits with the line for equal flange Z and 
channel sections of Fig. 7. Although this relatively 
simple method of analysis requires substantiation by 
considerably more test data than presented, the avail- 
able test data indicate that two corner elements tend 
to behave in the same manner regardless of their shape. 


Increased Corner Properties 

The increased corner properties of sections which are 
formed in the -T condition have been discussed briefly at 
various points. Now that the values of 8 in Eq. (11) 
have been established for formed Z and channel sections 
of several different materials, it is possible to examine 
quantitatively the effects of the increased corner prop- 
erties on crippling. 

The average compressive yield strength, é.,, of the 
corners of Z and channel sections as compared to the 
unformed sheet has been investigated by Heimerl and 
Roy!*® for 17S-T, Heimerl and Woods! * for 24S-T 
and Alclad 75S-T, and Gallaher’ for FS-lh. The ratios 
of the average compressive yield strength of the corner 
&-,, for these materials are presented in Table 3, to- 
gether with the value of 8 obtained herein on formed 


and extruded sections. 
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Also given in Table 3 are some recent test data on 
titanium alloy channel sections as presented in Fig. 9. 
These data are particularly valuable since they include 
machined channel sections with no fillet radius in the 
corner, Which then serve as a lower bound for the 
corner efiect. 

It can be observed that all of the 8 values for two 
corner sections given in Table 3 lie between 3.2 and 
1.0. The data are plotted as a function of é., o-, in 
Fig. 10 and indicate a fairly consistent trend. By using 
3 3.2 as a lower limit, a correction factor for the cor- 
ner effect is also given in Fig. 10. The increased 
properties in the formed corners can apparently raise 
the value of 6 for two corner sections by as much as 
25 per cent. 

Also given in Table 3 are the available data for 
angles. The data are too few to provide other than the 
estimate for the corner effect given in Fig. 10. 

The raised yield effect in the corners of formed sec- 
tions can result in significant increases in crippling 
strength. If this effect is considered to be of impor- 
tance in design, a systematic test program should be 
conducted to investigate the influence of forming and 
heat treatment upon the 6 value. 

It has been observed that the corner fillet in extru 
sions can result in an increase in 8. This also re- 


quires a systematic investigation. 


Multicorner Elements 


For convenience, all sections with more than two 
corners are designated as multicorner elements. Avail- 
able test data of Needham’ on multicorner sections of 
Alelad 24S-T3 formed in the T-condition, and of 
Crockett? on formed Alclad 24S-T and Alclad 24S-RT, 


TABLE 4 
Test Data on Multicorner Sections 


A (: 
Ref Shape o G; a ot? sf 


AC y 
y \l $2 ksi 0.624 0.913 S 
AZ 12 0.545 1.028 
F3 }2 0.839 0.615 
Gl $2 0.624 0.916 
G2 12 0.424 1.445 
5 Ci 42 0.412 1.478 1] 
ber 4 12 0.574 0.950 
C3 }? 0.655 0.760 
C4 42 0.748 0.664 
C5 12 0.571 0.940 
Fl 12 0.569 0.942 
F2 42 0.567 0. 980 
D1 }2 0.414 1.428 
D2 $2 0.671 0.798 
D3 $2 0.555 0.961 
D4 12 0.452 1.188 

3) D5 42 0.659 0.793 

} X 785-D $] 0.537 0.975 

4 X 785-D 51 0.760 0.639 

} X 785-D | 0.477 1.088 

} X 785-D 51 0.710 0. 7i2 

3) H 2 0.595 0.928 14 
} LS-152 51 0.775 0.678 

4 LS-153 51 0. 864 0.55 

5 Bl 42 0.495 1.128 17 
5) B2 42 0.645 0.806 


‘ 
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formed multicorner sections (references 4, 5 


The available test data are shown 
0.85, fit the data well 


are given in Table 4. 
in Fig. 11. Lines of slope, m 
according to Eq. (S) and the value of 8 0.55 is in 
close agreement with that obtained for angles, V-groove 
plates and square tubes. The data are only for the 
24S-T material and therefore no conclusions concerning 


the raised yield effects in the corner can be drawn. 


STIFFENED PANELS 


In analyzing the crippling strength of stiffened panels 
which consist of formed Z and hat or extruded Y sec 
tions attached to a flat sheet, certain additional con 
siderations arise beyond those encountered in the an 
alysis of the stiffener section alone. To avoid con 
sideration of the method of attachment of sheet to 
stiffeners, it is assumed that the panel is monolithic. 
Rivet strength requirements to achieve this condition 


are presented in reference 2. 


Z-Stiffened Panels 

Test data on the crippling strength of flat stiffened 
panels with Z-section stiffeners were obtained from three 
different sets of data listed in Table 5. In the first 
set, there was a systematic variation in the parameter 
t, t.. In the second set, nine different materials were 
used covering a wide range of / and o,, values. The 
third group covers potential strength estimates of the 
highest obtainable crippling strength of a large group 
of riveted panels. The potential strength of a riveted 
panel corresponds essentially to the strength of a mono 


lithic panel. In all tests the slenderness ratio of the 


panel L’ p = 20. 

The cross-sectional area of a stiffened panel is taken 
as the area of the stiffener plus the area of sheet corre 
sponding to the stiffener spacing. To account for the 
t, ts variation in the most simple and direct manner, 
the parameter A gf* for the formed element is replaced 
by A ¢t,,/, for the stiffened panel. Again this procedure 


arbitrary. It is relatively simple to 


is somewhat 
introduce some weighting factors which may result in a 
more realistic evaluation of the A gf* parameter for 
stiffened panels. However, such factors lead to addi- 
tional complexity and do not improve the correlation 


sufficiently to warrant their use. 
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CRIPPLING STRENGTHS OF 


TABLE 5 
Crippling Strength Test Data on Stiffened Panels 


Ref Material Stiffener f, Variation Remarks 
17 24S8-T Zz 1.00, 0.79, 0.63, strongly riveted 
0.51 panels 
lS 795-T6 Z 1.00 data for panels 
61S-T6 identified as 25 
525-1/4H 50-20 and 37.5 
75S-O 75-20 used since 
SAE 1010 other panels indi 
Copper cated evidence 
FS-1lh of wrinkling 
18-8-3 /4H 
Ti-1/4H 
248-T3 Zz 1.00, 0.63 potentialstrengt! 
798-T6 estimates 
20 248-T hat 1.25, 1.00, 0.63, h range of 0.6 
0.39 to 1.2 
21 248-T Y 1.00, 0.63, 0.40 
798-T 


The available data on the monolithic crippling 
strength of formed Z-stiffened panels are shown in 
Fig. 12 and represent the data from the three sources 
listed in Table 5 according tot, f, values. Satisfactory 
correlation is obtained with Eq. (S) in all cases using 
7.55 and 8, = 0.56. This correlation holds up to 


x 


values of &;/0,, = 1. 


Hat-Stiffened Panels 

The crippling strength analysis of formed hat- 
stiifened panels is essentially the same as that used for 
the Z-stiffened panels. The test data of reference 20 
cover a range of b, b, values (b, = width of top web, 
b,, = height of side webs) between 0.6 and 1.2 and a 
range of ¢,/t, as given in Table 5. 

rhe test data for each value of ¢, ¢, are shown in 
Fig. 13. All data correlate well with Eq. (S), although 
there is a small systematic variation of the 8, value with 
the parameter /,f, that may be due to the closed form 
of the hat section. 

For the particular riveting used for the hat-stiffened 
panels there is some evidence of wrinkling at é; o., 
values greater than 0.8, as shown by the dotted lines 
in Fig. 13. This may account for the somewhat anoma- 
lous behavior of the data shown in which the @, value 
for the b,/t, = 75 group is somewhat higher than for 
the b, t, group of 25, 35, and 50. 
picion of wrinkling, it is recommended that a cutoff of 
5;/o-, = O.S be placed on Eq. (8) for hat-stiffened 


Because of the sus- 


panels and that caution be exercised in using this 
equation for ¢,,/f; = 0.39 panels. 


Y-Stiffened Panels 


In analyzing the crippling strength of monolithic 
stiffened panels utilizing extruded Y stringers, certain 
additional considerations are encountered which have 
not arisen previously. For the particular Y extrusion 
used in reference 21, the thickness of the horizontally 
opposed flanges at the top of the section as shown in 
Fig. 14 is greater than the other elements in the cross 
section. Consequently, in determining the effective 
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thickness of the stiffener 7,,, the following method of 
weighting was used: 


i = Th t vh, 14) 


In Eq. (14), 6; and ¢; refer to the length and thickness 
respectively, of the cross-sectional elements 
Similarly, for the 75S-T panels, the compressive yield 
strength of the stiffener (c,, = 7S ksi) was significantly 
higher than that of the sheet (¢,, 67 ksi). The 
following weighting procedure was used 
Fey = SOcy. + Sey, ((b/ts) — 1) (it 15 
where / is the effective thickness of the stiffened panel 
Asa consequence of these weighting pr cedures, the 
generalized crippling relation for the Y-stiffened panel 
has the form 


G/F = Bo\(elt,/A) (E/Be,)*!?]" 16 


The test data of reference 21 for the 24S-T and 75S-T 
panels listed in Table 5 correlated according to the 
parameters of Eq. (16) for three different /,, f, ratios 
as shown in Fig. 14. As for the hat-stiffened panels 
the 3, values vary somewhat with the /,, ¢, parameter 

In view of the desirability of having one method of 
crippling analysis as far as possible, the results of the 
previous sections of this paper are summarized for all 
elements analyzed according to Eq. (Ss). All 3, values 
are given in Table 6, from which it can be observed 
that they are in substantially good agreement at /,, / 
1.0. The variation of 3, for hat- and Y-stiffened panels 
with /,, /, values other than unity are summarized in Fig 
15 in terms of a 8, correction. It can be observed that 
the correction is substantially the same for the hat 


and Y-stiffened panels. 


TABLE 6 
Summary of Empirical Values of 8, for Flange and Plate Elements 


and Stiffened Panels 


Element 


Disto ted Un loaded cd es,m =0 95 


1) extruded angles 2 0.56 
2) \ groove plates 3 0.56 
(3) extruded tubes 12 0. et 
$) formed multicorner sections 8-17 0.55 
5) formed Z panels 7.83 0.56 
6) formed hat panels 16.83 
1.25 0.59 
1.00 0.56 
0). 63 0 50 
0.39 0.48 
7) extruded Y panels 18.83 bw/t 
1.16 0. 56 
0.732 0.51 
0. 464 0.48 
Straight Unloaded Edges, m = 0.40 
(8) straight edge plates (theory 3 0.65 
9) extruded T 3 0.67 
10) extruded cruciform } 0.67 
11) extruded H  § 0.67 


For @;/oc¢y < O.8 
» For b./t. < 50 
Tentatively, for lack of sufficient test data, for ;/a-y < 
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perature. Because of thermal flight considerations, 


it is of interest to determine whether the method of an 


alysis can be applied to predicting the crippling strength 


of various sections under elevated temperature condi 


tions. At such temperatures, time is an important 
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short-time crippling strength as well as crippling under 
creep conditions, or “creep crippling.” 


Flanges and Plates 


A small amount of test data on short-time crippling 
of H extrusions at elevated temperatures were pre- 
sented by Heimerl and Roberts.** These data are 
shown in Fig. 16 in terms of the parameters of Eq. (8). 
Also shown is the line representing the room temper- 
ature behavior obtained from Fig. 5 and it can be ob- 
served that there is no significant diterence in the ele- 
vated temperature behavior of the H extrusion in this 
case. 

Some recent test data on the short-time crippling 
strength of 2024-T3 V-groove plates at elevated tem- 
peratures were presented by Mathauser and Deveikis.** 
These data are shown in Fig. 17 in terms of the param- 
eters of Eq. (8) and it is important to observe that 
there is no significant difference between the elevated 
temperature and room temperature data. 

However, when the value of 8, = 0.56 obtained 
from Fig. 2 on other V-groove plates is used, a con- 
siderable number of test points fall below the — 10 per 
cent limit. This is particularly true of the room tem- 
perature test points which should check with 3, = 0.56 
since this empirical value was established from other 
V-groove plate tests on 2024-T3. This suggests pos- 
sibly some defect in the test jig or technique. 

A value of 8, = 0.52 agrees well with the data and all 
points fall within the 10 per cent limits. Aside from a 
possible discrepancy resulting from the test technique, 
it may be concluded that there is no significant differ- 
ence between the room and elevated temperature test 


data. 


Formed Sections 

Short-time crippling data on formed titanium alloy 
C-110M angles and channel sections are shown in Figs. 
1S and 19 in terms of the parameters of Eqs. (S) and 
(11), respectively. These data were obtained by the 
Boeing Airplane Company and are given in reference 16. 

It can be noted that the room and elevated temper- 
ature data follow the same trend. The 8, and 8 values 
obtained from Figs. 1S and 19 are listed in Table 3 and 
it can be observed that the values for C-110M lie 
within the limits of the 8, and 8 value obtained for other 
materials. It is not possible to correlate these data 
according to the raised yield effect in the formed corners 
because of the absence of &,, ¢,, data for these sec- 
tions. ‘These sections were stress relreved after forming 
and therefore the raised yield effect is probably small 
in this case. 
Stiffened Panels 

Mathauser and Deveikis*! have recently obtained 
some data on the short-time crippling strength of Z- 
stiffened panels of ¢,, ¢, = 1.0 at elevated temperatures. 
These data are shown in Fig. 20 and again it can be ob- 
served that there is no significant difference between 


the room and elevated temperature data as correlated 
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according to the parameters of Eq. (8S). Also shown in 
Fig. 20 is the room temperature value of 8, = 0.56 ob- 


tained from Fig. 12. 


CREEP CRIPPLING 


The creep crippling phenomenon is concerned with 
the time required for failure or collapse under a con- 
stant creep stress of the various compression elements 
considered herein. It is distinct from creep buckling 
in the same sense that short-time crippling is different 
than short-time buckling when the latter occurs elas- 
tically. 

It would be desirable if the methods of crippling 
analysis developed for room temperature, which appear 
to be applicable under short-time elevated temperature 
This 
requires a time-dependent strength property of the 


conditions, could be applied to creep crippling. 


material which plays the same role in creep crippling 
that o,, does in short-time crippling. Furthermore, in 
creep crippling, the applied constant creep stress, &,, re- 
places the short-time crippling strength ¢, and the de- 
pendent variable is now the time for failure, f,. 


Method of Analysis 


In Part (1), Stowell’s analysis of the crippling of 


flanges indicate that failure occurs when the stress in- 
tensity reaches o,, or a permanent strain of 0.002 in. in. 
This apparently holds for various types of sections on 
the basis of the correlation obtained herein using o 
Therefore, in the creep crippling analysis, use will be 
made of this behavior to define a_ time-dependent 
strength property to replace o,,. 

An empirical equation which has been used to repre- 
sent the creep behavior of 2024-T3 aluminum alloy is 
of the following form :** 

«= 0 E+ Ae*t* (17) 
The terms A, B, and A are empirical constants at a 
given temperature. Denoting the creep stress corre- 
sponding to a permanent strain of 0.002 in. in. as 


g,,|t], as suggested by Mathauser and Deveikis,’ 
and noting that, under this condition, 
e — 0 /E = 0.002 (1S) 
Fq. (17) can be transformed into 
o-,{t] = (1B) In (0.002, A2*) (19) 
For creep crippling, Eq. (S) becomes 
Gy Oey (t] = By} (gt?/A) (E/ oe, {t})/24" (20) 


In Eq. (20), &, is the applied creep stress and a,,[t] is 
obtained from Eq. (19) or from any other creep equa- 
tion at a permanent strain of 0.002 in. in. at the time 
for creep crippling, f;. In a similar manner, Eq. (11) 
for two corner sections can be derived. 


V-Groove Plates 
Mathauser and Deveikis** have conducted tests on 


creep crippling of 2024-T3 V-groove plates. Their 
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test data include é, and the observed ¢, for individual 
plates as well as the pertinent values for the constants 
A, B, E, and K in Eq. (17). The data are shown in 
Fig. 21 in terms of the parameters of Eq. (20). 

With the exception of the numbered points on the 
left (the number indicates /,;), there is no significant 
difference among the creep buckling data and the room 
temperature test points. The several 500°F. points 
which are low are for crippling times of 9 hours or 
more. It is not known whether this discrepancy is due 
to Eq. (20) or possible inaccuracies in Eq. (17) at 500°F. 
for times greater than 9 hours. 

The value of 8, = 0.52 in Fig. 21 is directly com- 
parable with that obtained in the short-time crippling 
tests shown in Fig. 21 using the same test procedures. 
Thus, the remarks concerning the discrepancy between 
the room temperature data of Fig. 21 as compared to 
Fig. 2 pertain here also. 


Stiffened Panels 


In Fig. 22, some data on creep crippling of Z-stiffened 
panels tested by Mathauser and Deveikis** are shown 
according to the parameters of Eq. (20). Again it can 
be observed that there is no significant difference be- 
tween the room temperature and creep crippling data. 
Also shown is the room temperature value of 6, = 0.56 
obtained from Fig. 12. 

(3) SIGNIFICANCE OF RESULTS IN DESIGN 

After having analyzed a relatively large mass of data 
on various shapes and materials under various environ- 
mental conditions, it appears pertinent to summarize 
the essential results and view them from the broader 
aspects of the air-frame structure rather than as indi- 
vidual elements. 

The values of the crippling coefficients given in the 
following formulas do not include raised yield effects 
in the corners of formed sections or fillet effects in ex- 
trusion corners. The crippling coefficients are the 
lowest 8 and 8, room temperature values presented in 
Tables 3 and 6. For the sections with distorted un- 
loaded edges considered herein 
plates, multicorner sections and stiffened panels of 


angles, tubes, V-groove 


w/t; = 1-—the following crippling equation applies 
within +10 per cent limits: 
&;/Ocy = 0.56[(gt?/A) (E/oey)*!?}°-5 (21) 


For the sections with straight unloaded edges considered 


herein—plates, T, cruciform, and H_ sections—the 
following applies within +5 per cent limits for ¢;/o,, < 


O.S: 
Gy Oc, = 0.67 [(gt?/A) (E/o-,)'/?]?- (22) 


For two corner sections—Z, J, and channel sections 
the following applies within +10 per cent limits up to 
the values of ¢,/0,, shown in Figs. 7, 8, and 9: 


&// ey = 3-2[(t2/A) (E/oe,) 39-78 (23) 
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Beyond the cutoff values of &, ¢,,, plastic buckling 
generally occurs, and, therefore, 


Or = Cer (24) 


Corrections for Alclad materials are given by Eq 
(13) and 6-corrections for the raised yield effect in the 
corners of two corner sections are given in Fig. 10. For 
multicorner sections and stiffened panels, the raised 
yield effect does not appear to be of any significance 
in the data analyzed herein. For the /, ¢, effect on 
Y- and hat-stiffened panels, 8, correction factors are 
given in Fig. 15. In creep crippling problems, o,,, is 
the creep stress at a permanent strain of 0.002 in./in. 


Crippling in Air-Frame Structures 


In an air-frame structure, various stiffener shapes 
are invariably attached to skins, and, therefore, crip- 
pling strength estimates such as given by Eq. (23) for 
two-corner sections are of little practical value. The 
really significant crippling strength estimate is that 
of the stiffened panel or V-groove plate which is repre- 
sented by Eq. (21). This difference in the behavior of 
a two-corner section when tested alone as compared 
to when used in a stiffened panel, points up the fact 
that far too great an emphasis has been placed on 
crippling tests of Z and channel sections as representa- 
tive of the crippling behavior of the air-frame structure. 

Multicorner sections, square tubes, V-groove plates, 
and stiffened panels of f, f, = 1 eripple according to 
Eq. (21). Thus, it would appear that use of a multi- 
corner section, such as a hat section or a V-groove 
plate, to determine the crippling strength of a par- 
ticular material would be much more representative of 
the air-frame structure than use of a two-corner section. 

In addition, test data on the former elements indicate 
little, if any, &:,/o-, eect. In a practical sense, it is 
generally not possible to take advantage of this effect 
in design since stiffened sections are either extruded or 
formed in the annealed or as-quenched condition prior 
to aging, or else hot formed. These techniques tend 
to wash out the &,,/0,, effect as compared to cold form- 
ing in the fully heat-treated condition. 

From these considerations, it would appear that, in 
the design sense, only Eq. (21) is of significance. Since 
it applies to crippling of V-groove plates and stiffened 
plates as well as multicorner sections, tests on V-groove 
plates or simple hat sections would appear to have 
application to a wide variety of structural applications 
at room temperature as well as for short-time and creep 
crippling at elevated temperatures. Since the method 
is empirical in nature, the value of 8, may possibly vary 
somewhat from the value of 0.56 for materials and 
environmental conditions other than those considered 
herein. Eq. (21) offers the possibility that room tem- 
perature data may be extended to various materials 
and elevated temperature environments and thus re- 
duce the exhaustive tests associated with new materials 


and environments. 
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Comparative Crippling Efficiencies 

The remarks of the previous section focus particular 
attention upon Eq. (21) for structural design applica- 
tions. In such applications, the comparative crippling 
efficiencies of various configurations and materials par- 
ticularly at elevated temperatures are of interest. This 
concluding section, therefore, is devoted to a brief study 
of crippling efficiencies based on an extension of the 
methods of references 25 and 26 which were devoted 
to crippling of plates and multicorner sections. 

For efficiency studies, it is necessary to transform 
Eq. (21) so that it contains the structural index .V 6 
applicable to the crippling strength of V-groove plates 
and stiffened panels. By use of the following relation- 


ships: 


Eq. (21) can be transformed into 


[f(N/b) (E Tey) ol cia 26) 


Gr = (OgGc,)° 


g(t t.) (f./2)2 (b/6,) 27) 


where f = 


The coefficient fis a form factor where, for a V-groove 
plate, f = g. It can be observed from Eq. (26) that 
the form factor represents the efficiency of the cross- 
sectional configuration of stiffened panels in crippling 
applications. 

The weight per unit length of a stiffened panel or plate 
is simply 


w= obi 28) 
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Comparative crippling efficiencies of various material classes at elevated temperatures for short times 
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3y use of Eq. (25), Eq. (26) may be written as 


wt b? g Nb) &,]| 29) 
Upon substituting Eq. (26) into Eq. (2%), the following 
weight equation is obtained 


f/b? = [(fN/b)/B,] ¢/ E***) (E/e 1 (30) 


e 
“a 


If the allowable strength is limited to o,,, then the fol 


lowing weight equation applies 


w/b? = f(N/b) (¢/oe-,); for & 2 ta 31) 


It can be noted that the weight equations contain 
the loading index, .V }, and the form factor, /, which 
are loading and geometric properties. The pertinent 
material properties are ¢, /, and o,, for crippling 
strength [see Eq. (30)], and ¢ and o,, for compressive 
strength |see Eq. (31)]. 

By employing the procedures presented in reference 
26 for optimum stiffened panels, a wide variety of 
magnesium, aluminum, and titanium alloys and pre 
cipitation hardened stainless steels have been investi 
gated for crippling efficiencies. By considering the 
best material in its temperature range for a given alloy 
class, results are obtained which are representative of 
the best currently available properties. As discussed 
in reference 26, an idealized stress strain relation 1s 
used in this analysis for which Eq. (50) applies 1n the 


elastic range and Eq. (31) in the yield region. 
The results of this analysis are shown in Fig. 25, from 


which it can be seen that the aluminum alloys have a 
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Potential crippling efficiencies of the new high-strength heat-treatable titanium alloy of the Department of Defense sheet rolling program 


for short times at elevated temperatures. 
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The Characteristics of Two-Dimensional Sails 


in Hypersonic Flow 


WALTER DASKIN*® ann LEWIS FELDMAN** 
Gruen Applied Science Laboratories, Inc. 


SUMMARY 


The aerodynamic and geometric characteristics of two-dimen- 
sional sails in a hypersonic flow are investigated in this paper 
These surfaces would have very large lift and drag per unit 
Therefore, they appear attractive wherever large values 
W or CpA 
expressions are derived for the sail shape and the corresponding 
Charts of these quantities are 


weight 
of the parameters C,A W are needed. Simple closed 
lift, drag, and moment coefficients 
presented for the entire useful design range. The results can also 


be applied to two-dimensional parachutes. 


SYMBOLS 


A = sail, area, sq.ft 

B = total length of sail, ft. 

C = chord of sail, ft. 

Cy = drag coefiicient 

C, = lift coefficient 

Cy = moment coefficient 

D = drag, lbs. 

L = hiift, Ibs. 

\/ = pitching moment, ft.lbs 

7. = free-stream Mach Number 


= pressure, lbs. per sq-ft. 
= stagnation pressure behind normal shock, lbs. per sq.ft 
= dynamic pressure, (1/2)pV?, lbs per sq.ft. 


S = length along sail, ft 

7) = tension per unit span, lbs. per ft. 

V velocity, ft. per sec. 

Hi” = weight, Ibs. 

\ = longitudinal coordinate (Fig. 1), ft 

y = normal coordinate (Fig. 1), ft. 

a = angle of attack, rad., degree 

6 = local angle of sail, rad., degree 

vy = arbitrary family designation (Eq. 13), degree 
p = density, slugs per cu.ft. 


Subscripts 
7 = trailing edge 
L = leading edge 
= conditions far upstream of sail 


INTRODUCTION 


[ IS DESIRABLE, in a number of supersonic and hyper- 
sonic problems, to construct surfaces having large 
values of lift and or drag per unit weight. The possi- 
bility of using a sail of woven wire cloth seems attrac- 
tive since such a device would have extremely small 
weight per unit area and could presumably be made of a 
material that could withstand the high temperatures to 
which it would be exposed in hypersonic flight. The 
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aerodynamic characteristics of such a two-dimensional 
sail in a hypersonic flow are investigated in this paper. 

The problem originally arose as a simplification of the 
incompressible sail problem which will be reported else- 
where by Lewis Feldman. 


ANALYSIS 


Consider the coordinate system of Fig. 1, with uni 


form flow in the direction of the x axis. If we have a 


high ratio flexible sheet of uniform chord 


stretched with very high tension in the direction normal 


aspect 


to the paper (or if this sheet is stiff in the direction nor- 
mal to the paper but flexible in the plane of the paper), 
If p 


is the local pressure difference across the flexible sheet, 


we reduce the problem to a two-dimensional one. 


and 6 the local angle of inclination, then 


pcos@ds = Td(sin@) or p= T(d6 dS (1) 


where 7’ is the tension per unit span and dS is the dif- 
ferential arc length of the sail. 

Taking the modified Newtonian approximation! for 
the pressure on the lower surfaces, 


cos” @ (2a) 


Piower = Po sin? 6 + p 


In order to simplify the algebra and to make the re- 
sulting sail shapes independent of Mach Number, it is 
convenient to set Pupper = Px COs? @. 

This gives the required zero pressure difference at 
= 0 and it is a good approximation when 6 + 0 because 
the upper surface pressure is much smaller than po, and 


Po << po. Then the local pressure difference across the 
sail is 
p = posin? 6 (2b)* 
Equating (1) and (2b) we have 
po sin? 6 = T(dé dS (3) 


the basic equation of equilibrium of sail elements under 
tension and hypersonic loading. 

This is integrated to 
unmodified Newtonian been 


* Had the pressure equation 
)» ) 


used instead at Eq. (2) or (3), the normalization quantity for 
7 instead 


of po /T, and the factor 2 would have appeared in the aerodynamic 


lengths in all subsequent equations would have been 2g 


coefficients in place of po/g. Inclusion of an added term to cor 


rect for the finite-density ratio across the shock would have 


added no further difficulty. Use of the ‘Newtonian plus centrif 


ugal’’ form would have resulted in an _ integro-differential 


equation 
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Fic. 1. Hlustration of sail nomenclature 
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Fic. 2. Shape of hypersonic sails. 
poS/T = | (dé sin’ 6) = cot 6, — cot 6 (4) 
Jl 


where S is the are length of sail measured from the lead 
ing edge and 9, is the leading-edge angle. 
The total length of the sail, B, is given by 


(poB T) = cot 6, — cot 47 (5 


where 67 1s the trailing-edge angle. 
The xv and y coordinates are, respectively, 


pox T = [(1/sin 6,) — (1 /siné)]| | 

poy/T = In [tan(6/2) /tan (6,/2)]) ” 
giving the coordinates of the trailing edge as 

pouxr T = (1 sin 6,) — (1 sin 67) | a 

poy (T In [tan(@, 2) tan(6@, 2) |\ - 
The chord, C, of the sail is then 

poC/T V (poxr T)? + (poyr/T): (S 
and the geometrical angle of attack, a, is given by 

a = tan™! [(povr /T)/(porr/T) | (9) 


The shape of any such sail can be completely deduced 
from Eq. (6). This universal shape is plotted in Fig. 2. 
Superimposed on the actual sail shape are its local angles 
of inclination, 6. One can therefore determine the shape 
of any sail, given its leading- and trailing-edge angles, 
by simply using that portion of the universal curve be- 
tween these two angles. The chord may be drawn be- 
tween these two points, and the geometrical angle of 
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attack measured directly. The complete shape for sails 
with trailing-edge angles greater than 90° may be de 

termined by reflecting the curve of Fig. 2 about a ling 
parallel to the x axis passing through the point 6 = 90 

Thus, one could conceivably find the shape of slack 
two-dimensional sails which are secured through a fair 
lead or even two-dimensional parachutes, although it is 
doubtful whether the approximation of Eq. (3) is valid 
for such a shape if it is not ventilated near the region oj 
its vertical slope. 

One can also use Fig. 2 to determine the leading- and 
tailing-edge angles of a sail with fixed geometrical ang 
of attack and various values of tension by moving a ling 
representing the chord at the given angle of attack 
closer to the curve as the tension increases, until the 
curve and chord are tangent at the condition of infinite 
tension. 

Although the origin of the figure is at the leading edge 
of a 10° sail, it should be noted that this is quite arbi 
trary, and the origin could be shifted to any other point 
on the curve. Since, with the aerodynamic force nor 
mal to the surface the tension is uniform, we can find 
the total lift, drag, and leading-edge moment by simply 
considering the attaching forces and their angles at the 
leading and trailing edges. The results, identical with 
those found by integration of the pressure, are 


L = T(sin 67 — sin 6 } 
D = T (cos 6; — cos 07 (10 
M= T'(v7 COS 67 — X7 sin 07 \ 


and the corresponding coeflicients are 


Cc = L/(1/2)pV2C = (po/q)[(sin 67 — 

sin 0,) (poC T) 
Cp = D/(1/2 pV°C = (po g) X 

[(cos 0, — cos Or) (poC T 
Cy = M/(1/2)pV2C? = (po/a) {I(porr/T) X {| 
cos 07 — (poxr T) sin Or} (poC T)*; 

C, Cp = LD = [(sin 67 — sin 6; 

(cos 6; — cos 6; 


It should be noted that for the infinite tension (flat 


plate) case, these expressions reduce to 


c = (po q) sin? 6 cos 6, Cp = (fo q) sin® 6 


L/D = cet? f 


Cy = —(po’q)(sin? 6/2), ae 

The value of Pp») g is virtually independent of Mach 
Number for hypersonic flow, and therefore, as expected 
the coefficients are independent of Mach Number. The 
value Pp g = 1.839 corresponding to infinite Mach 
Number, and a specific heat ratio of 1.4 was used in all 
of the calculations. It can be shown that the results do 
not vary significantly for 5/3 > y > 1.2 

Table 1 lists the value of py) g as a_ function 
of the free-stream Mach Number, .1/.., for a perfect gas 
with specific heat ratio 1.4. It is seen that the asymp- 
totic value is essentially reached for Mach Numbers in 
excess of 4.* 


*The values in Table 1 may be used as Mach Number cor 


rections for the coefficients plotted in Figs. 3 and 4. 
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Fic. 3. Lift and drag characteristics of hypersonic sails 


One may choose, as independent variables, either @ 
and @; or any other pair of independent variables, such 
as chord pyoC 7 and angle of attack a. For convenience, 
the leading- and trailing-edge angles were chosen as the 


independent variables in the calculations. 


RESULTS OF CALCULATIONS 


and fC 7, 


the lift, drag, and moment coefficients and the lift-drag 


The dimensionless lengths pox 7, poy C, 


ratio were computed for the entire range of values of 6 
and @7 which might be encountered. 6, was taken from 
10° to 90° in 10° increments, and 67 was taken in 10 
increments from 6, to its complement. 

The results for C, and Cp are best represented by the 
lift-drag polar, Fig. 3. The envelope of the curve repre- 
sents the flat plate (infinite tension sail). Two families 
of curves appear in the figure, those of constant leading 
edge angle and those of constant angle v, where p is de- 
fined by 67 = 1SO° — 6; — ». 

The figure is considerably less confusing when plotted 
in terms of vy rather than 67, since the families plotted in 
this manner tend to intersect each other at relatively 
large angles, while the 67 method of plotting results 
in curves which are tangent, or nearly so, all along the 
envelope. 

Fig. 3 shows the value of »v which maximizes the lift 
for each leading-edge angle. The angles on the envelope 
of the curve represent the slope of the infinite tension 
sail at each point on the envelope. As can be seen from 
the figure, values of L D in excess of unity can be 
achieved for relatively tight sails with small leading- 
edge angles, even when moderate friction drag coeffi- 
cients are added to the drag coefficients here computed. 

Fig. 4 shows the aerodynamic characteristics of flat 
plates (infinite tension sails) as computed by Eqs. 

12). These values do not, of course, include the usual 
corrections due to the hypersonic boundary layer. 
Such corrections would be most important at low angles 


IN HVPERSONIC FLOW 











1 Lf for 1.4 
Ss 1. 850 1.84 S42 S S30 
) 1.006 1. OO4 )] 
| F 
3} 1.5} 
L/D Cp 
2} 1.0F 
~ | a 
= | lo, | 
li 
I} .5} NCL 
\ 
\ 
ol ol 1 4 4 4 1 1 4 
Oo 10 20 30 40 50 60 70 80 90 
@—s 
Fic. 4 \erodynamic characteristics of hypersoni ils wit 


infinite tension 


of attack and would result in an appreciable decrease 
of the computed values of L D. Nevertheless, Fig. 4 is 
included as representing the characteristics of the ele 
mental portions of the sails used in the calculations 
Values of IV .1 have been computed for several sail 
configurations using woven-wire cloth or uniform sheet 
material under tension. The configurations investi 
gated included swept-back delta as well as two-dimen 
sional high and low aspect r itio sails. The weight esti 
mates included the weight of spars necessary to prop 
erly support the sail. These sail weights were compared 
with the estimated weights of conventional wing struc 
tures, using accepted wing-weight estimating tech 
For identical materials, stress levels, and plan 
10 of the weights 


niques. 
forms, the sail weights were | 5 to | 
of the conventional wing structures. This significant 
decrease in weight is primarily due to the fact that the 
sail structure carries its load primarily in tension, while 
the conventional wing must be designed for bending 
loads. It should be noted that the weight of the sail 
itself was negligible in most of the cases investigated, 
the bulk of the weight of the sail structure being in the 
spars. 

Since the lift and drag coefficients of sails are roughly 
equal to those which can be expected from conventional 
surfaces at hypersonic speeds, see Fig. 3) and the sail 
weights per unit area are significantly less than conven 
tional Il’ .! values, the use of sails will result in large 
values of the parameters C,.1 IV and Cp.! HW. Large 
values of these parameters are often required in hyper 


sonic flight. 
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On Optimum Nose Curves for Missiles in the 
Superaerodynamic Regime 


H.S. Tan 

Professor of Aeronautical Engineering, 
Detroit, Mich. 

July 15, 1957 


University of Detroit 


IT REFERENCE I, it was shown that for a missile with a smooth 
surface in the superaerodynamic regime, the consideration 
of specular molecular reflection results in the following differential 
equation for the determination of optimum (minimum drag 


shape of the meridian nose curve 
y” = —(y’2/y) (1 + y’2)/(k — yy”? l 
where k = 3 for regular reflection and & = 2 for random reflection 
It was pointed out that there is little hope of finding an analytic 
solution to Eq. (1), and so, numerical integration was resorted to 
It is interesting to note, however, that a first integral of Eq 
1) can readily be obtained as follows: y'dy'/dy, 


Eq. (1 


Noting vy” = 
can be written 

ydy'/dy = —y'(1 + y’2)/(k — y”? 
for example, 


—dy/y = |k/y’ — (kR+ 1)y’/1 + y’?)] dy’ 2 


Now Eq. (2) can be integrated once to give 


yy"*®/A + y’2)@t02 = C 3 


where Cis the integration constant 
The-analytie behavior of the optimum nose curve can then be 


) 


qualitatively discussed, based on Eq. (3), as follows 

If the curve ever crosses the axis y = 0, or if its slope vanishes 
anywhere, then C must be identically zero. This results in the 
trivial case of 


y = 0, or constant (4 


which is physically evident, but has no practical interest 
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If C does not vanish, then the curve cannot reach the axis y = 0) 


but can be made to approach the axis as closely as we like 
From Eq. (1), it is evident that the curve has a cusp at y 


Vl Now writing Eq. (3) in the form 


y= C(l + y"?)¢ afy" 5 
nd considering y’ as a parameter, we can differentiate Eq. (5 
with respect to y’. To obtain an extremum value of y, we must 


solve the equation 
dy/dy’ = 0 6 


Substitution of Eq. (5) into Eq. (6) yields the solution 


y = tvi i 


Thus it is seen that at y’ = Vk, y has, indeed, an extremum 


Moreover, at y’ = Vk, the curve branches into two separate 


portions—one for which y’ < /k, and the other for which y 


\ } 
rhe asymptotic shapes of these two branches of the curve can 


further be obtained. For the branch with y Vk, when 
1, we have 
dy/dx y C/y 8 
Upon integration, this gives 
= Cy" Fe GC, GQ = k/(R+ 10C 9 
For the branch with y’ Vk, when y’ >1, we have 
dy/dx = y’ = 9/C 10 
Upon integration, this gives 
v Ciny+ C 1] 
C:, C; are both integration constants 
Since the differential Eq 1) is of second order, two conditions 


in be imposed on an initial point of the curve By choosing 


y = 1 12 


C is uniquely specified, and a numerical integration can be 
readily carried out using Eq. (3 

Phe two branches of curve have been qualitativ ely indicated 
in Fig. 1 It is interesting to note that the branch of the curve 


shown in reference 1 for which y’ > Vk is incorrect. But since 
this branch has no physical significance itself, validity of the 
result is not affected. On the other hand, the resulting curve for 
Vk obtained by numerical integration should be checked 


for its asymptotic behavior—i.e., Eq. (9 
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On Optimum Nose Shapes for Missiles in the 
Superaerodynamic Region 


|.D. Chang 

Graduate Student, Department of Aeronautics, Calitorr 
of Technology, Pasadena, Calif 

July 15, 1957 


a Institute 


N A RECENT PAPER in this JOURNAL,' W. J. Carter made a study 
I of the optimum nose shapes for missiles in the superaerody 
namic region. He found that for the case of specular molecular 
reflection from a smooth surface, the optimum nose shape is given 


by the equation 
y”" = —(y"*/y) [Cy + 18 — y” (1 


where y is the radius of the nose about the longitudinal axis 


We want to point out that Eq. (1) can be readily integrated in 


term of simple functions 


Multiplying both sides of Eq. (1) by dy/dx, we get 
p(dp/dx) = —(p?/y p? + 1)/(38 — p?)| (dy/dx 
which may be written as 
dy y p? — 3)dp/p(p* + 1 ; tf p> + | — (3/p {dp 
integrating, 
y ( 1)*/f ( const 2) 
\lso, from dx dy/p, using the result of Eq. (2 
dx ( — 2p? -3 
Hence 
x log + p /4 + } 
Eqs. (2) and (3) give the solution for Eq. (1) in parametric 
form For the case calculated by Carter, 2.50 and ¢ 
—(0.235. Eqs ind (3) then give the values for x and y shown 


} 


\s shown in Fig. 1, the curve y(x) has two branches, separated 
by a cusp at which y = 3 For the case calculated here, this 
| ; \ 
point corresponds to a 2.18 and y 7.70 
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Fic. 1 Nose curve for minimum drag with specular reflection 
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t 


The numerical results shown here differ considerably from that 


given by Carter \lso, in Fig. 3 in his paper, the nose curve 


corresponding to y V 5 Is erroneous 
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Supersonic Flow Around Blunt Bodies 


Hyman Serbin 

Physical Scientist, Aeronautics Department, The RAND Corporation 
Santa Monica, Calif 

June 7, 1957 


Kh NEWTONIAN THEORY Of impact has been shown to be 
useful for pressure calculations on the forward facing part 
of bodies moving at high speed. It is now a familiar practice 
to use this information to calculate nonviscous velocities at the 
wall and then to estimate rates of heat transfer. This pro 
cedure is perhaps open to question; heat-transfer rates depend 
on velocity gradients which are not given by the Newtonian 
analysis. Nor can one obtain information on boundary-layer 
stability or all the body stability derivatives. It seems, there- 
fore, inevitable that, as design proceeds with these hypersonic 
missiles, there will be a greater need for more accurate aerody 
namic theories either to predict what will happen in unfamiliar 
flight conditions or to effect an extrapolation from a known test 
result to the design condition 

\ start in this direction has been made by the writer! ? for 
two kinds of blunt bodies, one with a flat-dise face and the other 


with a hemispherical face. Our theory depends on the limiting 
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Fic. 1 (Top). Detachment distance versus density ratio A for 
discs. Fic. 2 (Bottom). Detachment distance versus density 
ratio A for spheres 
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Fic. 3. Comparison of experimental shock wave pattern 
V/ 


= 6.7 with theory 


case where y> the ratio of specific he tts, approac hes unity witl 
out being equal to 1. The limit 1 characterizes the New 
tonian theory; hence our theory may be characterized as a re 
formulation of the Newtonian theory on the basis of gas dynamics 

Following the scheme used by Prandtl for the simplification of 
the hydrodynamic equations when the Reynolds Number is 
large, we have simplified the gas-dynamie equations in the limit 
where .J/, the free-stream Mach Number, is large, and y —@ | 
The results obtained inelude the shock detachment distanee 
shock shape, streamlines, and velocities. These are to be re 
garded as asymptotic values appropriate to the limits described 
above 

For the purposes of designing at high speed, the results seem 
to be very satisfactory However, since there is substanti 
amount of data for blunt bodies at low supersonic speeds, it is of 
interest to see how our results can be extended into this speed 
range 

There seems to be an impression among some aerodynamicists 
that the case y ¥ 1, J < can be correlated by replacing bot! 


parameters by the single parameter 
K = density ratio across normal shocks at free-stream 1/. 


We do not know the origin of this remark; however, it may be 
easily derived from our equations! for large A, at least for the 
subsonic part of the flow field 

To translate our results to finite Mach Number, one need only 
make the formal substitution 


—1)/2—-1/K | 


In particular, the shock detachment distance A takes the form 
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Fic. 4. Comparison 


READE 


A/R = 1.03 A 2 dis¢ 2 
A/R = (2/3)K sphere 3 


R represents the radius of the dise and sphere respec 


\nother analysis of the flow around blunt bodies has been 





Probstein takes occasion to rewrite our asymptotic formula for 
jises in a similar fashion to that shown above but with A’ re 
placed by A — 1 

Within the framework of our theory (A large), there is no way 
) distinguish A from A — 1. The use of either one in Eqs. (2 
nd (3) is equally satisfactory when one examines test data at 
ige K. However, the use of A — 1 is attractive from a formal 
standpoint because when the free-stream Mach Number ap 
prow hes unity, the distance A so expressed approaches infinity 
This is consistent with our knowledge of transonic flow 

Thus, our formal rule is to make the substitution everywhere 


in references land 2 


n which the stand-off formulas become 





A/R = 1.03 (K — 1)~"? disc rs) 
= 2/3 (K — 1 sphere 6 


Plots of these relations are shown in Figs. 1 and 2, respectively 


ind compared with test data. Fig. 1 also shows theoretical re 


sults drawn from Probstein. The agreement between test and 
ur theory is very satisfactory at high A (= high 1/ What is 


surprising, and hardly to be anticipated, is the good agreement 
in the low supersonic range 

\ similar comparison could be made between the theoretically 
lerived and experimentally-observed shock. They appear to be 
imost identical at the higher Mach Numbers (Figs. 3 and 4 
The lower Mach numbers have not been checked 

Quite apart from the evidence presented, it appears that our 
theory is accurate (to principal order in A because all the 

juations defining the flow field are utilized. The shock shape ts 

consequence” of the analysis, not part of an assumption. To 
mention one item of interest, the velocity is constant in magni 
tude along each streamline almost up to the sonic line The 
direction of flow, however, changes continuously so that the 
flow is not the shear flow assumed by some writers. 

We have dwelt on the validity of the theory principally be 
cause of one theoretical consequence which is somewhat unex 
pected \ceording to our theory, the velocity along the stream 
line at the wall would be zero almost up to the shoulder of the 
yoy This means that the inviscid flow near the wall resembles 


Couette flow 


This is quite different from what one would 


of experimental shock wave pattern at 
VJ = 5.8 with theory. 


uted recently by Probstein \t the end of his paper, 
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expect on the basis of slender-body aerodynamics Phere the 
inviscid theory shows a finite velocity at the wall which must 
be reduced to zero by interposing a viscous boundary layer. On 
the other hand, our theoretical flow would exhibit a temperaturs 
discontinuity which can conceivably be bridged by a kind of 
boundary layer 

It might be argued that a velocity discontinuity would ay 
in a second-order analysis, and this is indeed the case Yet it 
would be worth while to see if the available heat-transfer data 
could not be satisfactoril\ predicted on the basis of the first-order 


inviscid shear flow already derived 
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Marin’s Strength Hypothesis* 


A. D. Topping 

Engineer Specialist (Stre 
Akron, Ohio 

June 17, 1957 


A HYPOTHESIS for describing the limiting conditions of stress 


in materials having different properties in different dire« 
tions and different strengths in tension and compression was re 
cently proposed by Marin It is the purpose of this note to 
show that the hypothesis is not only questionable, but untenabl 
Marin begins with the distortion energy or octahedral shear 
theory of strength, according to either of which the stress-space 
surface describing failure is a circular cylinder whose axis makes 
equal angles with each of the principal stress axes and is the axis 
of hydrostatic stress. Its equation may be written either as 
Marin’s Eq. (2) or his Eq. (7 Phis theory has been successful 
in predicting over a fairly wide range yield strengths of isotropic 


ductile metals having equal tensile and compressive strengths 


Marin speculates that it can be made to apply to all materials 
by generalizing Eq. (2) to the form of his Eq. (12)——-i« 
“n= <6 r w= & om 4 + 
q o —_ s o = T 7 —_ o ( 7 = ‘ x 
Oo 7 


where oj, 02, and o; are principal stresses 

Eq. (12) contains five constants, a, », c, g, and o, to be deter 
mined experimentally If one considers now orthotropic ma- 
terial having different tensile and compressive strengths along 
each of its three natural axes, it is apparent that any equation 
purporting to describe the limiting stress conditions for such a 
material must have at least six independent constants 

The possibility remains that Eq. (12) may prove useful for some 
less general case. Upon further examination, it can be seen that 
in a rectangular coordinate system with the coordinates o;, 2, ¢ 
Eq. (12) describes one of the following surfaces of revolution de 


pending on the value of g: 


g < — 1, hyperboloid of one sheet 
= —1, circular cylinder 
—-1< gq < 2, ellipsoid of revolution 
q = 2, two planes perpendicular to hydrostatic axis 


2, hyperboloid of two sheets 


Consider now the triaxial strength of an isotropic material 
having unequal tensile and compressive strengths, a problem to 


* The writer is indebted to his colleague, Mr. J. D. Marketos, for identi 


fication of the geometric surfaces described by Marin's equation 
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which Marin's ‘‘generalized strength theory’’ should apply as to 
a mere specific case 

Marin implies that his hypothesis is applicable to both yield 
and ultimate strengths. This is impossible because we know 
from experience that vielding is simply not the same phenomenon 
as rupture; rupture cannot always be described as a function of 
stress alone. Yielding is a shearing or sliding phenomenon while 
rupture may be either a sliding failure or a cleavage or separation 
The application of a hydrostatic tension may produce 
Thus it is conceivable that rupture 


failure 
rupture, but not vielding 
occurs before yielding. The nature of the yield surface for states 
of stress in this region is only of academic interest, of course, and 
is quite incapable of being determined experimentally. The 
limiting stress surface for rupture is closed at that end of the 


hydrostatic stress axis. At the compression end of this axis, 


however, we find experimentally neither vielding nor rupture of 
nonporous materials. Thus the limiting surface for either must 
be open at that end 

For an isotropic material having unequal tensile and com- 
pressive strengths, a = 6 = c in Eq. (12), and the axis of sym- 
metry of the surfaces which it describes becomes the hydrostatic 
axis. Considering the possible shapes of the surface, one may 
eliminate at once the values g = —1 (since then the simple tensile 
and compressive strengths would be equal); g = 2 and —1 < 
q < 2, since an ellipsoid is closed at both ends of the hydrostatic 
The hyperboloid of one sheet is obviously not applicable 
As for the 


axis 
to rupture, but might possibly describe vielding. 
hyperboloid of two sheets, clearly both sheets could not have 
physical significance. That sheet which is open toward the com- 
pression end of the hydrostatic axis, however, does have a plaus- 
3(1 + qg) so that a pure hydrostatic 
stress must be tensile to cause failure. Both g > 2 andg < —1 
are particular cases of the theory of strength first proposed by 
Nadai* * has shown can be stated as 


ible shape, provided a > o 


which 
is the octahedral shear and p is the hydro- 
static or octahedral differs from the Mohr 
theory in stress space only in that it is a surface of revolution, 
while Mohr limiting surface has a hexagonal cross section which 
can be inscribed inside the circular cross section of the Schleicher- 
Nadai surface at the same octahedral normal stress® ® 7, Both 
theories impose a higher degree of symmetry on the limiting sur- 
face than does the requirement of isotropy, and furthermore, 
careful triaxial tests of marble and other rocks by various in- 
vestigators, notably von Karman, B6éker, and Ros and Eichinger, 
vield results which cannot be described by a surface of revolu- 


Schleicher? and 


t, = f(p), where r, 


7] 


normal stress. It 


Expressions which do describe the results have been 


suggested.® * 
circular cylinder at 
® Such a surface could not be a hyperboloid. On 


tion.® § 
8 Further, the yield surface approaches asymp- 
totically a sufficiently high hydrostatic 
pressures.® 
the other hand, the rupture of a material should be described by 
a strain criterion rather than a stress criterion.” § 

Owing perhaps to a lack of data due to the difficulty of tri- 
axial testing of strong materials, metals having unequal tensile 
and compressive strengths do not appear to have been similarly 
investigated. Because brittleness is not so much a property of a 
material as a function of hydrostatic stress (among other things),’ 
and also because such metals—e.g., cast metals—tend to be 
brittle, it is reasonable to suppose in the absence of evidence to the 
contrary, that they behave much as the rocks do 

Having thus disposed of Marin’s hypothesis, a theory of 
strength for orthotropic materials proposed by Norris'® may be 
It also was generalized from the distortion energy theory 


The resulting equations 


noted 
but with the aid of a physical model. 
unlike Marin’s, could only be plotted with six-dimensional co- 
ordinates. Norris found his theory to be in fair agreement with 
the results of some biaxial tests; nevertheless, some similar ob- 


jections may be raised to it. 
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Me 


writer is not valid. 
would be necessary to make a comparison between test results 
This has not been done since there is 


that the theory proposed by the 


In order to prove his conclusion, it 


TOPPING considers 


and the proposed theory 
little available data on strength of nonisotropic materials sub- 
The theory proposed may only be 
The fact that the pro- 


jected to combined stresses. 
satisfactory for biaxial states of stress. 
posed nonisotropic theory reduces to a strength relation which 
has been found satisfactory for isotropic materials under biaxial 
stresses indicates that the proposed relation may provide an 
adequate representation for the biaxial strengths of noniso- 
tropic materials 

It is recognized that yield and ultimate strengths correspond to 
different phenomena. It does not follow, however, that similar 
mathematical expressions (relating the principal stresses and 
uniaxial strength values) cannot be used to define both vield and 


ultimate strengths under combined stresses. 


A Finite Vortex Method for Slender Wing-Body 
Combinations 


George S. Campbell 
Research Engineer, Hughes Aircraft Company, Culver City, Calif 


July 12, 1957 
INTRODUCTION 


T THE APPLICATION of the slender-body theory to aerodynamic 
problems, the base-plane geometry is frequently too com 


a 


plicated to permit the usual transformation of the body t 
circle. However, it has been found that such problems may be 
solved by using a reasonable number of two-dimensional vor 
tices and satisfying boundary conditions at appropriate ‘‘con 
trol points”’ The method is demon 
strated for a slender wing plus infinite body 
with the more exact theory.! Although this procedure has 


been used to predict loads on a missile with folding wing panels, 


in the cross-flow plane 
and compared 


details of this application are deferred to a later date 
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Fic. 1. Geometry of the wing-body problem 


rreichische 
Combine \NALYSIS 


The method for calculating the lift distribution of an arbitrary 
wing body is illustrated for the specific configuration shown in 
Fig. | \ slender wing is mounted on an infinitely long circular 
ylinder, the combination being at angle of attack a. The wing 
is represented by .V vortices across each panel, and the body is 
presented by a doublet at the origin. For each wing vortex 
there must be a corresponding image vortex within the body 

order to cancel induced velocities at the body surface The 
boundary condition of no flow through the wing is satisfied at 
N control points located midway between the vortices 


This procedure leads to a set of simultaneous equations that 


must be solved for the unknown vortex strengths > 
\ 

by the rz Cixi = 1 + (a/y,)?, =1,...N 1 
ision, it | 
results The coefficients C;, represent the downwash induced at the con 
there is trol points y; by the vortices at +y; and at the image points 
Is sub- +q2/y 
my be 
, l l l 
he pro- ( = me ds 

° » a , 4 , 
which - y J y y 
biaxial 
: l l 2 
ide an - » 

vy; — (a?/y y; + (a?/y 


10Oniso- 
The terms on the right side of Eq. (1) represent the potential of 


ond to the flow at infinity and the doublet at the origin \ll lengths in 

imilar these equations are expressed in semispans, and the vortex 

Ss and strength 4; is made nondimensional by +Vas. The positions of 
d and vortices and control points are, respectively, 

yi =at l1—a)/N | 2 

» 


For a given choice of V, the unknown vortex strengths 7; are 
btained by solving Eqs. (1) to (3 The wing loading is pro 


portional to the bound vorticity T 


yy 
r=); (4 
a 
alif The load distribution over the body may be obtained by evalu 
iting the velocity potential at the surface of the body due to all 
the vortices 
l vs dL dy y = 
\ 
Lmic 2 +» yi jw — 2tan 2y; ~/a* — y?)/(yi? — a*)]} 5 
om im 
Od 
be RESULTS 
sie In order to evaluate the accuracy obtainable with the approxi 
‘on mate method, results calculated from the preceding equations 
ss are compared with the ‘“‘exact’’ slender-body results of reference 
red 1. It is of course realized that the approximate method 1s 
has necessary only for more complicated cross sections 
els, The variation of wing and body lift with diameter to span 


ratio is shown in Fig. 2 for several values of V, the number of 
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vortices Lift has been divided by lift of the « xposed wing alone 


to give the usual interference factors of Morikawa.? The calcu 
lations show the expected improvement in accuracy with an in 
creased number of vortices V = corresponds to the exact 
theory.?) Calculations for a wing alone show that the error in 


lift prediction is inversely proportional to 


The load distribution for a wing-body combination of 0.3 


diameter to span ratio is shown in Fig. 3 rhe exact result was 
obtained from the potential function given in reference 1 so that 


for unit semispan 
1/gas) (dL/dy)\y = 4Y (1 — ¥? lL — (a*/#* / 


1/gas) (dL/dy)y = 4 (1 + a2)? — 4y? — 8a? — y?) 


The use of a finite number of vortices overpredicts the tip load 
ing of the wing but has relatively little effect on the body load 
distribution 


By way of summary, it has been found that a reasonably small 
number of two-dimensional vortices can be used to calculate the 
aerodynamic loads on slender wing-body combinations rhe 


1 


method is sufficiently general to permit arbitrary missile cross 
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sections provided computing machinery is available for the 


numerical calculation 
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Calculation of Local Heat-Transfer Coefficients 
on Slender Surfaces of Revolution by the 
Mangler Transformation 


Shao-Yen Ko 

Senior Research Scientist, Minneapolis-Honeywell Regulator 
Company, Hopkins, Minn 

August 16, 1957 


INTRODUCTION 


7 HEAT TRANSFER by convection from slender surfaces of 
revolution in axisymmetrical flow is of current interest with 
regard to projectiles, airplane fuselages, rocket motors, ete. At 
sufficiently high altitudes and speeds laminar boundary layers 
can exist at very high Reynolds Numbers. For a laminar bound 
ary layer, the velocity and temperature profiles, momentum and 
displacement thicknesses, and wall shear stress can be evaluated 
by using the Mangler transformation This enables direct use 
of the solution for flows over a two-dimensional body with a free 
stream velocity distribution for the determination of the bound 
ary-layer quantities in the three-dimensional case. This paper 
extends the Mangler transformation into the evaluation of the 
local heat-transfer coefficient on slender three-dimensional sur- 
faces of revolution in axisymmetrical flow. The mean heat 
transfer coefficient over the entire surface can be evaluated by 
integration. A typical example of this transformation technique 
is given for the case of laminar heat transfer from semi-ellipsoidal 


surfaces of fineness ratio 4 


MANGLER TRANSFORMATION 


The three-dimensional body with a given pressure gradient was 
transformed to a related two-dimensional body by means of the 
following transformation variables 


e =f 1K* 1.2(0)/L41 dx () 


y = [Kr,(x)/L] y (2) 


Where x was measured from the forward stagnation point 
along a meridian profile, y was measured from the surface of the 
body along a normal; 7,, was the radial coordinate of the body 
of revolution in a meridional plane and L was its total length. 
The barred variables are coordinates for the equivalent two- 
dimensional body. A is an arbitrary constant 


Furthermore, it was assumed that 


and T(x, y) = T(s, 9) (4 


CALCULATION OF LOCAL COEFFICIENTS OF HEAT TRANSFER 


The theoretical calculations predicted the distribution of the 
local coefficients of heat transfer assuming the constant proper- 
ties and incompressible flow. The method utilized the Mangler 
transformation as a basis. This enabled using directly the 
wedge-type flow solutions for the corresponding two-dimensional 
surface in the determination of the boundary-layer quantities for 
the three-dimensional case. 

The pressure distribution used in this example calculation for 
flow past a semi-ellipsoidal body with axis ratio of 1:4 was the 


theoretical approximation given by McNown and Hsu.’ The 
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Local coefficients of heat transfer from semi-ellipsoid 
surfaces of fineness ratio 4 for Prandtl Number of 0.7 


1 


corresponding velocity distribution at the outer edge of 
boundary came from the pressure distribution by means of 


Bernoulliequation 
U* = U/U, = 1 — [(p — p Uo?/2 


In the evaluation of local heat-transfer coefficients on t 
surface of the semi-ellipsoidal body, 25 calculating stations wer 
chosen along its major axis with finer steps near the stagnatior 
point. The dimensionless values of ** for the two-dimensional 
case were obtained from Eq. (1) by graphic integration at eacl 
calculating station. The Euler Number of the equivalent two 
dimensional surface was evaluated at each calculating statior 
from the definition 


m = (#*/0*) (d0*/dx* 6 

For each Euler Number m, there is a corresponding wedge. The 

graphic results by Brown! gave the heat-transfer parameter 

f(m) for wedge flow, here applied as 

(hs R)/N (OS vy) [ry t/S) Vv (0/0 = f(m 7 

where S was the total length of the body of revolution measure: 
from stagnation point along the meridian profile 

The last step of the calculations was to carry out the inverse 

transformation to get the heat-transfer parameter for the three 

dimensional case from the values for the two-dimensional case 

The following relations between the two-dimensional quantities 


and the three-dimensional quantities exist: 


For Two-Dimensional Case For Three-Dimensional Case 
§ = S 
h = h/Kr 
0 = l g 
U = l 


Z = [°K #2 dy 


Substitution of Eq. (8) into Eq. (7) yielded the local heat-transfer 
parameters at each calculating station on the surface of the 
semi-ellipsoidal body as 


S/L 


Nx hS/k 


Nr V (Uo5)/» #2 dx*) /p 


py. . 

- \ : =a 

VJ. 
Fig. 1 presents the final results of local heat-transfer parameters 
on the surface of the semi-ellipsoidal body Local heat-transfer 
parameters of a flat plate based on Eq. (10) also appear in this 
figure 

hS/k UoS/v)2? = 0.292/(x/S)'/? 10) 


1 


Integrating the theoretical local coefficients determines the mea 
coefficients 
A mean coefficient of heat transfer around a body of revolu 


tion is by definition 


The result of this integration is expressed as 


The 
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e value of Cw found to be 0.672 his value has In addition it is necessary to imsurt t ! ern 
- rt ( rom Eqs. (3) and (4) be identi I or é 
red with the m transl periment ] result b uined I Eq : . i 
} ind (5) the difference between these radial di ements 1 
Ki id nl rhe presented calculation gave re ss ‘ 
ver cent lower than their experimental values In view \ ij 7 
t 1 \ nature of this leulat gree ent ; 
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rpose of this note 1s to determine the thermal stre nan un 
d circular bulkhead (Fig. 1) cor ting of a web and a cay embly 
1 to a general axially symmetrical radial temperature distribution 
frequently occurring case of a parabolic temperature distribution, which 
linearly varying and uniform temperature considered in detail 
il case of a centrally located “hot spot ncluded in the general 
In actual practice the bulkhead is attached to the outer shell of 
we of the aircraft and additional stresses are induced WEB-ME mBERC) CAP -MEMBER(Q) 
. > = 
(GENERAL SOLUTION | 7 
o« i) 
2 
; DIFFERENTIAL EQUATION? of a disc subjected to an axially - 42 oe So 
symmetrical temperature distribution is : re . 
- 
. " 4 
V'¢ —akV?7 l ® 
; ; Fic. 1 Structure id temp ture distribution 
here @, a, E are respectively the Airy stress function, the co- 
ficient of linear expansion, and the flexural modulus 
Phe components of stress are determined from 
0 1 /r)(d@ dr 9 on = d*@ dr?), T 1p ) 2» 
I tress-strain relations are 1000 F— 
€ du/d 1/E)\io, — vo al ) 
ot(ib/in) 
€ 4/7 l E o — vo t al + 
is the radial displacement 500fF- 
It can be shown that 
1 log + Brrlogr +er?+ D— ak r dr/? T rd? 
5 
i ; : | 1 r 
s the complete solution of Eq. (1 Since the stresses and dis- 
placements are unaffected by D, it is taken to be zero .: ieee — F T 
Considering the case of a web without a central hole, the con 
' : ; : ; 
iti of finite stress at 7 0, requires | B=zQ@ Thus the , , Se i + = — 
tress function becomes i 
dy = ar? — aE; J (dr/r) J Trd? 65 -12900¢—__+__—_—+ "+ : 4. : 4 
1@] 15 30 45 
i » be determined from the associated stress or displacement r (inches) 
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u(a) = Ue(a) (15 
where the last condition can be replaced by 


€9,(a = €9,\a 16 


2c; (= E, t ) T*(a 
aEeT*(b) \arEs ti) T*(b 
te Bs 
Ge ar) ee eee se 
2 t, Ey 


T*(r) = (1/r?) fT rd 20) 


where 


NUMERICAL ILLUSTRATION 


For the purposes of illustration consider the structure of Fig. 1. 
The web and cap are assumed to be made of titanium and stain 
The 
properties, and temperature distribution are indicated in Fig. 1 


less steel respectively given geometrical and physical 


The resulting stress distribution is shown in Fig. 2 
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Boundary-Layer Correction in Supersonic 
Nozzle Scaling 


Abraham Kogan 


Department of Aeronautical Engineering, Israel Institute of 
Technology, Haifa, Israel 


September 3, 1957 


SYMBOLS 


= length of supersonic part of nozzle 
Vf Mach Number 


\ positive integer 

“u = velocity 

x = length parameter 

y = distance from wall 

56 = boundary-layer thickness 


displacement thickness 


@ = momentum thickness 

L = kinematic viscosity 

p = density 

to = local friction stress at the wall 
1 = value outside boundary layer 


= value corresponding to arithmetic mean temperature between wall 
and free stream 


DISCUSSION 
i — VON KARMAN’S Momentum integral for compressible 
flow, Tucker!) ? has obtained the following results for the 
boundary-layer development in a supersoniy tunnel: 


[6(x)//] = (1/F)} [6(0)/2] 4+ 
ex/l 
{ F(5/6) (p/p1) (ro/puy?)d(x/1 (1) 
vW 
and 
5*(x)/1] = (1/G) } [6*(0)/l] 4+ 


ey /| ) 
J G(6*/6) (p/ pr) (70 /puy*)d(x 1), (2) 
0 


* 


The ratios 6/6 and 6 


6 were calculated by Tucker as functions 
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(16), and (9) to (14 


coefficients are determined to be 


Using Eqs. (4), inclusive, the values of t 


A» beak. T*(b =|-—- 2 asEs7 *(b lj 


~ ™ 
| 
a, 
a 
~ > 
Bs | BS 


of the local Mach Number .V/ by postulating the power law 


They are tabulated in reference 2 as functions of .\/ for different 


values of LV. 


The parameters F, G, and p/p; in Eqs. (1) and (2) are also func 
tions of 17 only: 
, 7 Pree 
. 0/6 2 _ M? T 6*/6 
F= 3 exp ; , di 
(0/6 )yp—1 J 1 Mj1 t (¥ 1)/2] A1?} 
G = |(6/6*)/(6/6 ul i 
p/p = 1/11 + f(y — 1)/4). 


Only the skin friction coefficient 7)/(1/2)%,? varies with th 


Reynolds Number based on the distance x from the throat. W 
that the 


PF 1/7 
2)pu.? = 0.0262 
UX 


If we neglect the variation of 7 


observe, however variation of the friction coefficient 


with x is very slow 
7 /(1 


(1/2)pu,? with x and replace 
it by a suitable constant, the resulting approximations for Eqs 
(1) and (2) will contain integrands that depend on x only im 
plictly through 17. 

Let us consider a family of nozzles obtained by proportiona 
sealing of a given nozzle. The length parameter x// will be the 
MW for the whole We 


represent the boundary-layer development for any one of these 


same function of family can therefore 


» 


nozzles by the following approximate form of Eqs. (1) and (2 


fu uf / 
a(x)/I] = [1/FAD) 4 [a(0)/1) 4 J (Md Mt 
\ | } 
(6*/1) = [1/G(M , 6*(0 l| + J g( Wd My 
with 
f(.\1) = CF(6/6) (p/p d(x/l)/dM 
g(.\7) = CG(6*/0) (p/p) [d(x/1)/dM 


where C is the constant approximate value of the friction coeffi 
cient 

Thus, it appears that when a nozzle is developed from a givet 
prototype by geometrical scaling of the perfect-fluid-flow wall 
streamline, the boundary-layer correction may be calculated to a 
fair approximation by proportion from the prototype boundary- 


layer correction. 
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lues of th Hypersonic Flow Around a Sphere If v is assumed to be of the form Ll’. f(@)¢ , then according to 
Eq. (4) f = sin @ so 
R D nnell 1 
e staff Scientist, Convair Scientific Research Laboratory v= U. sin % 7 
San Diego, Calif F E 5) itis 1a 
: ‘rom Eq. (6) it is found that 
September 4, 1957 eX ere 
P “t= 2 COs O¢ : Ss) 
SUMMARY 
For g assume a three term series in 4 
‘ This note presents an estimate of the stagnation-point surface-velocity 
‘ 
gradient and the axial shock detachment distance for hypersonic flow around g =a t+ bb + C8? Q 
i sphere Real gas effects are taken into account rhe predicted values 
correlate well with available experimental data which are for perfect gas The constant a is zero since the radial velocity must be zero at the 
nditions surface of the sphere According to Eqs. (4) and (7 must be 
] e ] 
I! given by 
SYMBOLS 
b = (D/2) V2nk(1 — k) + R? 10 
constants, Eq. (9 
D diameter of the sphere Additional boundary conditions can be obtained at the shock, 
function of 4 on the axis where 6 = A From kl ind Eq. (8 
function of 7 
poof 2bAg + 2cAo? = k [Ap + (D/2)]? il 
pressure . 
Ano . ) id « O ned by matcl y 
radial coordinate, distance from center of sphere Another relation between Ay and c can be obtained | matching 
aw R shock radius, radial coordinate of shock 
radial component of velocity 
l free-stream velocity 1.0; T — 
circumferential component of velocity ; . 
different 3 2 D)v'g _¢, stagnation-point surface velocity gradient ho ! ' Ps os 
D 2), radial distance from sphere surface | Yl -* 
lso func ratio of specific heats [| > ii | 
eer 8 | 
mass density 
polar coordinate angle, measured from stagnation point , 690 — PERFECT GAS Xx 
1M Subscripts | 
conditions in free stream 6 4 7 | eines: 
conditions just behind the shock 6D 
stagnation point conditions 2U REAL GAS r 
0 conditions at the axis, @ 0 bad 
primes denote differentiation 4 | | 
. | | ” 
ith the | | 
t. We THE VELocIry GRADIENT 
ae : : © REF. 2 DATA 
: 7 DENSITY in the flow between the shock and the sphere is 
assumed to be constant at p;. This permits use of the in- 2 | | —— 
compressible Bernoulli equation along the axis | | 
p =p — (piu? 2) 1 | | 
replace Since the shock is normal at the axis, the normal shock equations @) mi | —— = 
x Eqs give u, = —kU., and p, = po +(1 — k) po U .? directly behind 0 2 at 6 8 1.0 
ly im the shock. These values are substituted into Eq. (1) to obtain 
the stagnation pressure k 
rtiona : , 2 . . ‘ ‘ 
b, = p + 2, U2 [21 — Rk) + RB) /20} (2 Fic. 1 Stagnation-point surface-velocity gradient parameter. 
be th } / vs ‘ ~ g I Ys I 
> Si 
retort If the surface pressure coefficient distribution is assumed to be 
these Newtonian, giving 1.4 . 
»)- * | ] | 
ee - PERFECT GAS | 
then the surface velocity obtained from the Bernoulli equation 
on the surface is 1.2 
= [ V 2AkK(1 — k) + Rk? sin 6 (4) r= | 
ind @ is 1.0 | 
B = (2/D)U. V2kM1 — k) + k? (9) kD. 
Pas laa , v REAL GAS 
cue: If \ is set equal to unity (p, = p,) in Eq. (5), the Li and Geiger! j O 
result is obtained. Values calculated from Eq. (5) are shown in 8 = : 
giver Fig. 1 along with test data points? for perfect gas conditions in a ‘e 
scatl wind tunnel. The correlation for perfect gas is improved by us- 8 O REF 4 DATA 
re ing perfect gas values of A for 7 = 1.4 
da ° . . . 
— The effect of using real gas values of A, obtained from Feldman,‘ cm —— | _| 
ary- : : es 
is to decrease 8 about 5 per cent ee 
THE DETACHMENT DISTANCE 
| | 














Z : > 1@) 1 — 
t in The shock detachment distance corresponding to the surface Oo 2 4 6 8 10 
ation conditions obtained above is obtained by starting with the con- ; ' : : : 


tinuity equation . 


sin 6 (Our?/Or) + 7 (Ov sin 6/08) = 0 (6 Fic. 2. Shock detachment distance parameter 
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dv/d0é tlong the shock evaluated at the axis The oblique shock 
relations provide 
1v/d0 l l l k) (Ro"/R 12 
where A,” is unknown. But for a sphere the shock is almost 
pherical near @ = O and therefore Ry” is assumed to be zero 
A fourth condition is then obtained by equating the shock rel 
tion value to the value obtained from Eq. (7 
Zod \ D/2 i) 
limination between Eq ind >) give he equation for 
I Iq 1] 13) g t juat f 
the shock d chment distance at the axi 
2h Pi) 2h/D tk(2b/D l 2b/D 2h 
14 
D 2(1 
where D is given by eq 10 
Caleulated values of the detachment parameter 2Ao/RD ir 
shown in Fig. 2 Phe correlation with the test data? obtained 
under perlect gas conditions is good for both 6 based on A l 
ind for perfect gas A Phe effect of using real gas values of A is 


to increase Ay about 2 per cent Phe theoretical result of Li 


ind Geiger! is also shown 
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Creep Buckling of Tubes in Torsion 


lain Finnie 
Materials Engineering and Corrosion Department, She 
Development Company, Emeryville, Calif 


Cor i 40 7 
September 5, 1957 


¢ bm RE Is little published information on the creep buckling of 
members other than pin ended columns. Even in this rela 
complete creep-buckling 
In the 


tively simple case the complexity of a 


analysis has led to the use of approximate solutions sim 


plest of these the effects of column eccentricity are ignored, and 
the tangent modulus of the inelastic buckling solution is replaced 
by a time-dependent tangent modulus The object of the present 
note is to indicate the extension of the time-dependent modulus 
approach to other types of buckling and, in particular, to com 
pare the prediction for torsional buckling with experiment 
Inelastic buckling solutions may be obtained by replacing the 
elastic modulus in the elastic buckling equations by an effective 
the 


modulus involves both 


the 


general, the effective 
E, = 0a/de 


stress and strain.t 


modulus In 


tangent modulus and secant modulus EF, = 


o/e where o and e€ are It is known that when 
longitudinal bending predominates, as in the familiar Euler col 
the When 


secant 


modulus 
the 


umn, the effective modulus is tangent 


longitudinal bending is small compared to twisting, 


} 


modulus is the effective modulus. Examples approaching this 


condition are a long evlinder in torsion, a cruciform section 


column which buckles by twisting, or a long flange with one un 
loaded edge simply supported. An example, given by Stowell, 
tf Under multiaxial stre o and e must be replaced by the invariant quan 


titie 


written here in terms of principal stresses and strains 
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of the more general case is the square plate loaded as a col 
for which the effective modulus is (0.114 a ¢ (L886 00/0 

When only the tangent modulus ts involved, we m vrit 7 
critical stress for inelastic buckling as @ C E,; where C is - 
function of the specimen dimension Writing the tot t ~ 
undet creep conditions of the column terial é c 
where ¢ is time and £ the elastic modulus, we find o ( F 
Of/Oe0) from which o,, is obtained funct f time F 
ipplication of this type of solution to pin ended colut 4 
discussed by Carlson c 

Following the same procedure ily t cant dt : 
involved, we write o C E, or Candt tr disap] C 
from the buckling equation Phis implies that for buckling : 
dominantly by twisting there is a critical strain indepen 
the ereep bel vior of the teria Phe time twiik It crit 
train is reached will, of course, depend on the stre ind mat 
properties. For the more general case in which bot ecant 
tangent modulus are involved, we would replace both t 
quantities by their time dependent lues iI 

A prediction of torsional buckling will illustrate the preceding 
nethod For hollow evlinders, length 1, wall thickness h, 
radius Rk, such that 50 h/R L/k 10 R/h, Batdort 
given the critical shear stress for elastic buckling 

T 0.747 E(h/R RIL | ia 
Gerard? has stated that the inelastic buckling stress may be e 
tained from Eq. (1) if & is replaced by E.(1 ue | ue é 
Taking w, and yp, the elastic and plastic values of Poisson’s rat . 
is 0.3 and 0.5, Eq 1) becomes Ma 
r O.865 Eth/R R/I - fail 

The secant modulus &, = o*/e* may be written in terms of t tu 
shear stress 7 and shear strain Remembering that for tors \ 
of a tube, the principal stresses and strains are 
Tr, € = e = 2, and using the invariants defined in the p S 


ceding footnote, we obtain EF, = 37 


leads to } 


E¢ 
for the critical shear strain at which torsional buckling shou 
occul 

Since unwanted buckling had troubled the writer during 
studies of multiaxial stress rupture, several tests were run t 
investigate the applicability of Eq. (8 Just as with ¢ 


] 


umns, the l 


creep buckling of tubes in torsion is observed as 


deformation increasing with time Phere is, therefore, son 


difficulty in defining a time or strain at which it can be said that 


One indication of torsional buckling 
Meas 


twist are 


creep buckling has occurred 


is the decrease in axial length of the tube being twisted 


urements of this quantity as a function of the angle of 
for 


three different twisting moments 


shown in Fig. 1 a plastic tube and for aluminum tubes under 


The angle of twist-time curves 


aluminum tubes are shown in Fig. 2 


he total 


for the 


plastic tube t time of test, 0.2 hour, was too short t 











enable this relation to be obtained. The specimen dimensions Oo 
are L = 5in., h = 0.100 in., R 0.75 in., and the corre 
sponding critical angle of twist 6 L/R predicted by Eq 
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Fic. 1 Axial shortening as a function of angle of twist for on 


plastic tube and three aluminum tubes 
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tected as a slight ellipticity of cross section (perpendicular D(x, y, 
é € é 
liameters of 1.600 in. differing by 0.010 in With increasing 
yt gle of twist the ellipticity and axial shortening increas¢ y YA € € et 
; As another example we may take the specimen dimensions The stream functions in a three-dimensional steadv flow were 
923 _ 25 » = YH j » } } } 
’s rati L 2.3 in., H = 0.03 in., Raverage = V.aOS in.) used by Johnson, used by some authors in the past Phe more proper name in 
Mathur, and Henderson’ for multiaxial rupture studies. Eq it nonsteady flow may be the particle function 
4 predicts > = ().058 while for a copper specimen which 
< failed in rupture after 5,442 hours these authors reported a rup- REFERENCI 
sof the ture shear strain of 7 = 0.144 Skobelkin, V. I.. Veriat p H at es 
torsior We would conclude that the general form of the curves in Figs. 1 ics JETP, Vol. 4, No. 1, pp 68-72, February, 1957 
ind 2 lends support to the hypothesis that torsional buckling 
the pre should occur at shear strains depending only on the specimen di 
Eg. (2 mensions and independent of stress or the material tested 
However, the critical strain predicted by the secant modulus 
underestimate the shear strain at which creep buckling Mass Transfer Cooling at Mach Number 4.87 
becomes significant 
shoul c ‘ amt 
ipl ; B. M. Leadon,* C. J. Scott, and G. E. Anderso 
) Pe 
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" Rosemount Aeronautical Laboratories, University of Minnesota 
during Carlson. R. L.. Time Dependent Tangent Modulus Applied to Colum» Creet Minneapolis, Minn 
run t B g, J. Appl. Mech., Vol. 23, p. 390, 1956 September 9, 1957 
th col Stowell, E. Z., A Unified Theory of the Plastic Buckling of Columns and 
1 ] NACA Technical Report No. 898, 1948 7 - 
renee Batdorf, S. B., A Simplified Method of Elastic Stability Analysi M* TRANSFER EXPERIMENTS on a 9 mil Wire porous cone 
some ylindrical Shells (1), Donnell’s ns, NACA TN 1341, 1947 of 20° total angle have been conducted at ./ +S using 
1 that . i eal biis f This led art ‘ — ‘ : : : 5 
ld th Gerard, G., ( snc teleseiaste ling of Thin Walled Cylin air and helium injection. Details of the experimental technique 
icklin Yield Region, NACA TN 3726, 1956 , ailtieds ae ak o : : , 
ng ‘ : ‘ S nees 1 and 2 In tl minar bo { 
rome ER a OR ROR Eee Me pee CE PTET eee ire described in refere ices 1 and le laminar undary 
2 eas f f Commercially Pure Copper at 250°C., The Engineer, Vol. 196 laver the recovery factors and heat-transfer coefficients measured 
ISt are pp. 261-265, 299-301, 1956 with zero injection agreed within +6.5 per cent with theory 
ley ah . , 1 " 
aes lransition Reynolds Numbers observed on the porous cone with 
curves i zero injection were half as large as observed on a smooth, imper 
; ; ; : ; ‘ ‘ 
of the meable model of identical geometry in the same channel, but 
wiih injection of large amounts of air or helium did not cause transi 


nsions On the Stream Functions in Nonsteady Three- tion to move forward from its zero-injection position on the 
— Dimensional Flow porous cone Distributed roughness of this type apparently does 





y Eq not disturb impermeable wall theory, but it masks whatever 
M. Z. v. Krzywoblock effective roughness may be caused by discrete pore injection 
] Professor, Panel on Fluid Dynamics, University of IIlinois The 12-in. model had a 2%/,-in. solid tip, 18/4 in. of increasing 
Irbana, Il porosity, and then 3 in. of fairly constant porosity nd laminar 
September 9, 1957 boundary layer) to the beginning of the transition region Local 
measurements of temperature recovery Tactors and heat-transfer 
— R the continuity equation coefficients plotted against an injection parameter * (Fig. 1 
: A indicated that the effect of injecting air is well predicted by 
a ; ie ier 0) theory® even though the actual injection distribution differed 
considerably from the assumed (a 2) distributi 
Che stream functions satisfying identically this equation are With helium injection the temperature recovery factors dis 
1) One-dimensional nonsteady flow (t, x played a tendency to rise with f..* before falling below the zer 
0 pay oa) = —%,- b= diz. i injection value (see Fig. 2 Phis tendency is been found 
P'wo-dimensional steady flow (x, y + This research wa Office of Research of 
r one the ARDC, USAF, under ¢ 00 2 
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Effect of air injection upon laminar temperature re 
covery factors and local heat-transfer coefficients 


Fic. 1 


theoretically under the assumption of constant Prandtl Number 


(P,),! but not when P, variations were considered.6 Local heat- 


transfer coefficients with helium injection appear to confirm 


variable ?, theories although experimental recovery factors were 


used in reducing the data. For the range of / 


measured concentrations of helium at the wall agreed with 
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Fic. 2. Effect of helium injection upon laminar temperature 


recovery factors and local heat-transfer coefficients. 


in Fig. 2 the 


and x as two independent variables, Eqs. (1 
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Variation of enthalpy recovery factors with wall con 
centration 


Fic. 3 


variation of enthalpy recovery factors with wall helium mass 


fraction is shown in Fig. 3 
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Boundary-Layer Slip Solutions for a Flat Plate* 


Hidenori Hasimoto 

Assistant Professor, Department of Aeronautics, University of 
Kyoto, Kyoto, Japan 

September 10, 1957 


open FLAT-PLATE boundary-layer problem with slip has been 
studied by several authors. In this problem we must solve 


the boundary-layer equations 
u(Ou/Ox) + v(Ou/Oy) = v(02u/Oy? (1 
(Ou/Ox) + (Ov/Oy¥) = OU (2 
taking into account the boundary conditions 


uw = h(du/oy), v=0 aty =0 (3 


u—>l1 atx =Oandasy—> (4 


where ft is the slip parameter proportional to the mean free 
path length 

Recently Schaaf and Sherman! discussed the qualitative na 
ture of the flow by solving the Rayleigh equation (or Oseen equa 
tion) corresponding to Eq. (1). Bell? examined the solutions of 
this simplified equation and showed the nonuniformity of the 
local skin friction as x ~ 0 and h > 0 

In this note I want to show a direct method of approach t 
Eqs. (1) to (4), using the von Mises transformation 


/ 


Taking the stream function y, such that 
u = Oy/oy, v= —OY/Ox s) 
to (4) are trans 


* The author wishes to express his cordial thanks to Professors I. Imai and 


kK. Tamada for their continual encouragement throughout this work 
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| b= QV wx/h ie) 
| é « Vv) ; : 
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4 n= WV/\2V vn 10 
| 
| js. (6) to (8) are reduced to another form 
O0/On) |ulOu/On + 2n(Ou/On) = 2ZEOu/OE (11 
——————————_ 
Ou/On = E— atyn = 0 (12 
u=1asn— 13 


wall cor 
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Plate* 


of 


as been 


st solve 
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n free 


ve na 
equa 
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ich to 





\ssuming that w can be expressed in a power series of & 


“i= l , ae Ul 
n=1 
Ve obtain from Eqs. (11) to (13) 
= f 4 2nf of, = O (j) = f = () (15 
= + 2nf if, = f fif f."(0) = 0 
f =() (16 
n) is found to be 
fi(n) = —Ei(n 17 


vhere E,(m) is the repeated integral of the error functior 


Il 


. 
I Y = 1" erfc yn = | E, n dn, E n 
J 1 


erfe n = (2/Y x) | é dyn (18 


Introducing f; into the right-hand side of Eq. (16), and taking into 
wcount the boundary conditions, f2 can be easily determined to 


nt 


(1/2) (E,\2 + E,xE.) + (7/4)E (19 


f(n) = 
Making use of these results the local skin-friction coefficient 
C; and the total skin-friction coefficient Cy, for a plate of length L 


ire given respectively by 


C; = [pu(x, O)/h| /(pto?/2) = (2v/hu x 
;1 \2/~ x) (ox hV tao) 4 [(5 $+) (2/a x 
(yx /h?ua) 4 } %)) 
a] 
ind ( ») = C dx ] 
J 0 
(2v/huax) }1 (4/3V 2) UV eL/hV uc) 
[(5/8 (l/x |} (ol h?ua) + } (21) 
C, is finite at x = O, in contrast with the Blasius solution corre- 
sponding to the case A = 0, and is larger than the value expected 


from the free-molecule flow theory Further it may be noticed 


that the value of Cp, obtained from the Rayleigh equation is given 


+ (vl /2h?u + 


fi — (4/3V 2) (VoL /hV cw ) 
which is in accord with Eq. (21) as far as the second term and is 
larger than Eq. (21) because of the third term. In this respect 
the Rayleigh solution may be considered to offer a very good 
ipproximation to the solution of Eqs. (1) to(4) if YW »L/(hy a9 
is small, although in this region Eqs. (1) to (4) must be used with 


some reservation. 
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A Simple Graphical Method for Constructing 
Two-Dimensional Supersonic Flows by 
Means of a Drafting Machine 


F. Edward Ehlers 
Mathematical Services Unit 
Boeing Airplane Company 


Physical Research Staff 
Wash 


Seattle 


September 10, 1957 


SUMMARY 


4 simplification of the Puckett field method of characteristics for super 


sonic two-dimensional flows is explained By means of special tables given 
here, an engineer can construct the Mach net rapidly using a drafting ma 


chine without any intermediate calculations 
INTRODUCTION 
P' CKETT® HAS DEVELOPED a simple method of designing two 
dimensional supersonic nozzles by the method of character 
of short, 
| 


h junction 


straight 
Con 
is introducing a Prandtl-Meyer expansion 


istics. The walls of the nozzle are made up 
sections having the same angle increment at ea¢ 
sidering each vertex 
through a given angle and resulting in a deflection of the stream 
of this same angle, one can find the direction of the lines of the 
Mach net by simple calculations using the tables of the Prandtl 

Phe method 
Puckett!» * but 


\ 


Meyer expansion angle in references 1, 2, and 3 


described herein is essentially the method of 


e necessity of at 


special tables are provided which eliminate tl 


intermediate calculations. The following discussion explains 


the use of the special tables given in this not 


THEORY 


Suppose a parallel uniform stream of air is flowing between 


two parallel walls at a Mach Number of unity At the point 


where the wall turns through an angle, A@, an expansion will 


occur The flow will be accelerated and turned through an angk 
Aé in a centered fan region of characteristics or Mach lines. If 
A@ is small, say 1°, then the centered fan may be approximated 
by a single Mach line across which the velocity changes discon 
ingle changes 


tinuously For an angle change of 1°, the Mach 


from wp = 90° to w = 67.574° while the flow angle changes from 
0° to 1° (see page 504.211-1, reference 2 rhe direction with 
respect to the horizontal of the last Mach line of the fan ts 67.574 


1° = 66.574 To approximate the centered fan by a singlk 
» 


Mach line, we choose the Prandtl-Mever angle of 1 which is 


the average of the Prandtl-Mever angles before and after the far 
The Mach angle corresponding to this Prandtl-Meyer angle is 


The direction of the mean Mach line then is found 


B= 72.099 
to be 











Mach net construction. 


Fic. 1 
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TABLE | Angular directions of the Mach lines running down 
ward to the right for positive flow angles and upward to the 
right for negative flow angles 





) 


TABLE 2. Angular directions of the Mach lines running 
upward to the right for positive flow angles and downward to 
the right for negative flow angles 


Me 6 = 72.099 71.599 


7 6 = 71°36’ 


Where this Mach line reaches the lower boundary (see Fig. 1 
it is reflected as an expansion wave in such a way that the angle 
The flow behind 
the reflected wave has expanded through a total Prandtl-Mever 


of flow behind the line is parallel to the wall 


angle of w = 2° while the flow before the Mach line has expanded 
through a Prandtl-Mever angle of » = 1°. The mean Mach 
line is given by the Prandtl-Meyer angle of w» = 1} which is 


found from the tables of reference 2 to be wp = 64.450 The 
angle of the mean Mach line with the horizontal is 


O° + 1 
uw +0 = 64.450° +4 = 64.95 


USE OF THE TABLES 


Tables for pu dand pw 6 are given for steps of 1° in w and 6 


by Tables 1 and 2. In constructing the Mach net, it is wise to 
put in each mesh the value of the Prandtl-Mever angle and of 
the flow angle (see Fig. 1). If the flow direction is measured 
positive in the counter-clockwise direction, the Mach line run- 
ning upward and to the right or downstream direction is the uw + 
@ characteristic for @ > 0. Since the table of w + 6 contains only 
positive flow angles, the same Mach line for @ < 0 is given by the 


uw — Otable. Similarly, the angle of the Mach line downward 


TICAL SBCLIEBNC HS JAN 


cae, RVSS 


and tu the left is found in the table of u 6 for 6 > O and int 


table of u 6 for 6 8) For example, the angle of Mac li 
\B in the figure is found in the third row and first column 


Table 1 to be 59°27’, and the angle of the Mack line BC is four 


; ra ‘ = ert 
in the third row and second column of Table 2 to be 61°27 . 
the tables, the Prandtl-Mever angle is designated by w and t = 
ertal 
flow angle by @ : 
sels . referé 
rhe tables give results for flow angles up to 10 Higher fl ae 
ingles can be found by adding (or subtracting) multip] ian 
10° to the angles in the tables <i 
| 1 
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A Note on the Reflection of Sound Waves at an “~~ 
Oblique Shock; “a 


C. T. Chang’ 

Aeronautics Department, The Johns Hopkins University 
Baltimore, Md ind e 

September 12, 1957 


, iw ONLY DISTURBANCES that can interact with a shock fror 
downstream are sound waves. Following the general pr 
cedure of reference (1), it is shown in this note that the nature Vol 
the reflected sound waves and the reflection coefficient can b Ch 
readily determined. If 8 is the shock angle and @ is the incl wide 
nation of the normal to the wave front with respect to the dow: . 
stream main-flow velocity Ul’, and A the speed of sound of the 
downstream main flow, the wave front drifts along the shock at 


speed of (see Fig. 1 
ce A + U cos 6)/cos (@ — B 


Referring to the coordinate frame OX J)’, where the Y axis is 


taken along the downstream effective flow velocity U’,, the prob 


lem can be reduced to a steady one The magnitude and i 





clination @ of U, is given by the law of sines 


C./sin a = U,/sin B l’/sin (a — B 2 





\s indicated previously,! disturbances corresponding to soun 
waves can be conveniently represented by the scalar potentia 


¢%, governed by a simple wave equation, thus 
? o(X + ¥ cot Me + (YY — VY cot pu, 5) 
where te = are sm (A/U,) | T 


is the effective downstream Mach angle. Clearly, ¢; represents 
the given incident wave while @¢, represents the unknown re 
flected wave 

The component of the perturbed velocity taken along the } 
axis, c's and the perturbed pressure p are related to the scalar 
potential by 4- 


= oy, p = —odx/sin pe ) 


\t the shock, thev are related to the local shock angle vy bv the eu 
condition 


t The work reported was sponsored by Project Squid, which is jointly 
supported by the Office of Naval Research, The Office of Scientific Research 
Air Force), and the Office of Ordnance Research (Army) under contract 
numbers N6-ORI-105 T.O. III, NR-098-038 

* Formerly, Research Assistant, Aeronautics Department, The Johns 
Hopkins University Now at the Division of Engineering and Applied 


Physics, Harvard University VIG 
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READE 


[?| [? | [ Vs VU cos B — C,/A 
D, ~ LB m1 (.M cos 8 — C,/A 
at X sin (a — B) 4 


here flow parameters with subscript (7) are those pertaining to 
reflected waves while parameters with subscript (2 


as indicated in 


are those 


ertaining to incident waves. 73), m1, and 74, 


and (2), are parameters which depend on the 


defined as the ratio of unperturbed pressures 


references l 
trengh x, 
cross the shock 


1 
shock 


Using Eqs. (3), 5) and (6), it can be deduced that? 


p >; = 


Pr/} = K(6, B, x 7) 
The parameter K can now be viewed as a reflection coefficient 
For a given medium, A generally is a of the shock 
strength x, the shock angle 8, and the inclination @ of the incident 
Its variation with respect to @ at a fixed value of 8 is given 


function 


vave 

in Fig. 1 
\ positive value of A indicates the reflected wave is of the same 

kind as the incident wave, namely, an expansion wave reflects as 

n expansion wave; a negative value of A indicates the reflected 

vave is of the opposite kind as the incident wave, namely, an 

For 


two inclinations of the incident wave 


xpansion wave reflects as a compression wave a given 


shock strength x, there are 


which no reflection occurs; one of them corresponds to the 

mal incidence, the other one depends on the particular value 

f the shock strength x under consideration 

Once the reflected sound wave is known, the shock deflection 
can be determined. Finally, the accompanying vorticity 

ind entropy disturbances generated at the shock can be evalu 


ted from the remaining conditions at the shock 
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FORUM 


Application of Gortler’s Series Method to the 
Boundary-Layer Energy Equation 


E. M. Sparrow 
Lewis Flight Propulsion Laboratory 
September 20, 1957 


NACA, Cleveland 


INTRODUCTION 


IT A RECENT ARTICLE of considerable importance, Goértler! has 
presented a new series method for calculation of stead) 
flows 


purpose 


conditions 
is to sketch 


laminar incompressible boundary-laver under 
of variable free-stream velocity Our her« 


out the application of his method to the heat-transfer problem for 


variable free-stream velocity 


the 


a two-dimensional body with 


unbient temperature 


Such an 


Both the surface temperature 7), and 
7... are taken to be uniform, but, of course, 7, + 7 
extension is easily carried out provided that full use is made of 


immediately able to provide 


GOortler’s results, and, in fact, we are 
numerical results for the heat transfer from a symmetrical body 
with a stagnation point for Pr = 1 
Detailed consideration was given by Gortler to free-stream v« 
locities of the type 
U(x) = x} \ \ i | 
U(x) = x} + Syyon x 
t ) 
; 2 
Numerical treatment was carried out for the cases 7 Qand m 


1 for the first series, and for m = 1 for the second series For 
these situations, universal functions have been computed? which 
used in determining the details of the boundary-layer 


may be 


flow on two-dimensional bodies 

Further, Gortler considered the case of an infinite cylinder in 
vaw subject to variable free-stream velocities of the type noted 
above. He formulated the equations for determining the cross 
parallel to cylinder axis) and gives 


of Eq. (1 


flow velocity component IV’ 


universal functions? corresponding to U(x for m = 1. 


THE HEAT-TRANSFER PROBLEM 


We start with the energy equation for incompressible steady 


flow over a two-dimensional body and then introduce the follow 
ing variables (due to Gortler 
. 
\ / 
v= tia) | Saef 7010) « 
/ vf \ 
e, 
= L/i / l x dx ) 
e 0 
r \ ea / 
I l Lia \ 
” ¥ fo es 
rhe result of this substitution ts 
1/Pr)e, + (F 2 Fro 2° FLO } 
where Py is the Prandtl number and 0 = (7 1 1 1 
The appropriate boundary conditions are 
ace, 0 = (), O(E, n) > las — tb 
But, except for the factor 1/P . Eqs fa) and (4b) are identical 
to Gértler’s Eq. (83) for the cross-flow velocity component VW of 
the yvawed cylinder problem. So, the mathematical formulation 
and results for the velocity problem on the yawed cylinder be 
come applicable, with only simple modifications, to our two 
dimensional heat-transfer problem In particular, for Pr = 1, 


there is a one-to-one correspondence in the results for 9 and WV’ 
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RESULTS For Pr = | able there. But, they may be easily obtained by following Gég 


: , ; , ler’s procedure. Introducing the expansions 
\s has already been noted, there are available? universal func I 5 I 


tions for the W problem corresponding to a free-stream velocity of 
Eq. (1) with m = 1—-1.e., for 


U(x) = sox + 35xX° + Sox? +. 5) 
(Eq. (5) corresponds to the flow about a symmetric body with a into Eq. (4), we find the following differential equations for cog 


stagnation point.) For this free-stream velocity, we can easil) puting the temperature distribution corresponding to the fre 


. . re > 4 "y (9 
write down an expression for the local heat flux g at the surface stream velocity of Eq. (2). 
of our two-dimensional body. First introduce the principal 


function B(£) according to the definition 


Pa 
= 2(7'(x) U(x) dx / U(x) (6) 
0 / 


and then set B(£) = By + Bi— + Bet? 4 for the present case 


With this definition, the heat-transfer result is 


0,/,.(0) = 0 
= ().5705 + (0.06212 B,)é 4 
0.02105 By? + 0.06845 B2)& + (0.00826 6; 
0.03880 BB. + 0.06273 Bs)E* + 0.00364 8)! 
().02222 B,?B2 0.03632 B18 0.01769 Bo? + 
0.06150 Bye + (0.00175 BY — 0.01299 8,38. + 0.02050 m= 1. 
B,7B; + 0.01994 8,82? 0.03437 6184 0.03300 B28; 4 
0.06016 B;)& 


Then, the functions 09,2 are written in terms of universal function 
in precisely the same fashion as is done by GOrtler for the Fi, 
his Eq. (114). Tabulated values for the F,,2 are available fg 


EXTENSION To Pr # 1 
where Nu. = gx/(T Rk, Re; The series method of solution for Eq. (4) with Pr # 

formulated without additional effort. For the free-stream ve 
The relationship between £ and x is provided by Eq locity of Eq. (1), it is only necessary to take Gértler’s Eqs. (108 
simple illustration, suppose that U(x) = sox, (s) = % = 1 ; , 
Then, & = Re,/2 and By = 1, 8B; = Be... = 0. 


the heat-transfer result becomes: Nu,/Re,"? = 0.5705 


“ ieneetane to (109); first replace all W’’s (or w’s) by thetas, and then muk 
For this situation, tiply all second derivatives by 1/P) 

‘ sins For the velocity variation of Eq. (2), all that is required is 
For the velocity variation of Eq. | 1) with m = 0 (flow around replace @” in Eqs. (9a) and (9b) by 6"/P) 
bodies with a cuspidal point at leading edge), Gortler does not The universal functions for F(z, 7) which are compiled by Gért 


provide universal functions for the W problem. Hence, results 


ler are, of course, usable as input data in the solution of Eq. (4 


for the temperature problem cannot be written without consider for any Prandt] Number 


able further computation. The necessary differential equations 
are given as Eqs. (104) through (109) in Gértler’s paper. It is 
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